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

Dit artikel heronderzoekt Majorana's methode om de Thomas-Fermi-vergelijking via schaaltransformaties te reduceren tot een eerste-orde differentiaalvergelijking, past deze toe op zowel neutrale als zwak geïoniseerde atomen, en berekent diverse belangrijke grootheden opnieuw met resultaten die overeenkomen met eerdere, meer tijdrovende numerieke procedures.

De kern

De Thomas-Fermi vergelijking opnieuw bekeken: Een verhaal over het vinden van een kortere weg

Stel je voor dat je een enorme, ingewikkelde berg moet beklimmen. Deze berg vertegenwoordigt een heel oud en belangrijk probleem in de natuurkunde: hoe gedragen zich de elektronen rondom een atoom? In de jaren '20 bedachten twee wetenschappers, Thomas en Fermi, een wiskundige kaart (een vergelijking) om deze berg te beschrijven. Maar deze kaart was lastig: het was een "tweede-orde" vergelijking. In het wilde woud van de wiskunde betekent dit dat je twee keer moet "stappen" om de top te bereiken, wat de berekeningen erg zwaar en traag maakt.

De mysterieuze Majorana en zijn geheime tunnel
In 1928 keek een genie genaamd Majorana naar deze kaart en dacht: "Wacht even, er is een trucje." Hij zag dat de berg bepaalde symmetrieën had (als je de kaart vergroot of verkleint, blijft de vorm hetzelfde). Hij bedacht een manier om de lange, moeilijke klim om te zetten in een kortere, rechte tunnel: een "eerste-orde" vergelijking.

Helaas vertelde Majorana zijn geheim aan niemand. Hij verdween mysterieus en zijn aantekeningen werden pas zestig jaar later gevonden. Pas toen konden wetenschappers zien hoe hij de tunnel had gebouwd.

Wat doet deze nieuwe paper?
De auteur van dit artikel, Berthold-Georg Englert, zegt: "Laten we Majorana's tunnel niet alleen gebruiken voor de standaard berg (een neutraal atoom), maar ook voor een iets andere berg (een atoom dat een beetje elektrisch geladen is, een 'geïoniseerd' atoom)."

Hij doet drie dingen:

1. Hij bouwt een tweede tunnel: Hij past Majorana's methode toe op het geval van de geladen atomen, waarvoor er nog geen duidelijke tunnel was.
2. Hij meet alles opnieuw: Met deze nieuwe tunnels (de eerste-orde vergelijkingen) berekent hij allerlei belangrijke getallen. Vroeger moesten wetenschappers in de jaren '80 urenlang met zware computers rekenen om deze getallen te vinden. Englert doet het nu veel sneller en eleganter.
3. Hij controleert het werk: Hij kijkt of zijn nieuwe, snelle resultaten overeenkomen met de oude, zware berekeningen. En ja, ze kloppen perfect!

De analogie: De rechte lijn versus de slingerweg
Stel je voor dat je een auto hebt die een steile heuvel moet oprijden.

