Thomas-Fermi equation revisited: A variation on a theme by Majorana

Cet article réexamine la méthode de Majorana pour transformer l'équation de Thomas-Fermi en une équation du premier ordre, en l'appliquant aux atomes neutres et faiblement ionisés afin de recalculer et de comparer des grandeurs physiques clés avec des résultats numériques antérieurs.

L'essentiel

Le Titre : Une vieille recette, un nouveau chef

Imaginez que vous essayez de cuisiner un gâteau très complexe (c'est l'atome) avec une recette mathématique vieille de presque un siècle, appelée l'équation de Thomas-Fermi. Cette recette est connue pour être difficile à suivre : elle demande de résoudre une équation de "deuxième ordre", ce qui est un peu comme essayer de deviner la trajectoire d'une balle de tennis en tenant compte non seulement de sa vitesse, mais aussi de comment sa vitesse change à chaque instant, tout en faisant des calculs infinis.

En 1927, deux physiciens ont créé cette recette. Mais un génie nommé Ettore Majorana (un physicien mystérieux qui a disparu dans les années 30) a vu quelque chose que personne d'autre n'avait remarqué. Il a découvert un "truc de magicien" : il pouvait transformer cette équation compliquée en une équation beaucoup plus simple, de "premier ordre". C'est comme si, au lieu de devoir calculer chaque virage de la route, il vous donnait une carte qui vous dit simplement : "Tourne à gauche, puis continue tout droit".

Le problème ? Majorana n'a jamais publié son astuce de son vivant. Ses notes ont été retrouvées plus de 60 ans après sa disparition.

Le but de l'article : Remettre l'outil au goût du jour

L'auteur de cet article, Englert, dit : "Attendez, reprenons l'outil de Majorana !"

Son objectif est double :

1. Réutiliser la méthode de Majorana pour les atomes "normaux" (ceux qui ont autant d'électrons que de protons, comme l'atome d'hydrogène ou d'hélium).
2. Étendre cette méthode à un nouveau cas : les atomes "faiblement ionisés" (ceux qui ont perdu un ou deux électrons, comme un ballon qui se dégonfle un tout petit peu).

L'analogie du "Miroir Magique"

Pour comprendre ce que fait Englert, imaginez que l'équation de Thomas-Fermi est un labyrinthe immense et sombre.

