Tensor renormalization group approach to the $O(2)$ models via symmetry-twisted partition functions

Dit artikel toont aan dat de tensorrenormalisatiegroep-methode met symmetrie-gedraaide partitiefuncties effectief kan worden gebruikt om spontane symmetriebreking in drie dimensies en het BKT-overgangspunt in twee dimensies te detecteren, en bevestigt bovendien de succesvolle toepassing op het gegeneraliseerde tweedimensionale $O(2)$-model voor het identificeren van overgangen tussen ferromagnetische, nematische en paramagnetische fasen.

De kern

Stel je voor dat je een gigantisch, driedimensionaal tapijt van magneetjes bekijkt. Elk magneetje kan in elke richting wijzen, maar ze willen graag met hun buren meedoen: ofwel allemaal in dezelfde richting wijzen (geordend), ofwel willekeurig rondzwieren (chaotisch).

Fysici proberen te begrijpen hoe dit tapijt zich gedraagt als je het verwarmt of afkoelt. Op een bepaald punt verandert het gedrag van het hele tapijt drastisch: dit noemen we een fase-overgang. Het is alsof water plotseling van vloeistof naar ijs verandert, maar dan met magneetjes.

Het probleem is dat het berekenen van precies wanneer en hoe dit gebeurt, voor computers extreem moeilijk is. Traditionele methoden (zoals Monte Carlo-simulaties) lopen vaak vast in een wiskundige "val" (het zogenaamde tekenprobleem), waardoor ze geen goed beeld krijgen van wat er gebeurt.

De Oplossing: Een "Symmetrie-Verdraaide" Spelregel

In dit paper gebruiken de onderzoekers een slimme nieuwe truc: de Tensor Renormalization Group (TRG). Denk hierbij niet aan een simpele rekenmachine, maar aan een superkrachtige manier om het tapijt stap voor stap te vouwen en te samenvatten, zonder de belangrijke details te verliezen.

Maar de echte innovatie is hun gebruik van symmetrie-verdraaide partitiefuncties. Laten we dit uitleggen met een analogie:

• Het Normale Tapijt: Stel je voor dat je een dansvloer hebt waar iedereen vrij kan dansen. Als het koud is, dansen ze allemaal synchroon (geordend). Als het heet is, dansen ze wild en willekeurig (chaotisch).
• De Twist (De Verdraaiing): De onderzoekers doen alsof ze de wanden van de danszaal een beetje "draaien". Ze zeggen: "Als je de vloer rondloopt en terugkomt bij je startpunt, moet je je danspas een klein beetje hebben gedraaid."
  • Als de dansers chaotisch zijn (hoge temperatuur), maakt deze draaiing hen niet uit. Ze dansen toch al willekeurig, dus het resultaat is hetzelfde als zonder draaiing.
  • Als de dansers synchroon zijn (lage temperatuur), is die draaiing in de wanden een enorm probleem. De dansers kunnen niet meer perfect synchroon blijven zonder botsingen. Het systeem "voelt" de spanning.

Wat levert dit op?

1. Het 3D-Model (De Magneetjes in de Ruimte):
   In drie dimensies breken de magneetjes op een bepaalde temperatuur plotseling hun symmetrie en kiezen ze allemaal één richting. De onderzoekers gebruikten hun "verdraaide wanden" om precies te meten wanneer dit gebeurt. Ze vonden de exacte temperatuur en hoe snel de orde ontstaat. Het is alsof ze een heel precies thermometer hebben gevonden die aangeeft op welk gradenpunt water vries.

2. Het 2D-Model (Het Vlakke Tapijt):
   Hier is het nog interessanter. In twee dimensies kunnen de magneetjes nooit perfect synchroon worden (ze blijven trillen). Maar ze kunnen wel een "quasi-orde" bereiken, een soort vloeibare kristal-achtige staat. Dit heet de BKT-overgang.
   De "verdraaide wanden" meten hier direct de stijfheid van het tapijt (de helicity modulus).

   • Analogie: Stel je voor dat je een rubberen mat trekt. Als het koud is, is de mat erg strak en veert hij terug (hoge stijfheid). Als het heet is, is hij slap. De onderzoekers zagen precies het punt waarop de mat van "strak" naar "slap" overgaat, zonder dat ze de magneetjes hoeven te "forceren" met externe krachten.

3. Het Geavanceerde Model (De "Nematische" Fase):
   Ze keken ook naar een complexere versie waar de magneetjes niet alleen in één richting willen, maar ook in "gebroken" richtingen (zoals 180 graden omgedraaid, maar niet 90). Dit creëert een nieuwe, vreemde fase genaamd nematisch.
   Met hun truc konden ze twee verschillende overgangen zien:

   • Van chaos naar nematisch.
   • Van nematisch naar geordend.
     Het is alsof ze een bril hebben opgezet die twee verschillende kleuren van een regenboog kan onderscheiden, terwijl andere methoden alleen de hele regenboog als één vage kleur zien.