• De oude methode (jaren '80): Je rijdt de hele steile, kronkelende weg op en neer. Het kost veel brandstof en tijd.
• Majorana's methode: Hij heeft een lift of een tunnel door de berg gebouwd. Je rijdt recht naar binnen, en je komt precies aan de andere kant uit. Het resultaat is hetzelfde, maar je bent veel sneller en je hebt minder brandstof (rekenkracht) nodig.

Englert heeft nu niet alleen de lift voor de normale berg gevonden, maar ook een nieuwe lift voor de "geïoniseerde" berg. Bovendien heeft hij laten zien dat je met deze liften heel precies kunt meten hoe hoog de berg is, zonder dat je de hele weg hoeft te rijden.

Waarom is dit belangrijk?
In de wereld van atoomfysica zijn deze getallen cruciaal. Ze vertellen ons hoe sterk atomen aan elkaar plakken (bindingsenergie) en hoeveel energie nodig is om een elektron los te maken (ionisatie-energie).

Door Majorana's slimme wiskundige truc te gebruiken, kunnen wetenschappers deze berekeningen nu veel makkelijker en nauwkeuriger doen. Het is alsof ze een oude, moeizame kaart hebben vervangen door een moderne GPS die de kortste route direct aangeeft.

Conclusie
Deze paper is een eerbetoon aan een genie uit het verleden (Majorana) en een bewijs dat slimme wiskundige ideeën, zelfs als ze al bijna 100 jaar oud zijn, nog steeds kunnen helpen om complexe problemen van vandaag sneller en beter op te lossen. Englert heeft de "Majorana-tunnel" uitgebreid en laten zien dat hij voor meer soorten atomen werkt dan eerder gedacht.

Technische samenvatting

Titel: Thomas-Fermi-vergelijking opnieuw bezocht: Een variatie op een thema van Majorana

Auteur: Berthold-Georg Englert
Doel: Heruitgave en uitbreiding van Majorana's methode voor het oplossen van de Thomas-Fermi-vergelijking, met name voor neutrale en zwak geïoniseerde atomen.

---

1. Het Probleem

De Thomas-Fermi (TF) vergelijking, geïntroduceerd in 1927/1928, is een niet-lineaire tweede-orde differentiaalvergelijking die het elektronendichteprofiel in atomen beschrijft. De vergelijking luidt:
$$f''(x) = x^{-1/2}f(x)^{3/2} \quad \text{voor } f(x) \geq 0$$
Hoewel deze vergelijking fundamenteel is voor de atoomfysica, is het analytisch oplossen ervan complex. Er zijn twee specifieke oplossingen van groot belang:

1. Neutrale atomen: De oplossing $F(x)$ met randvoorwaarden $f(0)=1$ en $f(\infty)=0$.
2. Zwak geïoniseerde atomen: De oplossing $\Phi(x)$ met randvoorwaarden $f(1)=0$ en $f'(1) = -\Lambda^2$.

Historisch gezien werden numerieke waarden voor belangrijke constanten (zoals de helling $B$ bij de oorsprong en integrals van de functie) in de jaren 80 berekend met tijdrovende numerieke methoden. De uitdaging is om deze waarden nauwkeuriger en efficiënter te berekenen door gebruik te maken van de schaal-eigenschappen van de vergelijking.

2. Methodologie

De kern van de methode is gebaseerd op het werk van Ettore Majorana, die ontdekte dat de schaal-invariantie van de TF-vergelijking kan worden gebruikt om de tweede-orde differentiaalvergelijking te reduceren tot een eerste-orde differentiaalvergelijking.

Stappen in de aanpak:

• Schaal-invariante variabelen: De auteur introduceert drie schaal-invariante combinaties van de oplossing $f(x)$ en zijn afgeleiden:
  • $P(x) = x^{3/2}f(x)^{1/2}$
  • $Q(x) = -x^4f'(x)$
  • $R(x) = -f(x)^{-4/3}f'(x)$
• Reductie tot eerste orde:
  • Voor het neutrale atoom (geval i) wordt gezocht naar $R$ als functie van $P$. Dit leidt tot Majorana's vergelijking voor een nieuwe functie $u(t)$:
    $$\frac{du}{dt} = -8 \frac{1 - t u(t)^2}{1 - t^2 u(t)}$$
    waarbij $t$ een gereduceerde variabele is ($P=12t^3$).
  • Voor het geïoniseerde atoom (geval ii) wordt $Q$ als functie van $P$ gezocht, wat leidt tot een analoog eerste-orde systeem voor een functie $v(s)$.
• Reeksontwikkelingen: In plaats van de oorspronkelijke TF-functie direct te benaderen, worden de nieuwe functies $u(t)$ en $v(s)$ benaderd via machtreeksen rondom de eindpunten ($t=1$ en $s=1$).
  • De reeksen voor $u(t)$ en $v(s)$ convergeren over het hele domein ($0 \leq t, s \leq 1$).
  • Dit staat in contrast met de reeksen voor de oorspronkelijke TF-functies, die slechts op beperkte, niet-overlappende intervallen convergeren.
• Numerieke Integratie: De fysieke grootheden worden berekend door integrals over het nieuwe domein ($t$ of $s$) uit te voeren, waarbij de integranden worden uitgedrukt in termen van de machtreeksen van $u(t)$ en $v(s)$.

3. Belangrijkste Bijdragen

• Herleving van Majorana's methode: De auteur formaliseert en verduidelijkt Majorana's techniek, die decennialang onbekend bleef totdat zijn notities werden gepubliceerd.
• Uitbreiding naar geïoniseerde atomen: Voor het eerst wordt Majorana's eerste-orde methode succesvol toegepast op de oplossing voor zwak geïoniseerde atomen (geval ii), wat een nieuwe eerste-orde differentiaalvergelijking oplevert.
• Efficiëntere berekeningsmethode: De auteur demonstreert dat het gebruik van deze eerste-orde vergelijkingen in combinatie met machtreeksen een veel eenvoudigere en robuustere manier is om numerieke waarden te berekenen dan de methoden uit de jaren 80.
• Nieuwe tabellen en coëfficiënten: Er worden nauwkeurige tabellen geproduceerd met coëfficiënten voor de machtreeksen van $u(t)$, $v(s)$ en de bijbehorende integralen.

4. Resultaten

De auteur berekent opnieuw diverse fundamentele constanten en integralen en vergelijkt deze met eerdere resultaten:

• Constanten voor neutrale atomen:
  • De helling bij de oorsprong: $B = 1.588\,071\,022\,61$.
  • De asymptotische coëfficiënt: $\beta = 13.270\,973\,848$.
  • Deze waarden worden bevestigd met hoge precisie (tot op de 13e decimaal of meer) via de nieuwe methode.
• Constanten voor geïoniseerde atomen:
  • De parameter $\Lambda = 32.729\,416\,116\,173$ en $\alpha = 1.040\,180\,657\,3862$.
• Integralen:
  • Verschillende integralen van de vorm $\int x^k F(x)^l dx$ worden nauwkeurig berekend.
  • Een specifieke correctie wordt opgemerkt in een eerdere publicatie (Ref. [4]) waar een typefout in de berekening van $\int (-F'(x))^3 dx$ was gemaakt.
• Energieformules:
  • De berekende integralen worden gebruikt om de coëfficiënten in de formules voor de bindingsenergie van neutrale atomen en de ionisatie-energie te verfijnen.
  • De ionisatie-energie voor grote $Z$ wordt berekend als $I \approx 0.1158$ Hartree ($\approx 3.15$ eV).

5. Betekenis en Conclusie

• Wiskundige elegantie: De paper toont aan dat de schaal-invariantie van de Thomas-Fermi-vergelijking een krachtig hulpmiddel is om de complexiteit van het probleem te reduceren. De overgang van een tweede-orde naar een eerste-orde vergelijking vereenvoudigt de numerieke implementatie aanzienlijk.
• Numerieke superioriteit: De methode van Majorana, ondersteund door machtreeksen die over het hele domein convergeren, is superieur aan traditionele numerieke integratie of het samenvoegen van reeksen met beperkte convergentiegebieden.
• Toekomstperspectief: De auteur suggereert dat deze methode ook kan worden toegepast op andere niet-lineaire vergelijkingen van het type $f''(x) = x^a f(x)^b$ (zoals de Lane-Emden vergelijking) en op de studie van atomen met een willekeurige mate van ionisatie ($0 < q < 1$), wat een nieuw onderzoeksgebied opent.

Kortom, dit artikel biedt een elegante, wiskundig onderbouwde en numeriek efficiënte herformulering van een klassiek probleem in de atoomfysica, bevestigend de waarde van historische inzichten (Majorana) voor moderne berekeningen.