• La méthode ancienne (années 80) : C'était comme essayer de sortir du labyrinthe en tâtonnant dans le noir, pas à pas, en vérifiant chaque mur. C'était long, fastidieux et sujet aux erreurs.
• La méthode de Majorana (et d'Englert) : C'est comme si Majorana avait accroché un miroir magique au centre du labyrinthe. Ce miroir transforme le labyrinthe 3D complexe en une simple ligne droite 2D.

Englert utilise ce miroir pour :

1. Redessiner les courbes : Il calcule avec une précision incroyable des nombres clés (comme la forme exacte du nuage d'électrons autour du noyau).
2. Vérifier les anciens calculs : Il compare ses nouveaux résultats avec ceux des années 80. Résultat ? Ses calculs sont plus précis et confirment les anciens, mais il y est arrivé beaucoup plus vite et plus élégamment. C'est comme si vous aviez trouvé une formule pour calculer $\pi$ en 10 secondes au lieu de 10 heures.

Les deux types d'atomes étudiés

L'article traite de deux situations, comme deux types de voyages :

1. Le voyage infini (Atome neutre) : Imaginez un atome qui s'étend à l'infini. Englert utilise la méthode de Majorana pour décrire comment l'électron se comporte depuis le centre jusqu'à l'infini. Il utilise une série de nombres (une suite mathématique) qui converge très vite, comme une voiture qui accélère doucement puis atteint sa vitesse de croisière sans jamais s'arrêter.
2. Le voyage borné (Atome ionisé) : Imaginez un atome qui a perdu un électron. Il est comme un ballon qui se dégonfle et s'arrête à une certaine limite. Englert a adapté le "miroir magique" de Majorana pour ce cas précis aussi, ce qui n'avait jamais été fait avec autant de clarté auparavant.

Pourquoi est-ce important ? (La recette du gâteau)

Pourquoi se soucier de ces calculs mathématiques abstraits ? Parce qu'ils nous aident à comprendre l'énergie.

Dans la dernière partie de l'article, Englert utilise ses nouveaux nombres précis pour recalculer l'énergie de liaison des atomes (l'énergie qu'il faut pour les casser) et l'énergie d'ionisation (l'énergie pour arracher un électron).
C'est comme si, grâce à sa nouvelle méthode de mesure, il pouvait dire exactement combien de sucre il faut dans le gâteau pour qu'il soit parfait, alors que les anciens cuisiniers devaient faire des essais et des erreurs.

En résumé

Cet article est un hommage à un génie oublié (Majorana) et une démonstration de puissance mathématique.

• Le problème : Les équations des atomes sont trop compliquées.
• La solution : Utiliser une astuce de transformation découverte par Majorana il y a 90 ans.
• Le résultat : On obtient des réponses plus précises, plus rapidement, et on découvre que cette astuce fonctionne aussi bien pour les atomes "morts" (neutres) que pour les atomes "blessés" (ionisés).

C'est une belle histoire de science où l'on redécouvre un vieux secret pour mieux comprendre l'univers qui nous entoure, le tout dédié à un ami, Goong Chen, qui partage cette passion pour les astuces mathématiques.

Résumé technique

1. Problématique

L'équation de Thomas-Fermi (TF), une équation différentielle du second ordre non linéaire, est fondamentale pour la modélisation des atomes à plusieurs électrons. Bien que la solution pour les atomes neutres soit bien connue, le calcul de certaines constantes physiques et intégrales associées a historiquement nécessité des procédures numériques lourdes et fastidieuses (notamment dans les années 1980).

De plus, il existe une seconde solution de l'équation TF, pertinente pour les atomes faiblement ionisés, qui n'avait pas été traitée avec la même élégance mathématique que le cas neutre. L'objectif de l'article est de réexaminer la méthode de Majorana (découverte dans ses notes privées publiées tardivement) pour transformer l'équation TF du second ordre en une équation du premier ordre, et d'appliquer cette méthode aux deux cas (atomes neutres et faiblement ionisés) afin de recalculer avec précision des grandeurs physiques clés.

2. Méthodologie

L'auteur s'appuie sur les propriétés d'échelle de l'équation TF pour réduire sa complexité.

• Transformation de Majorana : Au lieu de résoudre directement l'équation $f''(x) = x^{-1/2}f(x)^{3/2}$, l'auteur introduit trois fonctions invariantes d'échelle ($P, Q, R$) construites à partir de la solution $f(x)$ et de ses dérivées.
  • Pour le cas des atomes neutres (condition aux limites $f(0)=1, f(\infty)=0$), il exprime $R$ en fonction de $P$, ce qui conduit à une équation différentielle du premier ordre pour une fonction auxiliaire $u(t)$.
  • Pour le cas des atomes faiblement ionisés (condition aux limites $f(1)=0, f'(1)=-\Lambda^2$), il établit une équation analogue pour une fonction $v(s)$ reliant $Q$ à $P$.
• Développement en série de puissances : L'avantage majeur de cette transformation est que les fonctions $u(t)$ et $v(s)$ admettent des développements en série de puissances (autour de $t=1$ et $s=1$ respectivement) qui convergent sur tout l'intervalle de définition $[0, 1]$. Cela contraste avec les séries classiques de la fonction TF, qui ont des intervalles de convergence finis et non chevauchants.
• Calcul numérique : L'auteur utilise les relations de récurrence pour les coefficients de ces séries (déjà partiellement établies par Esposito pour le cas neutre) pour calculer numériquement les constantes avec une grande précision (arithmétique double précision).

3. Contributions Clés

1. Extension de la méthode de Majorana : L'article applique pour la première fois la transformation de Majorana au cas des atomes faiblement ionisés (solution (ii)), dérivant l'équation du premier ordre correspondante et sa série de puissances.
2. Unification des calculs : La méthode permet de traiter les deux régimes physiques (neutre et ionisé) avec la même approche mathématique élégante, évitant les méthodes numériques lourdes des années 1980.
3. Recalcul de constantes fondamentales : L'auteur recalcule les paramètres critiques $B$ (pente initiale pour l'atome neutre) et $\Lambda$ (paramètre pour l'atome ionisé), ainsi que les coefficients $\beta$ et $\alpha$ liés au comportement asymptotique.
4. Évaluation d'intégrales complexes : La méthode permet d'évaluer efficacement des intégrales impliquant des puissances de la fonction TF et de sa dérivée, nécessaires pour les formules d'énergie.

4. Résultats Principaux

Les calculs effectués confirment et affinent les valeurs obtenues dans les années 1980 :

• Constantes pour l'atome neutre :
  • $B = 1.588\,071\,022\,61$
  • $\beta = 13.270\,973\,848$
  • La fonction $u(t)$ est développée en série avec 20 coefficients listés dans le tableau 1.
• Constantes pour l'atome faiblement ionisé :
  • $\Lambda = 32.729\,416\,116\,173$
  • $\alpha = 1.040\,180\,657\,3862$
  • La fonction $v(s)$ est également développée en série avec une convergence garantie.
• Intégrales et Énergies :
  • Les intégrales de produits de puissances de $x$ et $F(x)$ (cas neutre) et $\Phi(x)$ (cas ionisé) sont recalculées avec une précision de 13 à 16 chiffres significatifs.
  • Énergie de liaison : Les coefficients $c_7$ et $c_5$ de l'expansion de l'énergie de liaison ($E \sim c_7 Z^{7/3} - \frac{1}{2}Z^2 + c_5 Z^{5/3}$) sont confirmés : $c_7 \approx 0.7687$ et $c_5 \approx 0.2699$.
  • Correction relativiste : Les coefficients de la correction relativiste ($c_{rel}$) sont recalculés, notamment $c_{rel}^1 \approx 2.075$.
  • Énergie d'ionisation : Pour les grands $Z$, l'énergie de première ionisation est recalculée à $I \approx 0.1158$ unités atomiques ($\approx 3.15$ eV).

5. Signification et Perspectives

• Validation et Efficacité : Ce travail démontre que la méthode de Majorana, bien que connue, offre une voie bien plus simple et efficace pour obtenir des résultats numériques de haute précision que les méthodes directes de résolution numérique de l'équation du second ordre.
• Précision théorique : La confirmation des valeurs des années 1980 renforce la fiabilité des modèles de Thomas-Fermi en physique atomique.
• Ouverture vers d'autres problèmes : L'article suggère que cette transformation pourrait être appliquée à d'autres équations non linéaires de type Lane-Emden ou à des solutions TF avec des conditions aux limites intermédiaires (atomes partiellement ionisés avec $0 < q < 1$), un domaine encore inexploré.

En résumé, l'article de Englert revitalise une méthode mathématique historique pour fournir une base numérique solide et élégante pour la physique atomique moderne, en particulier pour les calculs d'énergie et de structure électronique.