Waarom is dit belangrijk?

• Geen "Valse" Signalen: Traditionele methoden moeten vaak een externe kracht uitoefenen om de fase-overgang te zien, wat de resultaten kan verstoren. Deze methode kijkt puur naar de interne eigenschappen van het systeem.
• Precisie: Ze kregen veel nauwkeurigere resultaten dan eerder mogelijk was, zelfs voor de moeilijkste gevallen.
• Toekomst: Het bewijst dat deze "tensor-netwerk" methode een krachtig alternatief is voor de oude methoden, vooral voor problemen waar computers normaal op vastlopen.

Kortom:
De onderzoekers hebben een slimme manier bedacht om naar een complex systeem te kijken door het een beetje te "verdraaien". Door te kijken hoe het systeem reageert op deze draaiing, kunnen ze precies zien wanneer en hoe het van chaos naar orde (of naar een andere vreemde toestand) springt. Het is als het vinden van de perfecte temperatuur om een ijsblokje te laten smelten, alleen dan met magneetjes in een computerwereld.

Technische samenvatting

Titel: Tensor-renormalisatiegroepbenadering voor O(2)-modellen via symmetrie-gedraaide partitiefuncties

Auteurs: Shinichiro Akiyama et al.
Context: 42e Internationale Symposium voor Lattice Field Theory (LATTICE2025).

---

1. Het Probleem

Het bestuderen van kritische fenomenen en faseovergangen in roosterveldtheorieën, met name bij systemen met continue symmetrieën zoals het O(2)-model, vormt een uitdaging voor traditionele numerieke methoden.

• Beperkingen van Monte Carlo (MC): MC-methode lijdt vaak aan het "sign-probleem" bij complexe acties en heeft moeite met het direct berekenen van partitiefuncties en de grondtoestandsdegeneratie.
• Detectie van symmetriebreking: Het identificeren van spontane breking van continue symmetrieën (zoals U(1)) is lastig met conventionele ordeparameters, vooral in dimensies waar geen lange-afstandsorde optreedt (zoals in 2D door het Mermin-Wagner-theorema, hoewel er een BKT-overgang is).
• Gu-Wen ratio: Hoewel de Gu-Wen ratio (een verhouding van partitiefuncties) nuttig is voor discrete symmetrieën, is deze minder effectief voor continue symmetrieën of voor het onderscheiden van complexe fasen zoals nematiciteit.

2. Methodologie

De auteurs gebruiken de Tensor Renormalization Group (TRG) methode, een numerieke real-spaar renormalisatiegroepbenadering die geschikt is voor systemen met een sign-probleem en directe toegang biedt tot partitiefuncties.

• Symmetrie-gedraaide partitiefuncties ($Z_{g_0}$):
  In plaats van een externe veld te introduceren, wordt een "twist" (draaiing) in de randvoorwaarden geïmplementeerd. Dit wordt gedaan door een achtergrond-eenveld $A$ in te voeren, gekenmerkt door holonomieën $(g_1, \dots, g_d)$ langs de cyclus van de torus.
  De gedraaide partitiefunctie wordt gedefinieerd als $Z_{g_0} = Z[A=(1, \dots, 1, g_0)]$.

  • In de thermodynamische limiet geeft de verhouding $Z_{g_0}/Z_1$ informatie over de symmetrie-realiserende fase:
    • Symmetrische fase: $Z_{g_0}/Z_1 = 1$ voor alle $g_0$.
    • Symmetrie-gebroken fase: $Z_{g_0}/Z_1 = 0$ als $g_0$ niet in de resterende subgroep $H$ zit, en $1$ als $g_0 \in H$.
  • Deze discontinuïteit fungeert als een scherpe ordeparameter voor de faseovergang.

• Implementatie in TRG:

  • De auteurs construeren fundamentele tensoren met een eindige bond-dimensie (via een afgekapt karakter-expansie) die de globale symmetrie direct respecteren.
  • Ze gebruiken Symmetry Blocking om de twist in het TRG-algoritme te integreren.
  • Voor 3D wordt Anisotropic TRG (ATRG) gebruikt; voor 2D wordt Bond-weighted TRG (BTRG) gebruikt.
  • De twist wordt geparametriseerd als $g_0 = e^{i\alpha}$, waarbij $\alpha$ de draaiingshoek is.

• Modellen onderzocht:

  1. 3D O(2)-model: Spontane breking van U(1) symmetrie.
  2. 2D O(2)-model: Berezinskii-Kosterlitz-Thouless (BKT) overgang zonder spontane symmetriebreking.
  3. 2D Generalized O(2)-model: Een Hamiltoniaan met ferromagnetische en nematische termen (gecontroleerd door parameter $\Delta$ en integer $q$), wat leidt tot fasen: paramagnetisch, nematic en ferromagnetisch.

