Superfluid dome in the spatially modulated two-dimensional XY model

Technische Samenvatting: Superfluïde Koepel in het Ruimtelijk Gemoduleerde Twee-dimensionale XY-model

Probleemstelling
In sterk gecorreleerde elektronensystemen is de wisselwerking tussen supergeleidendheid (SC) en ladingsdichtheidsgolven (CDW) een centraal thema, wat zich vaak manifesteert als coëxisterende of verstrengelde orden. Recent theoretisch werk suggereert dat de supergeleidende ordeparameter kan koppelen aan modulaties in de elektronische ladingsdichtheid, wat leidt tot fenomenen zoals pair density waves en gefragmenteerde supergeleidende condensaten. Hoewel numerieke studies op Hubbard-modellen de begin fase van deze toestanden hebben verkend, moet een minimaal model dat de specifieke mechanisme vastlegt van hoe de ruimtelijke modulatie van de koppelingssterkte de supergeleidende overgangstemperatuur ($T_c$) beïnvloedt, nog volledig worden gekarakteriseerd. Specifiek vereist de afhankelijkheid van $T_c$ van de golflengte van dergelijke modulaties, en de onderliggende topologische mechanismen die deze afhankelijkheid drijven, nader onderzoek.

Methodologie
De auteurs onderzoeken een twee-dimensionaal XY-model op een vierkant rooster met een ruimtelijk gemoduleerde koppelingssterkte $J_{ij}$ langs de $x$-richting, gedefinieerd als $J_{ij} = J_0 [1 + \alpha \sin(k x_i)]$, waarbij $\alpha$ de modulatie-amplitude is en $\lambda = 2\pi/k$ de modulatiegolflengte. De koppeling langs de $y$-richting blijft uniform ($J_0$). Dit model is afgeleid van een verstrengelde Ginzburg-Landau-theorie waarbij de supergeleidende amplitude wordt gemoduleerd door een periodieke ladingsdichtheid.

Om de fasediagram en kritische eigenschappen te analyseren, maken de auteurs gebruik van een combinatie van:

1. Tensornetwerk-methoden: Specifiek wordt het Variational Uniform Matrix Product State (VUMPS) algoritme gebruikt om de rij-naar-rij transferoperator voor het oneindige systeem op te lossen. De partitiefunctie wordt gerepresenteerd als een contractie van lokale tensoren afgeleid via een duale transformatie. Kritische temperaturen worden geïdentificeerd met behulp van de verstrengelingsentropie ($S_E$) van de fixed-point matrix product state (MPS), en correlatiefuncties worden berekend met behulp van onzuiverheidstensors.
2. Monte Carlo (MC) Simulaties: Metropolis- en Wolff-algoritmen worden gebruikt voor eindige systemen ($L \times L$ tot $L=128$) om vortexconfiguraties, interne energie en heliciteitsmoduli te bestuderen.
3. Gemoduleerde $q$-state Clock Modellen: Deze worden bestudeerd om onderscheid te maken tussen domeinwand-gedreven en vortex-gedreven transities, in het bijzonder voor $q \ge 5$.
4. Renormalisatiegroep (RG) Analyse: Een hybride benadering die een enkele Migdal-Kadanoff (MK) RG-stap combineert met Kramers-Wannier dualiteitsvergelijkingen wordt gebruikt om kritische temperaturen te schatten en de schaling van effectieve koppelingen te begrijpen.

Belangrijkste Resultaten

