Ah-SCDFT:A general approach for superconductivity with an-harmonic corrections

Technische Samenvatting: Ah-SCDFT – Een Algemene Benadering voor Supergeleiding met Anharmonische Correcties

Probleemstelling
Eerstelijnsstudies naar supergeleiding vertrouwen vaak op de harmonische benadering, waarbij anharmonische effecten (AHE) worden genegeerd. Hoewel de harmonische benadering computationeel efficiënt is, slaagt deze er niet in om cruciale fysische verschijnselen te vatten in systemen met zachte vrijheidsgraden, dynamisch instabiele roosters (geïndiceerd door imaginaire fononfrequenties), of grote atomaire verplaatsingen bij hoge temperaturen. Het weglaten van AHE kan leiden tot significante kwantitatieve fouten en kwalitatieve misclassificaties, zoals het onjuist voorspellen van thermodynamische instabiliteit, het verkeerd identificeren van ladingsdichtheidsgolven, of het substantieel onderschatten van de elektron-fonon-koppelingssterkte (EPC) en de supergeleidende kritische temperatuur ($T_c$). Bestaande methoden zoals de McMillan/Allen-Dynes formules vertrouwen op empirische parameters, terwijl standaard anisotrope Migdal-Eliashberg vergelijkingen en Superconducting Density Functional Theory (SCDFT) doorgaans binnen het harmonische kader opereren.

Methodologie: Ah-SCDFT
De auteurs introduceren "anharmonische supergeleidende dichtheidsfunctionaaltheorie" (ah-SCDFT), een verenigde computationele benadering die systematisch anharmonische correcties integreert in de standaard SCDFT. De methodologie verloopt via een iteratieve zelfconsistente lus:

1. Initiële Berekening: Een grondtoestandsgeometrie ($S_1$) wordt verkregen bij 0 K met standaard Density Functional Theory (DFT) binnen de harmonische benadering.
2. Supergeleidende Evaluatie: Supergeleidende eigenschappen, inclusief een initiële kritische temperatuur ($T_{c1}$), worden berekend met gevestigde formalismen (anisotrope Migdal-Eliashberg vergelijkingen of SCDFT).
3. Rooster-renormalisatie: Een temperatuurafhankelijke evenwichtsstructuur wordt bepaald voor de voorspelde $T_{c1}$. Dit wordt bereikt door de Helmholtz vrije energie ($F = F_0 + F_{vib}$) te minimaliseren, waarbij $F_{vib}$ de vibratievrije energie is die berekend is met evenwichtsstructuren bij eindige temperatuur via de quasi-harmonische benadering (QHA) en de thermo_pw code.
4. Iteratie: Supergeleidende berekeningen worden herhaald met deze gerenormaliseerde geometrie om een bijgewerkte kritische temperatuur ($T_{c2}$) te verkrijgen.
5. Convergentie: Dit proces iteratert totdat het verschil tussen opeenvolgende kritische temperaturen voldoet aan $|T_{cn} - T_{cn-1}| < \epsilon$ (bijv. 0,01 K), wat een definitieve zelfconsistente $T_c$ oplevert.

De berekeningen maken gebruik van de QUANTUM ESPRESSO (QE) package en de SUPERCONDUCTING-TOOLKIT, waarbij de PBE generalized gradient approximation en de Sanna2020 exchange-correlation kernel worden toegepast.

Belangrijkste Resultaten: MgB2 Benchmark
De auteurs valideren de ah-SCDFT benadering met behulp van de prototypische s-golf supergeleider MgB2 onder zowel atmosferische condities als externe druk.

• Atmosferische Condities:

  • Harmonische Benadering: Voorspelt een $T_c$ van 36 K, consistent met eerdere rapporten maar lager dan de experimentele waarde.
  • Ah-SCDFT: Door het integreren van anharmonische rooster-renormalisatie bij 39 K, stijgt de voorspelde $T_c$ naar 39 K, wat overeenkomt met de experimentele gegevens.
  • Mechanisme: De verbetering wordt toegeschreven aan het feit dat anharmonische effecten de elektronische dichtheid van toestanden (DOS) bij het Fermi-niveau verhogen. Dit vergroot de EPC-constante ($\lambda$) van 0,82 (harmonisch) naar 0,88 (anharmonisch). Specifiek verzachten de anharmonische correcties de fononvertakkingen geassocieerd met Mg- en B-atomen, met name nabij het A-punt waar de in-plane B-vibraties de sterkste koppeling vertonen.

• Onder Externe Druk:

  • De studie behandelt conflicterende literatuur over de drukafhankelijkheid van de $T_c$ van MgB2 (sommige rapporten suggereren een monotone afname, andere een niet-monotone trend met een omslagpunt).
  • Ah-SCDFT Bevindingen: De benadering reproduceert een niet-monotone trend waarbij $T_c$ eerst afneemt en daarna herstelt, maar verschuift het omslagpunt naar aanzienlijk hogere drukken (rond 200 GPa) vergeleken met eerdere op harmonische principes gebaseerde voorspellingen (~100 GPa).
  • Drukeffecten: Naarmate de druk toeneemt, verstijft het rooster, waardoor het verschil tussen harmonische en anharmonische $T_c$-voorspellingen afneemt. De methode voorspelt ook structurele instabiliteit boven 340 GPa, gesignaleerd door imaginaire fononfrequenties.

Betekenis en Claims
Het artikel claimt dat ah-SCDFT een krachtig, efficiënt en breed toepasbaar kader vestigt voor de kwantitatieve studie van supergeleiding. De primaire betekenis ligt in:

1. Nauwkeurigheid: Het bereikt hoogwaardige voorspellingen van supergeleidende eigenschappen (specifiek $T_c$) door expliciet rekening te houden met temperatuurafhankelijke rooster-renormalisatie, waardoor discrepanties tussen theorie en experiment in materialen zoals MgB2 worden opgelost.
2. Efficiëntie: De methode brengt slechts een bescheiden toename van de computationele kosten met zich mee, aangezien een beperkt aantal convergentiestappen voor de iteratieve lus vereist is, wat de toepassing voor bredere toepassingen haalbaar maakt.
3. Betrouwbaarheid: Het biedt een consistente benadering voor het begrijpen van drukafhankelijke supergeleiding en biedt een robuust hulpmiddel voor de voorspelling en het ontwerp van nieuwe supergeleidende materialen waar anharmonische effecten een cruciale rol spelen.