First-principles calculation of coherence length and penetration depth based on density functional theory for superconductors

De auteurs presenteren een eerste-principes raamwerk binnen supergeleidende dichtheidsfunctionaaltheorie om de coherentielengte en magnetische penetratiediepte parameterloos te berekenen, wat leidt tot nauwkeurige voorspellingen voor diverse supergeleiders en een microscopische interpretatie van de Uemura-plot.

De kern

De Superhelden van de Wereld: Hoe Wetenschappers de "Onzichtbare Krachten" van Supergeleiding Berekenen

Stel je voor dat je een groep mensen hebt die door een drukke stad lopen. Normaal gesproken stoten ze elkaar, botsen ze tegen lantaarnpalen en raken ze in de war. Dit is hoe elektriciteit zich normaal gedraagt in een gewone draad: er is weerstand, en energie gaat verloren als warmte.

Maar in een supergeleider gebeurt er iets magisch. De mensen (elektronen) vormen paren en dansen perfect synchroon door de stad. Ze botsen niet meer, ze glijden als één enkel, onzichtbaar stroompje. Er is geen weerstand. Dit is het geheim van supergeleiding.

Deze nieuwe studie van Kawamura en zijn team is als het bouwen van een superkrachtige simulator om precies te begrijpen hoe deze dans werkt, zonder dat we eerst een echt lab nodig hebben.

Hier is wat ze hebben gedaan, vertaald in alledaags taal:

1. Het Probleem: De "Onzichtbare Liniaal"

Wetenschappers weten al hoe ze de temperatuur kunnen voorspellen waarop een materiaal supergeleidend wordt (de kritieke temperatuur). Maar er zijn twee andere cruciale maten die ze nog niet goed konden berekenen:

• De "Dansgrootte" (Coherentie Lengte): Hoe groot is de dansvloer waarop één koppel elektronen samenwerkt? Is het een kleine dansvloer in een kamer, of een gigantische dansvloer over de hele stad?
• De "Schildkracht" (Penetratiediepte): Als je een magneet op de supergeleider legt, hoe diep dringt het magnetische veld erin? Is het als een schild dat alles buiten houdt, of dringt het diep door?

Zonder deze maten is het moeilijk om nieuwe supergeleiders te ontwerpen voor echte toepassingen, zoals snellere treinen of krachtige MRI-scanners.

2. De Oplossing: Een Nieuwe Simulator

De onderzoekers hebben een nieuwe manier bedacht om dit te berekenen, gebaseerd op DFT (Dichtheidsfunctionaaltheorie). Denk hierbij aan een extreem krachtige computer die de atomen en elektronen van een materiaal naait, tot in de kleinste details.

Maar ze hadden een probleem: om de "dansgrootte" te meten, moeten ze kijken naar paren die een heel klein beetje "schuiven" of bewegen. In de computerwereld is dit als proberen een auto te meten terwijl hij met een snelheid van 0,0001 km/u rijdt. Dat is heel lastig om numeriek stabiel te houden; de computer wordt snel duizelig.

Hun truc: Ze hebben een slimme wiskundige "bril" op de computer gezet. Hierdoor kunnen ze die heel kleine bewegingen (de schuifbeweging van de elektronenparen) stabiel berekenen. Ze noemen dit het introduceren van "Cooper-paren met eindige impuls". Klinkt ingewikkeld, maar het betekent simpelweg: "Laten we kijken wat er gebeurt als de dansparen een heel klein beetje vooruitlopen."

3. De Resultaten: Van Aluminium tot H3S

Ze hebben hun simulator getest op verschillende materialen:

• De Klassiekers: Materialen zoals Aluminium en Niobium. Hier klopten hun berekeningen perfect met de echte meetresultaten uit het lab. Het was alsof hun simulator een perfecte voorspelling deed van hoe groot de dansvloer was.
• De Uitdaging: Ze keken ook naar H3S (waterstofsulfide). Dit is een heel speciaal materiaal dat alleen onder enorme druk (zoals in het binnenste van een planeet) supergeleidend wordt. Het is bijna onmogelijk om daar in een lab de "dansgrootte" te meten.
  • Het resultaat: De simulator gaf een antwoord dat perfect overeenkwam met de beperkte metingen die er wel waren. Dit betekent dat we nu kunnen voorspellen hoe supergeleiders werken in omstandigheden die voor mensen te gevaarlijk of onbereikbaar zijn.

4. De Grote Ontdekking: De "Uemura-Plot"

Het mooiste deel van hun werk is dat ze een grote puzzel hebben opgelost. Er is een bekende grafiek (de Uemura-plot) die laat zien dat er een verband is tussen hoe heet een materiaal moet zijn om te supergeleiden en hoe "stijf" de dans van de elektronen is.