1. Niet-monotone "Superfluïde Koepel": De kritische temperatuur $T_c$ vertoont een niet-monotone, koepelvormige afhankelijkheid van de modulatiegolflengte $\lambda$. Voor een vaste modulatie-amplitude $\alpha$ neemt $T_c$ aanvankelijk toe met $\lambda$, bereikt een piek, en neemt vervolgens af wanneer $\lambda$ zeer groot wordt.
2. Verschuiving van de Piekpositie: Naarmate de modulatie-amplitude $\alpha$ toeneemt, verschuift de piek van de $T_c$-koepel naar grotere golflengtes.
3. Kritisch Gedrag:
   • Onder $T_c$ vertoont het systeem quasi-lange-afstandsorde (QLRO) met algebraïsche afname van correlatiefuncties ($G(r) \sim r^{-\eta}$) en een centrale lading $c=1$, consistent met de Berezinskii-Kosterlitz-Thouless (BKT) universaliteitklasse.
   • Boven $T_c$ nemen correlaties exponentieel af. De afname is anisotroop, waarbij correlaties langs de modulatierichting ($x$) sneller afnemen dan de correlaties loodrecht daarop ($y$).
   • De pieken in de specifieke warmte verbreden met toenemende $\lambda$, wat wijst op een uitgebreid proces van vortex-ontbinding.
4. Vortex Pinning Mechanisme: MC-simulaties laten zien dat vortices bij voorkeur lokaliseren in de "dalen" van de modulatie (regio's met zwakke koppeling).
   • Bij kleine $\lambda$ zijn de dalen smal, wat de vorming van vortices beperkt en leidt tot een lagere $T_c$ door sterke anisotropie.
   • Naarmate $\lambda$ toeneemt, worden de dalen breder, waardoor er meer vortex-entropie mogelijk is. Echter, de effectieve interactiekracht wordt ook gerenormaliseerd. De koepelvorm ontstaat uit een competitie: aanvankelijk versterkt het vergroten van $\lambda$ de effectieve interactie (verhoging van $T_c$), maar verdere toenames faciliteren vortex-proliferatie en entropie (verlaging van $T_c$).
5. Vergelijking met Clock Modellen: In gemoduleerde $q$-state clock modellen neemt de lagere overgangstemperatuur $T_{c1}$ (gedreven door domeinwanden) monotoon toe met $\lambda$. In contrast hiermee vertoont de hogere overgangstemperatuur $T_{c2}$ (gedreven door vortices) voor $q \ge 5$ dezelfde koepelvormige afhankelijkheid als het XY-model, wat bevestigt dat de vortex-fysica de oorsprong is van de koepel.
6. Beperkingen van de RG-Analyse: De hybride MK-RG benadering slaagt erin de monotone toename van de effectieve koppeling met $\lambda$ en de onderdrukking van orde met toenemende $\alpha$ te vangen. Het overschat echter $T_c$ in de XY-limiet en slaagt er niet in de koepelstructuur te vatten zonder expliciet rekening te houden met vortex-excitaties, wat de noodzaak van topologische defecten benadrukt om het fenomeen te verklaren.

Betekenis en Claims
Het artikel claimt een microscopisch begrip te bieden van de "superfluïde koepel" die wordt waargenomen in de context van verstrengelde supergeleidende en ladingsgeordende toestanden. Door het probleem te mappen naar een gemoduleerd XY-model, demonstreren de auteurs dat de niet-monotone afhankelijkheid van $T_c$ van de modulatiegolflengte een generiek kenmerk is dat voortvloeit uit de effectieve pinning van vortices in de zwakke-koppelingsregio's van de modulatie.

De bevindingen suggereren dat het complexe fasedrag in materialen waar CDW en SC coëxisteren (zoals cupraten en kagome-metalen) begrepen kan worden via de thermodynamische competitie tussen vortex-energie en entropie in een ruimtelijk inhomogene omgeving. Het werk biedt een theoretisch kader waarin de modulatiegolflengte (gerelateerd aan de dopingdichtheid) fungeert als een controleparameter voor het optimaliseren van supergeleidendheid, waarbij een optimale golflengte voortkomt uit de wisselwerking tussen pinning en thermische fluctuaties. De auteurs stellen voor dat experimentele technieken zoals scanning tunneling microscopy de voorspelde vortex-pinning aan ladingsdichtheidsstrepen kunnen verifiëren.