Oddities in the Entanglement Scaling of the Quantum Six-Vertex Model

Probleemstelling

De auteurs onderzoeken de verstrengelingseigenschappen van het Quantum Six-Vertex Model (6VM) op een cilinder, met name in de limiet van de Rényi-entropie van orde $n = \infty$ (de min-entropie of Shannon-Rényi entropie). Hoewel het gedrag van verstrengeling bij kritieke systemen (beschreven door Conformale Veldtheorie of CFT) goed begrepen is, is er weinig bekend over de invloed van de pariteit van de systeemgrootte (of de lengte $L$ even of oneven is) op deze schaling in twee dimensionale systemen.

In eerdere studies aan 1D kritieke systemen bleek dat oneven lengtes tot oscillaties in verstrengeling kunnen leiden door geometrische frustratie. De centrale vraag in dit werk is: Heeft de pariteit van de diameter van de cilinder (oneven vs. even) een meetbaar effect op de verstrengelingsschaling van het 6VM, en zo ja, wat onthult dit over de onderliggende fysica?

Methodologie

De studie combineert analytische theorie, numerieke simulaties en het gebruik van oneindige Matrix Product States (iMPS).

1. Model en Koppeling:

   • Het 6VM wordt gedefinieerd op een vierkant rooster op een cilinder met omtrek $L$ en hoogte $2h$.
   • De grondtoestand wordt beschreven door een Rokhsar-Kivelson (RK) golf functie, een superpositie van alle microscopische configuraties gewogen door $e^{-E(c)/2}$.
   • Er wordt een equivalentie gebruikt tussen de Schmidt-eigenwaarden van de 2D RK-golf functie en de waarschijnlijkheidsverdeling van de grondtoestand van een 1D XXZ spin-1/2 keten met periodieke randvoorwaarden (PBC).
   • De berekening van de Rényi-verstrengeling $S_\infty$ van het 2D-systeem komt neer op het berekenen van de Shannon-Rényi entropie van de 1D XXZ-keten.

2. Analyse van Pariteit:

   • Even $L$: De grondtoestand heeft totale magnetisatie $S_z = 0$. De maximale toestand (Néel-toestand) is tweevoudig ontaard.
   • Oneven $L$: Door PBC kan niet elke spinpaar antiparallel zijn. Er moet een paar parallelle spins zijn, wat fungeert als een spinon (een $S_z = \pm 1/2$ excitatie). Dit leidt tot een extensieve ontaarding (groeiend met $L$) omdat het spinon zich op $L$ verschillende posities kan bevinden.

3. Technieken:

   • Exacte Diagonalisatie: Berekening van $S_\infty$ voor kleine systemen ($L$ tot 23) bij verschillende anisotropieparameters $\Delta$.
   • iMPS Ansatz: Gebruik van een oneindige Matrix Product State Ansatz gebaseerd op vertex-operatoren van een chirale bosonveld om de schaling voor grotere systemen te analyseren en de relatie met de Luttinger-parameter te leggen.
   • CFT Argumenten: Een heuristisch argument gebaseerd op de schaling van primaire velden in Conformale Veldtheorie om de oorsprong van de logaritmische correctie te verklaren.

Belangrijkste Bijdragen en Resultaten

1. Ontdekking van een Logaritmische Correctie:
   De auteurs ontdekken dat bij oneven systeemgroottes ($L$ oneven) er een extra logaritmische term ($\propto \log L$) optreedt in de schaling van de min-entropie $S_\infty$, die niet aanwezig is bij even $L$.
   De schalingsformule voor oneven $L$ is:
   $$S_\infty(L) \approx aL + b \log L + c$$
   Waarbij de coëfficiënt $b$ direct gerelateerd is aan de onderliggende CFT.

2. Relatie met de Luttinger-parameter en Compactificatieradius:
   De coëfficiënt $b$ van de logaritmische term is niet willekeurig, maar wordt bepaald door de compactificatieradius $R$ (of equivalent de Luttinger-parameter $K$) van de lage-energie bosonische veldtheorie die het kritieke punt beschrijft:
   $$b = \alpha = \pi R^2 = \frac{1}{4K}$$
   Hierbij is $\alpha$ de variatieparameter in de iMPS Ansatz. Dit betekent dat de verstrengelingsschaling direct informatie geeft over de spectrum van de onderliggende CFT.

3. Fysische Oorsprong (Spinon):
   De extra logaritmische term wordt toegeschreven aan de aanwezigheid van een spin-1/2 excitatie (spinon) in de grondtoestand bij oneven $L$. In CFT-termen correspondeert de oneven grondtoestand niet met het vacuüm, maar met een primaire toestand met conformale dimensie $\Delta_{1/2}$. De logaritmische correctie is een schalingsfactor die voortkomt uit deze extra excitatie.

4. Numerieke Bevestiging:

   • Bij specifieke punten ($\Delta = 0$ en $\Delta = 1/2$) werd de coëfficiënt $b$ numeriek bepaald als respectievelijk $1/4$ en $1/3$, wat perfect overeenkomt met de theoretische voorspellingen ($\Delta = -\cos(2\pi b)$).
   • De resultaten voor andere waarden van $\Delta$ in het kritieke gebied ($-1 < \Delta \leq 1$) tonen een goede overeenkomst tussen exacte diagonalisatie, iMPS-berekeningen en de theoretische curve.

5. Gedrag voor $n \neq \infty$:
   Numerieke simulaties tonen aan dat deze logaritmische correcties ook bestaan voor Rényi-entropieën met eindige $n$ (niet alleen $n=\infty$), hoewel de sterkte ervan afhangt van zowel $n$ als $\Delta$. Er wordt een discontinuïteit waargenomen bij de overgang naar de von Neumann-entropie ($n \to 1$).

Significantie en Conclusie

Dit werk biedt een nieuw inzicht in de relatie tussen geometrie, frustratie en kritieke verstrengeling. De belangrijkste conclusies zijn:

• Pariteit als Proef: De pariteit van de systeemgrootte is een cruciale factor die de verstrengelingsschaling beïnvloedt in 2D kritieke systemen. Een oneven grootte introduceert een fundamentele geometrische frustratie die leidt tot een extra logaritmische term.
• Toegang tot CFT-data: De coëfficiënt van deze logaritmische term fungeert als een directe "sonde" voor de onderliggende Conformale Veldtheorie. Het maakt het mogelijk om de Luttinger-parameter ($K$) en de compactificatieradius ($R$) af te leiden uit verstrengelingsmetingen, zonder dat men de volledige spectrum hoeft te kennen.
• Universeel Gedrag: Het effect is universeel voor kritieke RK-golf functies en is niet afhankelijk van de topologie van de variërende manifold of de gladheid ervan, maar puur van de interne spin-configuratie en de grootte-pariteit.

Samenvattend demonstreert het artikel dat de "onevenheid" van een systeem niet slechts een kleine correctie is, maar een fundamenteel mechanisme dat de verstrengelingsschaling herschrijft en directe toegang biedt tot de lage-energie veldtheoretische beschrijving van het materiaal.