Analytical Charge Density Profile of Vortex Core in Weak-Coupling Superconductor

Technische Samenvatting: Analytisch Profiel van de Ladingsdichtheid in de Vortexkern van een Zwak-gekoppelde Supergeleider

Probleemstelling
De elektronische structuur van vortexkernen in type-II supergeleiders staat erom bekend discrete quasipartikel-gebonden toestanden (Caroli–de Gennes–Matricon of CdGM-toestanden) te bevatten. Hoewel zelfconsistente Bogoliubov–de Gennes (BdG) berekeningen gekoppeld aan de Poisson-vergelijking numeriek hebben aangetoond dat de ladingsdichtheid binnen een vortexkern oscillerende tekens wisselt met een periode van $\pi/k_F$, heeft een eenvoudige analytische verklaring voor dit fenomeen lang ontbrak. Eerdere numerieke studies (bijv. Machida en Koyama) observeerden dat de afgeschermde ladingsdichtheid van teken verandert, maar het mechanisme dat deze tekenomkeer drijft—specifiek de interactie tussen de ladingsverdeling van de gebonden toestand en de elektrostatische afscherming—was nog niet vanuit eerste principes in een gesloten vorm afgeleid. De uitdaging ligt in de complexiteit van de BdG-vergelijkingen, waarbij normalisatieconstanten doorgaans mode-afhankelijk zijn en de golffuncties gekoppelde nietlineaire differentiaalvergelijkingen bevatten, wat analytische sommatie bemoeilijkt.

Methodologie
De auteur hanteert een elementaire, verklarende benadering met behulp van het meest economische model dat de essentiële fysica van de vortexkern behoudt:

1. Model Systeem: Een enkele vortexlijn in een zwak-gekoppelde s-golf supergeleider met een stapvormig gap-profiel ($\Delta(r) = 0$ voor $r < \xi$ en $\Delta(r) = \Delta_\infty$ anders). Het spectrum wordt genomen als de standaard CdGM-resultaat, $E_\mu = \mu \Delta^2/E_F$.
2. Constructie van de Golffunctie:
   • Binnenregio ($r < \xi$): De gap verdwijnt, waardoor de BdG-vergelijkingen ontkoppelen. De auteur gebruikt exacte Bessel-functieoplossingen ($J_{\mu \mp 1/2}$) in plaats van de WKB/BKJT-ansatz, die slecht gedefinieerd is nabij de vortexlijn vanwege de topologische winding.
   • Buitenregio ($r > \xi$): De BKJT-envelope wordt gelineraliseerd, waarbij de coëfficiënten worden vastgesteld door de matchingscondities bij de kernrand.
3. Normalisatie-analyse: Een cruciale stap is het berekenen van de normalisatieconstante $C_\mu$. De auteur demonstreert dat, ondanks de variërende klassiek toegestane vensters en vervalratio's, de normalisatieconstante bijna onafhankelijk is van het impulsmomentum kwantumgetal $\mu$ voor het grootste deel van de CdGM-ladder. Dit maakt het mogelijk om $C_\mu$ te behandelen als een constante prefactor $\langle C^2 \rangle$ in de mode-som.
4. Mode-sommatie:
   • Binnen de Kern: De ladingsdichtheidssom wordt gereduceerd tot een gedeeltelijke Bessel–Parseval-som. Door de volledigheidsidentiteit voor Bessel-functies toe te passen en op te merken dat de topologische beperking het $\nu=0$ kanaal uitsluit, stort de som in tot een gesloten vorm.
   • Buiten de Kern: De som levert een vervallende $1/r$ envelope op met een residuele rimpeling.
5. Elektrostatica: De ladingsdichtheid van de gebonden toestand dient als bronterm voor de Poisson-vergelijking. De auteur lost de afgeschermde Poisson-vergelijking op in de impulsruimte, waarbij de gladde achtergrond en de $2k_F$ oscillerende componenten apart worden behandeld.

Kernbijdragen en Resultaten