Vroeger was dit alleen een observatie: "Oh, als het zo is, dan is dat zo."
Met hun nieuwe simulator kunnen ze dit voorspellen vanuit de basiswetten van de natuurkunde. Ze laten zien:

• Gewone supergeleiders (zoals Aluminium) hebben een grote, losse dansvloer en een zwakke dansstijfheid.
• Krachtige supergeleiders (zoals H3S of koper-oxides) hebben een heel strakke, sterke dans (sterke koppeling) en een heel stijve dansvloer.

Waarom is dit belangrijk?

Stel je voor dat je een architect bent die een brug wil bouwen. Je kunt de brug niet bouwen als je niet weet hoe sterk het staal is of hoe diep de fundering moet zijn.

Voor supergeleiders is dit onderzoek hetzelfde. Nu kunnen we:

1. Nieuwe materialen vinden: We kunnen in de computer zoeken naar materialen die misschien nog nooit zijn ontdekt, maar die wel supergeleidend zijn bij kamertemperatuur.
2. Extreme omstandigheden begrijpen: We kunnen voorspellen wat er gebeurt in sterren of onder extreme druk, zonder dat we daarheen hoeven te reizen.
3. Technologie verbeteren: Van snellere treinen tot betere medische apparatuur, alles begint met het begrijpen van deze "onzichtbare linialen".

Kortom: Deze onderzoekers hebben een nieuwe, krachtige bril op de natuurkunde gezet. Hierdoor kunnen we niet alleen zien dat iets supergeleidend is, maar precies begrijpen hoe het werkt, zelfs in de meest extreme hoeken van het universum.

Technische samenvatting

Hier is een gedetailleerde technische samenvatting van het artikel "First-principles calculation of coherence length and penetration depth based on density functional theory for superconductors", geschreven in het Nederlands.

Titel: Eerste-principes berekening van coherentielengte en penetratiediepte voor supergeleiders op basis van dichtheidsfunctionaaltheorie

Auteurs: Mitsuaki Kawamura, Takuya Nomoto, Niklas Witt, en Ryotaro Arita.
Datum: 6 maart 2026 (voorgesteld).

---

1. Het Probleem

Hoewel superconducting density functional theory (SCDFT) al succesvol wordt gebruikt om de kritieke temperatuur ($T_c$) van supergeleiders voorspellen, ontbreekt er tot nu toe een robuust eerste-principes kader voor het berekenen van fundamentele lengteschalen die de supergeleidende toestand karakteriseren:

• De coherentielengte ($\xi_0$): Een maat voor de grootte van de Cooper-paren.
• De magnetische penetratiediepte ($\lambda_L$): Een maat voor hoe diep een magnetisch veld in de supergeleider doordringt.

Deze grootheden zijn cruciaal voor het classificeren van supergeleiders in Type-I en Type-II, het bepalen van kritieke magnetische velden, depairing-stromen en supergeleidende stijfheid. Bestaande methoden (zoals dynamische mean-field theorie) zijn vaak niet direct toepasbaar binnen het SCDFT-kader of vereisen fenomenologische input. Vooral bij materialen onder extreme omstandigheden (zoals hoge druk) is experimentele bepaling van deze waarden uiterst moeilijk, wat een betrouwbare theoretische methode noodzakelijk maakt.

2. Methodologie

De auteurs ontwikkelen een nieuw microscopisch kader binnen SCDFT dat Cooper-paren met een eindige impuls (finite-momentum) expliciet incorporeert.

• Theoretische Basis:

  • De methode begint met de Kohn-Sham-Bogoliubov-de Gennes (KS-BdG) vergelijkingen.
  • In plaats van alleen Cooper-paren met impuls $Q=0$ te beschouwen, wordt een eindige impuls $Q$ geïntroduceerd. Dit leidt tot een modulerende ordeparameter: $\chi(r+R, r'+R) = e^{iQ\cdot R}\chi(r, r')$.
  • Door een generalisatie van het Bloch-theorema en een decoupleringsbenadering, worden de golfvectoren verschoven naar $k \pm Q/2$.

• Numerieke Stabilisatie:

  • Omdat de coherentielengte vaak veel groter is dan de roostervaste, moet $Q$ extreem klein zijn ten opzichte van de $k$-punt-grid. Dit veroorzaakt numerieke instabiliteit.
  • De auteurs introduceren een Taylor-ontwikkeling voor de Kohn-Sham-energieën bij de verschoven punten ($k \pm Q/2$) om deze nauwkeurig te schatten zonder de grid te hoeven verfijnen.
  • Een hulp-energieafhankelijk gap-functie wordt gebruikt om de $k$-punt-integratie te stabiliseren, vooral in de buurt van het Fermi-niveau waar de integrand snel verandert.

