Disorder averaging in random lattice models with periodic boundary conditions: Application to models with uncorrelated and correlated disorder

Deze paper ontwikkelt technieken voor het middelen van wanorde in willekeurige roostermodellen met periodieke randvoorwaarden om de variatie en hogere momenten van polarisatie te berekenen, en valideert deze methoden op zowel volledig gelokaliseerde als power-law-gecorreleerde wanordemodellen.

De kern

Stel je voor dat je een enorm, oneindig groot tapijt hebt. Dit tapijt is niet egaal; het is een wirwar van verschillende patronen, knopen en vlekken. In de natuurkunde noemen we dit een ongevormd systeem (een systeem met wanorde). De vraag die fysici zich al decennia stellen, is: Is dit tapijt een goede geleider van energie (zoals elektriciteit), of is het een isolator die de energie vasthoudt?

In een perfect, geordend tapijt (zoals een kristal) is het makkelijk om te zien hoe de deeltjes zich gedragen. Maar zodra je wanorde toevoegt, wordt het een chaos. De auteurs van dit artikel, Balázs Hetényi en zijn collega's, hebben een nieuwe manier bedacht om deze chaos te bestuderen, zelfs als je het tapijt in een eindige doos stopt (wat in computersimulaties noodzakelijk is).

Hier is een simpele uitleg van wat ze hebben gedaan, met behulp van alledaagse vergelijkingen:

1. Het Probleem: De "Eindige Doos" en de "Grote Wereld"

Wanneer wetenschappers computersimulaties doen, kunnen ze geen oneindig groot tapijt maken. Ze moeten werken met een eindig stukje. Om dit stukje te laten lijken op een oneindig tapijt, gebruiken ze vaak periodieke randvoorwaarden.

• De Analogie: Stel je een video-game wereld voor (zoals Pac-Man). Als je naar rechts loopt en de rand bereikt, verschijn je aan de linkerkant. Het is alsof het tapijt in een cirkel is gerold.
• Het Probleem: In zo'n cirkel is het lastig om te meten hoeveel "lading" (elektriciteit) ergens zit, omdat er geen begin en eind is. De auteurs gebruiken een slimme wiskundige truc (de moderne theorie van polarisatie) om dit op te lossen. Ze kijken niet naar de lading zelf, maar naar een geometrische fase.
• De Vergelijking: Denk aan het meten van de "draaiing" van een kompasnaald in een storm. Je meet niet hoe hard de wind waait, maar hoe de naald reageert op de windrichting. Dit geeft een veel nauwkeuriger beeld van de stormkracht.

2. De Nieuwe Tool: De "Statistiek van de Chaos"

De auteurs hebben een methode ontwikkeld om te middelen over duizenden verschillende versies van dit "tapijt".

• De Analogie: Stel je voor dat je 1000 keer een dobbelsteen gooit, maar elke keer is de dobbelsteen een beetje anders (sommige kanten zijn zwaarder, andere lichter). Je wilt weten: Is de dobbelsteen eerlijk (geordend) of is hij zo scheef dat hij nooit op een bepaalde kant landt (gevangen/lokaal)?
• Ze gebruiken een maatstaf genaamd de Binder-cumulant.
  • Als de waarde laag is, betekent dit dat de deeltjes vastzitten (geïsoleerd). Ze kunnen niet vrij bewegen, net als mensen in een drukke menigte die vastlopen in een smalle gang.
  • Als de waarde hoog is, betekent dit dat de deeltjes vrij bewegen (gedelokaliseerd). Ze kunnen over het hele tapijt zwermen, zoals mensen in een groot, open plein.

3. De Twee Experimenten

De auteurs testen hun methode op twee modellen:

A. Het Anderson-model (De "Willekeurige Muur")

• Het Systeem: Een rij huisjes waarbij elke deur een willekeurige vergrendeling heeft.
• Het Resultaat: Ze bevestigen wat we al wisten: als de vergrendelingen willekeurig genoeg zijn, blijven de deeltjes in één huisje hangen. Ze kunnen niet naar het volgende huisje. De "geometrische fase" toont dit perfect aan.

B. Het de Moura-Lyra-model (De "Gereguleerde Chaos")

• Het Systeem: Hier is de wanorde niet helemaal willekeurig. Er zit een patroon in, zoals een rimpeling in water die langzaam afneemt. Dit wordt bepaald door een parameter genaamd $\alpha$.
• De Vraag: Is er een punt waar de deeltjes plotseling van "vastzitten" naar "vrij bewegen" gaan? Dit wordt een mobiliteitsrand genoemd.
• De Ontdekking:
  • De auteurs vinden dat er inderdaad een overgang is bij een bepaalde waarde van $\alpha$.
  • Maar er is een raar fenomeen in het midden van het tapijt (bij een specifieke energie). Hier sluiten de "gaten" tussen de energieniveaus zich in paren.
  • De Analogie: Stel je voor dat je een dansvloer hebt. Normaal gesproken kunnen mensen overal lopen. Maar in dit specifieke gebied (waar $\alpha > 2$) vormen de mensen plotseling koppels die perfect in elkaars armen passen. Als je een even aantal mensen hebt, kunnen ze allemaal paren vormen en stilstaan (geïsoleerd). Maar als je een oneven aantal mensen hebt, blijft er één persoon over die geen partner heeft. Die ene persoon kan dan wel vrij rondrennen!
  • Dit betekent dat het gedrag van het systeem sterk afhangt van of je een even of oneven aantal deeltjes hebt.

