Variational optimization of projected entangled-pair states on the triangular lattice

Technische Samenvatting: Variationele Optimalisatie van Projected Entangled-Pair States op het Triangulaire Rooster

Probleemstelling
De simulatie van gefrustreerde kwantumveeldeeltjessystemen, met name op triangulaire en kagome-roosters, blijft een grote uitdaging in de gecondenseerde materie fysica. Hoewel Projected Entangled-Pair States (PEPS) een standaardraamwerk zijn geworden voor het simuleren van tweedimensionale kwantumsystemen, is hun variationele optimalisatie op triangulaire geometrieën historisch gezien beperkt geweest. Conventionele benaderingen brengen het triangulaire rooster doorgaans in kaart op een vierkant rooster, wat diagonale interacties tussen naburige plaatsen (next-nearest-neighbor) introduceert. Deze mapping dwingt correlaties langs de diagonale richting om via de resterende roosterrichtingen te worden gerouteerd, wat een inefficiënte representatie is van de intrinsieke verstrengelingsstructuur van het systeem. Bovendien vertrouwen bestaande methoden vaak op indirecte mappings die er niet volledig in slagen de rotationele symmetrie en de specifieke connectiviteit van het triangulaire Bravais-rooster te vatten, wat de nauwkeurigheid van de variationele energiediminimalisatie bij een vaste bindingsdimensie kan beperken.

Methodologie
De auteurs introduceren een algemeen Corner Transfer Matrix Renormalization Group (CTMRG) algoritme dat specifiek is afgestemd op oneindige PEPS op het triangulaire rooster. De belangrijkste methodologische componenten zijn:

• Natuurlijke Triangulaire Formulering: In plaats van te mappen naar een vierkant rooster, opereert het algoritme direct op de triangulaire geometrie. De PEPS-ansatz wordt gedefinieerd met zeven-index tensoren (één fysieke index en zes virtuele indices die verbinding maken met alle naburige plaatsen), wat de coördinatiegetal van het triangulaire rooster reflecteert.
• Automatische Differentiatie (AD): De aanpak integreert AD om energiegradiënten efficiënt te berekenen, wat directe variationele energie-optimalisatie mogelijk maakt zonder de noodzaak van handmatige afleiding van complexe update-regels.
• Radiale CTMRG-procedure: In tegenstelling tot de vierkante rooster CTMRG, die orthogonale richtingen gebruikt, hanteert de triangulaire variant een radiaal absorptieschema. Het omgeving bestaat uit zes hoektensoren ($C$) en twaalf randtensoren ($T$). Het algoritme voert iteratieve machtsmethoden uit om fixed-point omgevingstensoren te berekenen.
• Absorptie en Renormalisatie: De absorptie van PEPS-tensoren in de omgeving wordt radiaal uitgevoerd. Vanwege de niet-orthogonale roostervectoren zijn updates aan hoek- en randtensoren onderling afhankelijk. De methode maakt gebruik van een tweestaps renormalisatieproces waarbij twee sets projectoren ($P$ en $Q$) worden gebruikt om de vergrote bindingsdimensies terug te schalen naar een vaste omgeving-bindingsdimensie $\chi_E$, wat numerieke stabiliteit garandeert.
• Symmetrie-afhandeling: Het algoritme gaat geen rotatiesymmetrie aan van de input-tensoren, waardoor het willekeurige toestanden kan accommoderen, inclus� including symmetrie-gebroken en spiraalvormige toestanden. Het kan ook worden uitgebreid om fysieke symmetrieën te integreren voor verbeterde efficiëntie.

Belangrijkste Bijdragen