3. Belangrijkste Bijdragen en Resultaten

A. 3D O(2)-model

• Doel: Bepalen van het kritieke punt ($T_c$) en de kritieke exponent $\nu$.
• Methode: Analyse van de verhouding $Z_\pi/Z_0$ (met twist $\alpha=\pi$) voor verschillende systeemgroottes $L$.
• Resultaat:
  • Er wordt een schaal-invariant kruispunt waargenomen dat het kritieke punt markeert.
  • Via finite-size scaling analyse wordt $T_c = 2.2017(2)$ bepaald.
  • De kritieke exponent wordt bepaald als $\nu = 0.663(33)$. Dit is een van de eerste precieze TRG-bepalingen van $\nu$ voor dit model.

B. 2D O(2)-model (BKT-overgang)

• Doel: Bepalen van de BKT-overgangstemperatuur ($T_{BKT}$).
• Innovatie: De gedraaide partitiefunctie levert direct de helicity modulus ($\Upsilon_\alpha$) op via de relatie:
  $$\Upsilon_\alpha = -\frac{T}{\alpha^2} \frac{L_2}{L_1} \log \left( \frac{Z_\alpha}{Z_0} \right)$$
• Methode: Toepassing van het Nelson-Kosterlitz (NK) criterium. In de BKT-fase vertoont de helicity modulus een "universele sprong" naar $2T/\pi$ bij de overgang.
• Resultaat:
  • Door de kruising van $\Upsilon_{\pi/6}$ met de NK-curve te extrapoleren naar $L \to \infty$, wordt $T_{BKT} = 0.8927(5)$ gevonden.
  • Dit resultaat komt uitstekend overeen met eerdere studies, maar wordt bereikt zonder externe velden.

C. 2D Generalized O(2)-model ($q=2$)

• Doel: Identificeren van overgangen tussen paramagnetisch, nematic en ferromagnetische fasen.
• Methode: Gebruik van verschillende twist-hoeken $\alpha$:
  • $\alpha \in \pi\mathbb{Z}$ (bijv. $\pi$) detecteert de breking van de discrete $\mathbb{Z}_2$ symmetrie (ferromagnetisch vs. nematic).
  • $\alpha \notin \pi\mathbb{Z}$ (bijv. $\pi/6$) detecteert de BKT-overgang (nematic vs. paramagnetisch).
• Resultaten:
  • Ferromagnetisch-Nematic overgang: Een sprong in $Z_\pi/Z_0$ van 0 naar 1 duidt op een tweede-orde overgang (Ising-universaliteit). $T_c = 0.35253(7)$, wat nauwkeuriger is dan eerdere MC-schattingen op basis van specifieke warmte.
  • Nematic-Paramagnetisch overgang: De helicity modulus toont een universele sprong, maar aangepast voor half-vortex excitaties ($q=2$). De NK-criterium wordt aangepast naar $2q^2T/\pi$.
  • $T_{BKT}$ wordt bepaald als $0.7581(4)$.
• Voordeel: Deze methode vereist geen expliciet symmetrie-brekend veld, wat de rekenkosten verlaagt en de nauwkeurigheid verhoogt ten opzichte van traditionele TRG-benaderingen die magnetische susceptibiliteit analyseren.

4. Significantie en Toekomstperspectief

• Universeel Toepasbaar: De paper demonstreert dat symmetrie-gedraaide partitiefuncties een krachtige, universele ordeparameter zijn voor zowel discrete als continue symmetriebreking, en voor topologische overgangen zoals BKT.
• Efficiëntie: De methode elimineert de noodzaak van externe velden, wat de numerieke stabiliteit en efficiëntie van TRG-berekeningen aanzienlijk verbetert.
• Nieuwe inzichten: Het biedt een directe route om de helicity modulus en de renormalisatie van de compactificatieradius te berekenen in roostertheorieën.
• Toekomst: De auteurs plannen om deze techniek uit te breiden naar generalisaties met $q > 4$ (waarbij ongebruikelijke ferromagnetische fasen worden verwacht) en naar 3D generalisaties van het O(2)-model, gebieden die tot nu toe weinig zijn onderzocht met TRG-methoden.

Conclusie:
Dit werk markeert een belangrijke stap in de toepassing van Tensor Networks op kritische fenomenen. Door de combinatie van TRG met symmetrie-gedraaide randvoorwaarden, slagen de auteurs erin om faseovergangen in O(2)-modellen met hoge precisie te karakteriseren, inclusief het kwantificeren van kritieke exponenten en het lokaliseren van BKT-overgangen, zonder de beperkingen van conventionele Monte Carlo-methoden.