4. Waarom is dit belangrijk?

Vroeger dachten wetenschappers dat ze alleen naar één deeltje hoefden te kijken om te begrijpen hoe een materiaal geleidt. Dit papier laat zien dat je naar het hele gezelschap (veel deeltjes tegelijk) moet kijken.

• De auteurs tonen aan dat hun nieuwe methode (het meten van de geometrische fase en de statistiek van de chaos) een zeer krachtige manier is om te zien of een materiaal een isolator of een geleider is, zelfs in complexe, onregelmatige situaties.
• Ze hebben ook laten zien dat het "de Moura-Lyra" model, dat eerder als raar en onbetrouwbaar werd beschouwd, eigenlijk heel interessant is omdat het deze vreemde "paar-dans" vertoont.

Samenvatting in één zin

De auteurs hebben een nieuwe "thermometer" bedacht om te meten of elektronen in een rommelig materiaal vastzitten of vrij rondzwerven, en ze ontdekten dat in sommige rommelige systemen, elektronen zich gedragen als dansers die alleen kunnen bewegen als er een "alleenstaande" in de groep is.

Technische samenvatting

Probleemstelling

Het onderscheid tussen geleiders en isolatoren in kwantumsystemen wordt fundamenteel bepaald door de localisatie van de ladingsdragers. Hoewel de moderne theorie van polarisatie (MPT) een krachtig raamwerk biedt om dit onderscheid in kristallijne systemen te maken door polarisatie te beschouwen als een geometrische fase (in plaats van een operatorverwachtingswaarde), is de toepassing van deze theorie op ongestoorde systemen met periodieke randvoorwaarden (PBC) problematisch.

In systemen met PBC kunnen polarisatie en zijn cumulanten niet direct worden berekend als operatorverwachtingswaarden. Bovendien is het moeilijk om de effecten van ongecorreleerde en gecorreleerde wanorde (disorder) kwantitatief te analyseren in de thermodynamische limiet zonder de gebruikte methoden aan te passen. Bestaande methoden zoals de inverse participatieverhouding of de localisatielengte zijn soms onvoldoende of moeilijk toe te passen op veeldeeltjessystemen. Het artikel richt zich op het ontwikkelen van technieken voor het middelen over wanorde binnen het kader van de MPT om localisatieovergangen nauwkeurig te detecteren.

Methodologie

De auteurs ontwikkelen een geavanceerde numerieke aanpak die de moderne theorie van polarisatie combineert met statistische momenten en cumulanten, specifiek aangepast voor systemen met periodieke randvoorwaarden en wanorde.

1. Geometrische Fase en Karakteristieke Functie:

   • In plaats van de polarisatie als operator te behandelen, wordt deze beschreven via de "polarisatie-amplitude" $Z_q = \langle \Psi | \exp(i \frac{2\pi q}{L} \hat{X}) | \Psi \rangle$, waarbij $\hat{X}$ de totale positie-operator is.
   • $Z_q$ fungeert als een discrete karakteristieke functie voor de waarschijnlijkheidsverdeling van de totale positie.

2. Middelen over Wanorde en Statistische Momenten:

   • De auteurs berekenen $Z_q$ voor vele realisaties van wanorde en middelen deze over energie-intervallen en wanorde-realisaties.
   • Om de schalingsinformatie te behouden en logaritmische afgeleiden te vermijden (die divergeren in metalen fasen), wordt de verdeling "gecentreerd" door de absolute waarde $|Z_q|$ te nemen.
   • Op basis hiervan worden benaderde tweede en vierde momenten ($\tilde{M}_2$ en $\tilde{M}_4$) gedefinieerd via eindige-differentie-derivaten.

3. Geometrische Binder Cumulant en Schalingsexponent:

   • De geometrische Binder cumulant ($U_4$) wordt afgeleid als een maat voor localisatie: $U_4 = 1 - \frac{1}{3} \frac{\tilde{M}_4}{\tilde{M}_2^2}$.
   • De grootte-schalingsexponent ($\gamma$) van de variantie wordt bepaald via de relatie $M_2 \propto L^\gamma$.
     • $\gamma = 2$ duidt op een gedelokaliseerde (metallische) toestand.
     • $\gamma \le 1$ (vaak $\approx 0$ of $1$) duidt op een gelokaliseerde (isolator) toestand.