1. Algoritmische Ontwikkeling: Het artikel presenteert de eerste computationeel efficiënte en nauwkeurige methode voor het contracteren en variationeel optimaliseren van PEPS direct op het triangulaire rooster, waarbij de noodzaak voor vierkante rooster-mappings wordt omzeild.
2. Verbeterde Representatiekracht: Door de natuurlijke geometrie te gebruiken, vangt de triangulaire PEPS (TL-PEPS) de kwantumcorrelaties langs drie onafhankelijke roosterrichtingen getrouwer op dan de vierkante rooster (SL-PEPS) tegenhangers. Theoretische analyse toont aan dat voor een gegeven bindingsdimensie $\chi_B$, de bovenste grens van de verstrengelingsentropie voor TL-PEPS aanzienlijk hoger is (schalend als $(8L-2)\log\chi_B$) vergeleken met SL-PEPS (schalend als $4L\log\chi_B$).
3. Toename van Variationele Parameters: Bij een vaste bindingsdimensie bezit de TL-PEPS ansatz een groter aantal variationele parameters ($p\chi_B^6$) vergeleken met de SL-PEPS ($p\chi_B^4$), wat meer expressieve kracht biedt voor dezelfde computationele middelen in termen van bindingsdimensie.
4. Benchmarking Framework: De auteurs bieden een rigoureus benchmarking-protocol dat native TL-PEPS vergelijkt met gemapte SL-PEPS, waarbij een eerlijke vergelijking wordt gewaarborgd door het matchen van de variationele parameters of het analyseren van convergentie met betrekking tot de omgeving-bindingsdimensies.

Resultaten
De methode werd toegepast op twee paradigmatische antiferromagnetische spin-1/2 Heisenberg-modellen:

• Triangulair Rooster Heisenberg Model:

  • De TL-PEPS benadering leverde lagere variationele energieën op dan de SL-PEPS bij dezelfde bindingsdimensies.
  • Bij $\chi_B = 5$ bereikte de TL-PEPS een grondtoestandsenergie van $E = -0.55136$, de laagste gerapporteerde variationele tensornetwerkwaarde tot nu toe voor dit model, wat verbetering biedt ten opzichte van eerdere geëxtrapoleerde DMRG en iPEPS resultaten.
  • Extrapolatie naar de oneindige bindingsdimensie-limiet leverde $E = -0.55164$ op, consistent met SL-resultaten maar verkregen met verbeterde convergentie-eigenschappen.
  • De magnetische ordeparameter $m$ vertoonde een systematische reductie in TL-PEPS ($m \approx 0.157$) vergeleken met SL-PEPS ($m \approx 0.182$), toegeschreven aan de verbeterde gevoeligheid van de ansatz voor de verstrengelingsstructuur.

• Kagome Rooster Heisenberg Model:

  • De TL-PEPS bood consequent lagere variationele energieën over alle toegankelijke bindingsdimensies vergeleken met de SL-formulering.
  • De geëxtrapoleerde grondtoestandsenergie was $E = -0.43929$, wat de eerdere DMRG en exact diagonalisatie resultaten licht heeft verbeterd.
  • De TL-PEPS data vertoonden een schoner schaalgedrag met betrekking tot de inverse correlatielengte ($1/\xi$), wat meer betrouwbare extrapolaties faciliteert, zelfs bij lagere bindingsdimensies waar SL-PEPS data onnauwkeurigheden vertoonden.
  • De resultaten ondersteunen een niet-magnetische grondtoestand die consistent is met een kwantum spin vloeistof (QSL), hoewel de vraag of deze een gap of gapless is, een open vraag blijft in de literatuur.

Significantie en Claims
Het artikel claimt dat de native triangulaire roosterformulering een significant voordeel biedt in variationele nauwkeurigheid voor gefrustreerde kwantummagneten. De verbeteringen worden toegeschreven aan:

1. Getrouwe Verstrengelingsrepresentatie: De ansatz respecteert van nature de roostergeometrie en rotationele symmetrieën, wat leidt tot een nauwkeurigere representatie van de verstrengelingsstructuur.
2. Parameter Efficiëntie: Het grotere aantal variationele parameters bij een vaste bindingsdimensie stelt de ansatz in staat om de Hilbertruimte effectiever te verkennen.
3. Haalbaarheid: Ondanks de hogere computationele kosten per stap (schalend als $O(\chi_B^{13})$ vergeleken met het vierkante rooster), blijft de methode haalbaar en competitief, vooral wanneer deze wordt gematcht tegen vierkante rooster-simulaties met een gelijk aantal variationele parameters.

De auteurs concluderen dat deze methode een onbevooroordeeld perspectief biedt op open problemen in de gefrustreerde kwantummagnetisme en de weg vrijmaakt voor het bestuderen van andere op driehoeken gebaseerde roosters (bijv. honingraat, dice, maple-leaf) en modellen met next-to-nearest-neighbor interacties zonder complexe routing. De implementatie is beschikbaar in de open-source variPEPS bibliotheek.