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[1] est le physicien allemand Carl Friedrich von Weizsäcker (-) qui l'a proposée en [2],[3] dans un article publié dans le Zeitschrift für Physik[4],[5]. Les physiciens Hans Bethe (-) et Robert Bacher (-) en ont simplifié l'expression en [2],[3],[6],[7]. Le physicien Eugene Wigner (-) l'a généralisée en [2],[3],[6],[8]. En , Niels Bohr (-) et John Wheeler (-) y ont introduit le terme δ {\displaystyle \delta } [9],[10].

Expression

B ( A , Z ) = a v A a s A 2 3 a c Z ( Z 1 ) A 1 3 a a ( A 2 Z ) 2 A ± a p A 1 2 {\displaystyle 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 {\displaystyle a_{v}A} est le terme de volume[11],[12],[13] ;
  • a s A 2 / 3 {\displaystyle a_{s}A^{2/3}} est le terme de surface[11],[14],[15] ;
  • a c Z ( Z 1 ) A 1 / 3 {\displaystyle a_{c}Z\left(Z-1\right)A^{-1/3}} est le terme coulombien[11],[16],[17] ;
  • a a ( A 2 Z ) 2 A 1 {\displaystyle a_{a}\left(A-2Z\right)^{2}A^{-1}} est le terme d'asymétrie[3],[17],[18],[19] ;
  • a p A 1 / 2 {\displaystyle a_{p}A^{-1/2}} est le terme de parité[17] ou d'appariement[3],[20].

Remarques :

Dans le terme coulombien, Z ( Z 1 ) {\displaystyle Z\left(Z-1\right)} est souvent remplacé par Z 2 {\displaystyle Z^{2}} [11],[13],[21] car Z ( Z 1 ) Z 2 {\displaystyle Z\left(Z-1\right)\approx Z^{2}} [22].

Dans certaines formules, le terme d'appariement varie en A 3 / 4 {\displaystyle A^{-3/4}} plutôt qu'en A 1 / 2 {\displaystyle A^{-1/2}} [23].

Les valeurs des constantes utilisées sont (en MeV) :

  • av = 15,56
  • as = 17,23
  • ac = 0,7
  • aa = 23,6
  • ap = 11,2
Courbe représentant l'énergie de liaison par nucléon en fonction du nombre de nucléons dans le noyau.

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 {\displaystyle \,B=E_{v}+E_{s}+E_{c}+E_{a}+E_{p}} .
Illustration des termes de la formule de Bethe-Weizsäcker dans le modèle de la goutte liquide du noyau atomique. Les protons sont représentés par les ronds rouges et les neutrons par les ronds blancs.
Illustration des termes de la formule de Bethe-Weizsäcker dans le modèle de la goutte liquide du noyau atomique. Les protons sont représentés par les ronds rouges et les neutrons par les ronds blancs.

Pour expliquer ces différents termes, il faut supposer que le noyau est sphérique, de rayon R 0 {\displaystyle R_{0}} . Et comme il est compact (son volume est proportionnel au nombre de nucléons A), alors R 0 {\displaystyle R_{0}} est proportionnel à A 1 / 3 {\displaystyle 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 {\displaystyle 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 {\displaystyle \,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[24] :

r = r 0 A 1 / 3 {\displaystyle \,r=r_{0}A^{1/3}}
Avec r 0 {\displaystyle \,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 π ( r 0 A 1 / 3 ) 2 = 4 π r 0 2 A 2 / 3 {\displaystyle \,S=4\pi (r_{0}A^{1/3})^{2}=4\pi r_{0}^{2}A^{2/3}}
D'où,
A 2 / 3 S {\displaystyle \,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 {\displaystyle 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 = 3 5 ( 1 4 π ϵ 0 ) Q 2 R {\displaystyle E={\frac {3}{5}}\left({\frac {1}{4\pi \epsilon _{0}}}\right){\frac {Q^{2}}{R}}}

Q {\displaystyle Q} est la charge totale, R {\displaystyle R} le rayon de la sphère. En identifiant Q {\displaystyle Q} à Z e {\displaystyle Ze} et en prenant le rayon proportionnel à A 1 / 3 {\displaystyle 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 {\displaystyle Z^{2}} devient Z ( Z 1 ) {\displaystyle Z(Z-1)} . La valeur de a c {\displaystyle a_{c}} peut être calculée approximativement en utilisant l'équation ci-dessus :

R = r 0 A 1 3 {\displaystyle R=r_{0}A^{\frac {1}{3}}}
Q = Z e {\displaystyle Q=Ze}
Z 2 = Z ( Z 1 ) {\displaystyle Z^{2}=Z(Z-1)}

E c = 3 5 ( 1 4 π ϵ 0 ) Q 2 R = 3 5 ( 1 4 π ϵ 0 ) ( Z e ) 2 ( r 0 A 1 3 ) = 3 e 2 Z 2 20 π ϵ 0 r 0 A 1 3 = 3 e 2 Z ( Z 1 ) 20 π ϵ 0 r 0 A 1 3 = a c Z ( Z 1 ) A 1 / 3 {\displaystyle 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 = 3 e 2 Z ( Z 1 ) 20 π ϵ 0 r 0 A 1 3 {\displaystyle E_{c}={\frac {3e^{2}Z(Z-1)}{20\pi \epsilon _{0}r_{0}A^{\frac {1}{3}}}}}

La constante a c {\displaystyle a_{c}} du terme de répulsion électrostatique est :

a c = 3 e 2 20 π ϵ 0 r 0 {\displaystyle a_{c}={\frac {3e^{2}}{20\pi \epsilon _{0}r_{0}}}}

Une autre valeur de a c {\displaystyle a_{c}} peut être obtenue en utilisant la constante de structure fine :

a c = 3 5 ( c α r 0 ) {\displaystyle a_{c}={\frac {3}{5}}\left({\frac {\hbar c\alpha }{r_{0}}}\right)}

α {\displaystyle \alpha } est la constante de structure fine, r 0 A 1 / 3 {\displaystyle r_{0}A^{1/3}} , le rayon du noyau avec r 0 {\displaystyle r_{0}} qui vaut approximativement 1,25 femtomètres. Cela donne une valeur théorique de a c {\displaystyle a_{c}} de 0,691 MeV ce qui est peu éloigné des valeurs mesurées.

a c = 0 , 691 M e V {\displaystyle 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 ( N Z ) 2 A = a a ( A 2 Z ) 2 A {\displaystyle 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 {\displaystyle 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 { + A 1 2 cas pair-pair (A pair, Z pair, N pair)  0 cas pair-impair (A impair, Z et N de parité différente) A 1 2 cas impair-impair (A pair, Z impair, N impair)  {\displaystyle 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 1 3 4 a a A 1 ) Z 2 + ( a c A 1 3 + 4 a a ) Z + ( a v a a ) A a s A 2 3 ± a p A 1 2 {\displaystyle 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 ) {\displaystyle B_{A}(Z)} . Ainsi en cherchant les valeurs de Z {\displaystyle Z} qui annulent la dérivée B A ( Z ) {\displaystyle B_{A}(Z)} par rapport à Z {\displaystyle Z} , on peut obtenir une équation donnant les noyaux de la vallée de stabilité.

B A ( Z ) Z = 0 Z = a c A 1 3 + 4 a a 2 ( a c A 1 3 + 4 a a A 1 ) 1 2 A 1 + A 2 / 3 a c 4 a a {\displaystyle {\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

Article détaillé : modèle de la goutte liquide.

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[25].

Modèle de Duflo–Zuker

Proposé initialement en 1995 par Jean Duflo et Andres Zuker[26] puis raffiné en 1999[27].

Notes et références

  1. Taillet, Villain et Febvre 2018, s.v.Weizsäcker (formule de), p. 785-786.
  2. a b et c Pfalzner 2006, B.4, p. 217.
  3. a b c d et e Taillet, Villain et Febvre 2018, s.v.Weizsäcker (formule de), p. 786, col. 1.
  4. Wapstra 1958, § 15, p. 26, n. 1.
  5. Weizsäcker 1935.
  6. a et b Wapstra 1958, § 15, p. 26.
  7. Bethe et Bacher 1936.
  8. Wigner 1937.
  9. Wapstra 1958, § 15, p. 26, n. 3.
  10. Bohr et Wheeler 1939.
  11. a b c et d Le Sech et Ngô 2014, chap. 2, § 2.9, p. 30.
  12. Mayet 2017, chap. 2, § 2.1, p. 41, 1).
  13. a et b Reuss 2003, chap. 2, § 2.2.4, p. 49.
  14. Mayet 2017, chap. 2, § 2.1, p. 42, 2).
  15. Reuss 2003, chap. 2, § 2.2.4, p. 49-50.
  16. Mayet 2017, chap. 2, § 2.1, p. 42, 3).
  17. a b et c Reuss 2003, chap. 2, § 2.2.4, p. 50.
  18. Le Sech et Ngô 2014, chap. 2, § 2.9, p. 30-31.
  19. Mayet 2017, chap. 2, § 2.1, p. 42-43, 4).
  20. Mayet 2017, chap. 2, § 2.1, p. 43, 5).
  21. Taillet, Villain et Febvre 2018, s.v.Weizsäcker (formule de), p. 785, col. 2.
  22. Khater 2009.
  23. Le Sech et Ngô 2014, chap. 2, § 2.9, p. 31.
  24. « Le noyau atomique », sur www.sciences-en-ligne.com (consulté le )
  25. Guy Royer, Contribution à la description macroscopique des phénomènes nucléaires de fusion, fission et fragmentation (thèse de doctorat en physique nucléaire), (lire en ligne [PDF]), p. 14
  26. (en) J. Duflo et A.P. Zuker, « Microscopic mass formulas », Physical Review C, vol. 52,‎ , R23 (DOI 10.1103/PhysRevC.52.R23)
  27. (en) J. Duflo et A.P. Zuker, « Microscopic mass formulas », Physical Review C, vol. 59,‎ , R2347 (DOI 10.1103/PhysRevC.59.R2347)

Voir aussi

Bibliographie

Document utilisé pour la rédaction de l’article : document utilisé comme source pour la rédaction de cet article.

Publications originales

  • [Weizsäcker 1935] (de) Carl Friedrich von Weizsäcker, « Zur Theorie der Kernmassen » [« Sur la théorie des masses nucléaires »], Zeitschrift für Physik, vol. 96, nos 7-8,‎ , p. 431-458 (OCLC 5653580744, DOI 10.1007/BF01337700, Bibcode 1935ZPhy...96..431W).
  • [Bethe et Bacher 1936] (en) Hans A. Bethe et Robert F. Bacher, « Nuclear physics : A. – Stationary states of nuclei » [« Physique nucléaire : A. – États stationnaires des noyaux »], Reviews of Modern Physics, vol. 8, no 2,‎ , p. 82-229 (OCLC 4643731788, DOI 10.1103/RevModPhys.8.82, Bibcode 1936RvMP....8...82B, lire en ligne [PDF]).
  • [Wigner 1937] (en) Eugene Wigner, « On the structure of nuclei beyond oxygen » [« Sur la structure des noyaux au-delà de l'oxygène »], Physical Review, vol. 51, no 11,‎ , p. 947-958 (OCLC 4644011967, DOI 10.1103/PhysRev.51.947, Bibcode 1937PhRv...51..947W, résumé).
  • [Bohr et Wheeler 1939] (en) Niels Bohr et John Archibald Wheeler, « The mechanism of nuclear fission » [« Le mécanisme de fission nucléaire »], Physical Review, vol. 56, no 5,‎ , p. 426-450 (OCLC 4644057022, DOI 10.1103/PhysRev.56.426, Bibcode 1939PhRv...56..426B, résumé, lire en ligne [PDF]).

Manuels d'enseignement supérieur

  • [Basdevant 2016] Jean-Louis Basdevant, La physique quantique et ses applications, Paris, Vuibert, coll. « LMD / Physique », , 1re éd., 1 vol., IX-421, ill., fig. et graph., 17 × 24 cm (ISBN 978-2-311-40039-7, EAN 9782311400397, OCLC 957581835, BNF 44272737, SUDOC 194834077, présentation en ligne, lire en ligne), § 5.2.2 (« La formule de masse de Bethe-Weizsäcker »), p. 102-105.
  • [Basdevant, Rich et Spiro 2002] Jean-Louis Basdevant, James Rich et Michel Spiro, Énergie nucléaire, Palaiseau, École polytechnique, hors coll., (réimpr. ), 1 vol., 336, ill., 17 × 24 cm (ISBN 978-2-7302-0901-4, EAN 9782730209014, OCLC 64768914, BNF 38862997, SUDOC 069243158, présentation en ligne, lire en ligne), § 2.2 (« La formule de masse de Bethe-Weizsäcker »), p. 67-71.
  • [Le Sech et Ngô 2014] Claude Le Sech et Christian Ngô, Physique nucléaire : des quarks aux applications, Paris, Dunod, coll. « Sciences Sup / Physique », , 2e éd. (1re éd. ), 1 vol., XVI-271, ill., 24 cm (ISBN 978-2-10-070541-2, EAN 9782100705412, OCLC 881576971, BNF 43822194, SUDOC 178846775, présentation en ligne, lire en ligne), § 2.9 (« Le modèle de la goutte liquide »), p. 29-31. Ouvrage utilisé pour la rédaction de l'article
  • [Pfalzner 2006] (en) Susanne Pfalzner, An introduction to inertial confinement fusion [« Une introduction à la fusion nucléaire par confinement inertiel »], New York et Londres, Taylor & Francis, coll. « Series in plasma physics », (réimpr. ), 1re éd., 1 vol., 232, ill., 16 × 24 cm (ISBN 978-0-7503-0701-7 et 978-0-3673-9109-6, EAN 9780750307017, OCLC 62679899, BNF 41099603, DOI 10.1201/9781420011845, SUDOC 112896030, présentation en ligne, lire en ligne), B.4 (« Constants of the semiempirical mass formula ») [« Constantes de la formule semi-empirique de masse »], p. 217. Ouvrage utilisé pour la rédaction de l'article
  • [Reuss 2003] Paul Reuss, Précis de neutronique, Les Ulis, EDP Sciences et INSTN, coll. « Génie atomique », (réimpr. ), 1re éd., 1 vol., 533, ill., fig. et graph., 17 × 24 cm (ISBN 2-86883-637-2, EAN 9782868836373, OCLC 470236789, BNF 39042390, SUDOC 074172719, présentation en ligne, lire en ligne), § 2.2.4 (« Modèle de goutte »), p. 49-51. Ouvrage utilisé pour la rédaction de l'article

Dictionnaires et encyclopédies

  • [Taillet, Villain et Febvre 2018] Richard Taillet, Loïc Villain et Pascal Febvre, Dictionnaire de physique, Louvain-la-Neuve, De Boeck Supérieur, hors coll., , 4e éd. (1re éd. ), 1 vol., X-956, ill., fig. et graph., 17 × 24 cm (ISBN 978-2-8073-0744-5, EAN 9782807307445, OCLC 1022951339, BNF 45646901, SUDOC 224228161, présentation en ligne, lire en ligne), s.v.Weizsäcker (formule de), p. 785-786. Ouvrage utilisé pour la rédaction de l'article
  • [Wapstra 1958] (de + en) Aaldert Hendrik Wapstra, « Atomic masses of nuclides », dans Siegfried Flügge (éd.), Handbuch der Physik / Encyclopedia of physics [« Encyclopedie de physique »], t. XXXVIII, 1re part. : Äussere Eigenschaften der Atomkerne / External properties of atomic nuclei [« Propriétés externes des noyaux atomiques »], Berlin, Göttingen et Heidelberg, Springer, coll. « Kernphysik / Nuclear physics » (no 8 / 38 / 1), (réimpr. ), 1re éd., 1 vol., VI-471, ill., fig. et graph., 17 × 24,4 cm (ISBN 978-3-642-45902-3, OCLC 467895019, BNF 37353867, DOI 10.1007/978-3-642-45901-6, SUDOC 015349322, présentation en ligne, lire en ligne), chap. 1er [« Masses atomiques des nucléides »], p. 1-37 (OCLC 7322443398, DOI 10.1007/978-3-642-45901-6_1). Ouvrage utilisé pour la rédaction de l'article

Articles connexes

Liens externes

  • [Khater 2009] Antoine Khater, « Introduction à la physique nucléaire », dans Ressources numériques pédagogiques pour la licence physique, (lire en ligne), § 2.10 (« Relation empirique pour l'énergie de liaison par nucléon ») (lire en ligne). Ouvrage utilisé pour la rédaction de l'article
  • [Valentin n.d.] Luc Valentin, « Nucléaire (physique) : noyau atomique », dans Encyclopædia Universalis, n.d. (lire en ligne), § 2 (« Masse et énergie de liaison ») (lire en ligne).
  • icône décorative Portail de la physique