• Gesloten Vorm Binnenste Dichtheid: De auteur leidt een analytische expressie af voor de ladingsdichtheid van de gebonden toestand binnen de kern: $\delta n(r < \xi) \propto 1 - J_0^2(k_F r)$. Dit resultaat voldoet aan de randvoorwaarde $n(0)=0$ (door de uitsluiting van het $\nu=0$ Bessel-kanaal) en onthult dat de dichtheid bestaat uit een gladde achtergrond plus Friedel-achtige oscillaties met golfvector $2k_F$.
• Afschermingsmechanisme: Het artikel verheldert waarom de afgeschermde ladingsdichtheid van teken wisselt. De ladingsverdeling van de gebonden toestand integreert niet naar nul; daarom vereist algehele ladingneutraliteit dat uitgebreide toestanden dit compenseren.
  • Thomas–Fermi-afscherming is zeer efficiënt bij lange golflengten (lage momentumoverdracht), waardoor de gladde achtergrondcomponent van de lading volledig wordt geneutraliseerd.
  • Echter, afscherming is inefficiënt bij grote momentumoverdrachten ($q \approx 2k_F$). Bijgevolg overleeft de $2k_F$ oscillerende component van de ladingsverdeling van de gebonden toestand het afschermingsproces.
• Analytische Expressie voor de Overlevende Lading: De uiteindelijke afgeschermde ladingsdichtheid wordt gegeven door:
  $$\delta n_{\text{total}}(r) \simeq -\langle C^2 \rangle \frac{4k_F^2}{4k_F^2 + k_{TF}^2} \frac{\sin(2k_F r)}{\pi k_F r}$$
  Deze expressie laat zien dat de overlevende oscillatie een van teken wisselende $2k_F$ modulatie is. De amplitude wordt verminderd met een factor $4k_F^2/(4k_F^2 + k_{TF}^2)$, die tussen de $1/2$ en $3/4$ ligt voor typische metalen.
• Validatie: De analytische resultaten worden geverifieerd tegen directe numerieke oplossingen van de Poisson-vergelijking met behulp van de volledige CdGM mode-som. De analytische gesloten vormen komen overeen met de numerieke gegevens tot op 1% (waarbij de resterende fout wordt toegeschreven aan de mode-gemiddelde benadering van de normalisatieconstante).

Betekenis en Claims
Het artikel beweert een elementaire account te bieden van de van teken wisselende ladingsdichtheidsprofielen die in numerieke simulaties worden waargenomen. De primaire betekenis ligt in:

1. Mechanistische Helderheid: Het identificeert dat de oscillatie geen delicaat effect is, maar een direct gevolg van ladingneutraliteit en de inefficiëntie van Thomas–Fermi-afscherming bij grote momentumoverdrachten ($2k_F$).
2. Analytische Benchmark: Het biedt een gesloten vorm analytische benchmark (Gelijkheden 18, 22 en 28) waartegen toekomstige zelfconsistente numerieke berekeningen kunnen worden gecontroleerd, met name in het zwak-gekoppelde regime waar het gap-profiel rigide is.
3. Onderscheid van Onzuiverheidsverstrooiing: De auteur verduidelijkt dat hoewel de periodiciteit ($2k_F$) hetzelfde is als Friedel-oscillaties rond onzuiverheden, het mechanisme omgekeerd is. Bij onzuiverheidsverstrooiing is de bron vormloos en ontstaat de oscillatie door de niet-analyticiteit van de responsfunctie. In de vortexkern is de bron zelf gepiekt bij $2k_F$, en de respons (afscherming) filtert slechts de gladde componenten weg.

De auteur stelt expliciet dat het doel verklarend is en niet voorspellend, en claimt niet de kwantitatieve details in het sterk-gekoppelde regime of voor complexe gap-symmetrieën te hebben opgelost, hoewel hij deze als natuurlijke uitbreidingen merkt op. Het werk bevestigt dat de "rimpels" in de ladingsdichtheid robuuste kenmerken zijn van de vortexkern in zwak-gekoppelde supergeleiders.