• Berekening van Grootheden:

  • Coherentielengte ($\xi_0$): Wordt afgeleid uit de $Q$-afhankelijkheid van het gemiddelde supergeleidende gap-functie $\langle \Delta(Q) \rangle$. Volgens de Ginzburg-Landau-theorie daalt dit kwadratisch met $Q$: $\langle \Delta(Q) \rangle \propto \sqrt{1 - \xi_0^2 Q^2}$.
  • Penetratiediepte ($\lambda_L$): Wordt berekend uit de superstroomdichtheid $j_{sc}(Q)$ die door de Cooper-paren wordt gedragen. $\lambda_L$ is gerelateerd aan de helling van $j_{sc}$ bij $Q \to 0$ (London-benadering) of via de depairing-stroom (het maximum van $j_{sc}$).

3. Belangrijkste Resultaten

De methode is toegepast op een reeks materialen, waaronder elementaire supergeleiders (Al, Nb, Sn, In, Ta, Pb), het A15-verbinding $V_3Si$, en het hoge-temperatuur supergeleider $H_3S$ onder hoge druk.

• Overeenkomst met Experiment:

  • Voor Niobium (Nb) wordt een coherentielengte van 34 nm en een penetratiediepte van 34 nm berekend, wat uitstekend overeenkomt met experimentele waarden (respectievelijk ~39 nm en 27-39 nm).
  • Voor $H_3S$ (bij ~200 GPa) voorspellen de auteurs $\xi_0 \approx 3.0$ nm en $\lambda_L \approx 19-22$ nm. Dit bevestigt dat $H_3S$ een Type-II supergeleider is en komt overeen met experimentele schattingen die variëren tussen 20 en 37 nm.
  • De methode reproduceert correct de classificatie van Type-I ($\xi > \sqrt{2}\lambda$) versus Type-II ($\xi < \sqrt{2}\lambda$) voor alle onderzochte materialen.

• Depairing-stroom:

  • De berekende depairing-stroom ($J_{dp}$) voor $H_3S$ is extreem hoog ($697 \times 10^7$A/cm²), wat aanzienlijk hoger is dan bij koper-oxide supergeleiders. Dit suggereert dat$H_3S$ in staat is uitzonderlijk grote superstromen te dragen.

• De Uemura-plot:

  • Een uniek inzicht is de mogelijkheid om de Uemura-plot ($T_c$ versus Fermi-temperatuur $T_F$) volledig uit eerste-principes te construeren.
  • De resultaten tonen aan dat conventionele elementaire supergeleiders systematisch een kleine verhouding $T_c/T_F$ hebben.
  • Hogere-$T_c$ systemen (zoals $V_3Si$ en $H_3S$) worden gekenmerkt door een gelijktijdige realisatie van sterke koppeling (korte $\xi_0$) en grote fase-stijfheid (korte $\lambda_L$, hoge $T_F$). Dit biedt een microscopische verklaring voor empirische correlaties die eerder alleen fenomenologisch werden beschreven.

4. Bijdragen en Significantie

• Eerste-principes Voorspellend Vermogen: Dit werk vestigt een volledig parameter-vrije route om supergeleidende lengteschalen ($\xi_0, \lambda_L$) en kritieke stromen te voorspellen zonder empirische aanpassing.
• Extreme Omstandigheden: De methode is bijzonder waardevol voor materialen onder extreme druk (zoals $H_3S$), waar experimentele metingen van magnetische eigenschappen technisch zeer uitdagend of onmogelijk zijn.
• Unificatie: Het koppelt voor het eerst de berekening van $T_c$, $\xi_0$ en $\lambda_L$ binnen één consistent theoretisch raamwerk (SCDFT).
• Fundamenteel Inzicht: Door de Uemura-plot te verklaren vanuit microscopische principes, versterkt het artikel het begrip van de relatie tussen paringsterkte en fase-stijfheid in supergeleiders.

Conclusie:
De auteurs hebben een krachtige, numeriek stabiele SCDFT-methode ontwikkeld die de berekening van fundamentele supergeleidende lengteschalen mogelijk maakt. De uitstekende overeenkomst met beschikbare experimentele data, zelfs voor complexe systemen onder hoge druk, valideert de methode als een essentieel hulpmiddel voor het ontwerpen van nieuwe supergeleiders en het begrijpen van hun gedrag onder extreme omstandigheden.