4. Degeneratie-indicator:

   • Een nieuwe indicator wordt ontwikkeld die gebruikmaakt van de degeneratie van de grondtoestand in eindige systemen. Door een Peierls-fase ($\Phi = \pi/L$) toe te passen (wat overeenkomt met het wisselen tussen periodieke en antiperiodische randvoorwaarden), kunnen gelokaliseerde systemen worden onderscheiden van gedelokaliseerde systemen op basis van de verandering in $|Z_1|$. In gedelokaliseerde systemen kan dit leiden tot een tweevoudige degeneratie, terwijl gelokaliseerde systemen hier niet gevoelig voor zijn.

5. Toepassing op Modellen:

   • Anderson-model: Een 1D-model met ongecorreleerde wanorde.
   • de Moura-Lyra-model (dMLM): Een 1D-model met machtswet-gecorreleerde wanorde (gecontroleerd door parameter $\alpha$), bekend om zijn pathologische eigenschappen en de controverse rondom een "mobiliteitsrand" (mobility edge).
   • Berekeningen worden uitgevoerd voor zowel enkele-deeltjes toestanden als veeldeeltjes systemen (Slater-determinanten van fermionen), waarbij deeltjesdichtheid ($\rho$) een variabele is.

Belangrijkste Resultaten

1. Validatie met het Anderson-model:

   • De methode bevestigt dat het 1D Anderson-model volledig gelokaliseerd is voor alle wanordestrengths ($W > 0$).
   • Voor enkele-deeltjes toestanden daalt de schalingsexponent $\gamma$ sterk onder de eenheid (zelfs negatief in bepaalde intervallen), wat "hyperlocalisatie" aangeeft.
   • Voor veeldeeltjes systemen (met fermionen) is $\gamma \approx 1$ in de gelokaliseerde fase, wat overeenkomt met een eindige diëlektrische susceptibiliteit, en springt naar $\gamma = 2$ bij delokalisatie.

2. Analyse van het de Moura-Lyra-model (dMLM):

   • Globale Delokalisatie: De resultaten ondersteunen eerdere studies (Santos Pires et al.) die aangeven dat er een globale delokalisatie-overgang plaatsvindt bij $\alpha \approx 1$. Voor $\alpha > 1$ zijn alle toestanden gedelokaliseerd; er bestaat dus geen mobiliteitsrand (mobility edge) in de traditionele zin.
   • Specifiek Gedrag voor $\alpha > 2$: In het gebied $\alpha > 2$, nabij het bandcentrum (half-vulling), vertonen de energieniveaus een uniek gedrag: paren van niveaus sluiten hun gap (degenereren) en blijven gesloten.
   • Invloed van Deeltjesaantal: In veeldeeltjesberekeningen leidt deze degeneratie tot een sterk afhankelijkheid van het aantal deeltjes. De localisatiegevoelige grootheden (zoals $U_4$ en de degeneratie-indicator) bereiken hun maximale waarden (indicatie van maximale delokalisatie) alleen wanneer het aantal deeltjes oneven is (wat overeenkomt met een halfgevulde geleidingsband). Bij een even aantal deeltjes is het effect minder uitgesproken.

3. Fasediagram:

   • De auteurs stellen een nieuw fasediagram voor het dMLM op, gebaseerd op $\alpha$ en deeltjesdichtheid $\rho$. Dit diagram onderscheidt een gelokaliseerde fase, een volledig gedelokaliseerde fase, en een "intermediate delokaliseerde fase" rond $\alpha > 2$ en $\rho \approx 0.5$, gekenmerkt door de paar-degeneratie.

Bijdrage en Significantie

• Methodologische Vooruitgang: Het artikel biedt een robuust raamwerk om de moderne theorie van polarisatie toe te passen op wanordelijke systemen met periodieke randvoorwaarden. Door het gebruik van de geometrische Binder cumulant en de schalingsexponent $\gamma$ kunnen localisatie-overgangen worden gekwantificeerd zonder de noodzaak van het vergelijken van systemen met verschillende groottes op een complexe manier.
• Oplossing van een Controverse: De studie draagt bij aan het oplossen van de langdurige discussie over het bestaan van een mobiliteitsrand in het de Moura-Lyra-model. De resultaten tonen aan dat de schijnbare mobiliteitsrand in eerdere studies waarschijnlijk een artefact was van eindige systeemgroottes of specifieke meetmethoden, en dat in plaats daarvan een globale overgang optreedt bij $\alpha \approx 1$.
• Nieuw Inzicht in Correlatie: Het artikel benadrukt het belang van het onderscheid tussen enkele-deeltjes en veeldeeltjesberekeningen. Fermionische correlaties (Pauli-principe) spelen een cruciale rol in hoe wanorde localisatie beïnvloedt, en de nieuwe methoden maken het mogelijk om deze effecten expliciet te bestuderen.
• Toepasbaarheid: De ontwikkelde technieken (geometrische cumulanten, degeneratie-indicator) zijn van bredere waarde voor het bestuderen van andere disordered systemen, inclusief topologische isolatoren en systemen met veeldeeltjeslocalisatie (MBL).

Kortom, dit werk verenigt de theorie van polarisatie met statistische analyse van wanorde, levert een nauwkeurige diagnose van localisatie in 1D-systemen en verheldert de fysische eigenschappen van het complexe de Moura-Lyra-model.