Spin subdiffusion in perturbed infinite-U Hubbard chain

Titel: Spin-subdiffusie in een verstoord oneindig-U Hubbard-keten

Auteurs: J. Rękas, M. Mierzejewski, Z. Lenarčič, en P. Prelovšek

---

1. Probleemstelling en Context

Het artikel onderzoekt het transport van spin en lading in het één-dimensionale $t$-model, wat de limiet $U \to \infty$ van het Hubbard-model voorstelt. In dit regime zijn dubbele bezettingen van roosterplekken verboden, wat leidt tot een "gefrragmenteerde" Hilbertruimte.

• Integrabiliteit: Het ongestoorde $t$-model ($t_\uparrow = t_\downarrow$) is integraal. De ladingsdynamica gedraagt zich als die van vrije spinloze fermionen, maar de spinvolgorde is "bevroren" (behalve bij translatie).
• Uitdaging: Hoewel de spinvolgorde gefragmenteerd is, is de transportdynamica van spin niet triviaal. De vraag is hoe spintransport zich gedraagt in het integrabele geval versus het geval waar de integrabiliteit wordt verbroken (door verstoringen zoals $t_\uparrow \neq t_\downarrow$ of $J_z \neq 0$), terwijl de Hilbertruimte-fragmentatie behouden blijft.
• Specifiek doel: Het begrijpen van de overgang van ballistisch transport naar diffusie en subdiffusie, en het ontrafelen van het mechanisme achter spin-subdiffusie in een systeem zonder wanorde of dipoolbehoud.

2. Methodologie

De auteurs combineren analytische afleidingen met uitgebreide numerieke simulaties:

• Modellen:
  • Het $t$-model (integrabel).
  • Het verstoord $t$-model met spin-afhankelijke hopping ($t_\uparrow \neq t_\downarrow$).
  • Het $t-J_z$-model (met spin-spin interactie zonder uitwisseling).
• Analytische technieken:
  • Thouless-argument: Gebruik van spin-afhankelijke flux (twisted boundary conditions) om de gevoeligheid van energieniveaus te meten en zo de spinstroom te bepalen.
  • Mazur-ongelijkheid: Om de spinstijfheid te relateren aan de overlap tussen spin- en ladingsstromen.
  • Lineaire respons theorie: Berekening van de dynamische spin-diffusieconstante $D_s(\omega)$.
• Numerieke methoden:
  • Microcanonical Lanczos Method (MCLM): Voor het berekenen van dynamische correlatiefuncties en diffusie in gesloten systemen (tot $L=20$ sites).
  • Open systeem dynamica: Oplossing van de Lindblad-masterequatie voor systemen gekoppeld aan rand-baden om stationaire stromen te analyseren.

3. Belangrijkste Resultaten

A. Koppeling tussen Lading en Spin

Een fundamentele bevinding is dat spintransport in dit model intrinsiek gekoppeld is aan ladingstransport vanwege de Hilbertruimte-fragmentatie.

• Bij magnetisatie $m=0$ (in een canoniek sector) is de diagonale matrixelement van de spinstroom nul. Er is geen stationaire spinstroom mogelijk, zelfs niet in een open systeem, tenzij er ook ladingsuitwisseling mogelijk is.
• Bij $m \neq 0$ is er een niet-nul overlap tussen de spinstroom ($j_s$) en de behouden ladingsstroom ($j_c$) in het integrabele geval.

B. Transport in het Integrabele $t$-model

• Bij $m \neq 0$: Het transport is ballistisch (beperkt door de behoudswet van de ladingsstroom).
• Bij $m = 0$ (Grand Canonicaal Gemiddelde): Hoewel elke individuele sector met $m=0$ geen stroom heeft, leidt het gemiddelde over alle magnetisatiesectoren tot een dissipatieloze diffusie. De diffusieconstante $D_s$ is eindig en schaalt met de temperatuur en deeltjesdichtheid, vergelijkbaar met het easy-axis XXZ-model.

C. Transport in het Verstoord (Niet-Integrabel) $t$-model

Wanneer de integrabiliteit wordt verbroken (bijv. door $t_\uparrow \neq t_\downarrow$ of $J_z \neq 0$), maar de fragmentatie behouden blijft:

• Bij $m \neq 0$: Het transport verandert van ballistisch naar normale diffusie. De diffusieconstante schaalt als $D_s \propto m^2 / \Delta t^2$.
• Bij $m = 0$ (Grand Canonicaal Gemiddelde): Dit is de meest opvallende bevinding. Het systeem vertoont subdiffusie.
  • De diffusieconstante $D_s(\omega \to 0)$ gaat naar nul in de thermodynamische limiet ($L \to \infty$).
  • De geïntegreerde diffusie $I_s(\omega)$ vertoont een schaling die afwijkt van normaal diffuus gedrag.

D. Het Mechanisme van Subdiffusie

De auteurs identificeren het mechanisme voor subdiffusie als zijnde verschillend van dat in wanordelijke systemen of systemen met dipoolbehoud.

• Omdat de ladingsdiffusieconstante $D_c$ evenredig is met $m^2$ (de lokale magnetisatie), wordt de effectieve diffusiecoëfficiënt voor spinexcitaties lokaal evenredig met het kwadraat van de magnetisatie: $D_{eff}(x) \propto u(x)^2$.
• Dit leidt tot een niet-lineaire subdiffusie-vergelijking (vergelijkbaar met de vergelijking van het poreuze medium):
  $$\frac{\partial u}{\partial \tau} \propto \frac{\partial}{\partial x} \left( u^2 \frac{\partial u}{\partial x} \right)$$
• De oplossing hiervoor (Barenblatt-oplossing) voorspelt een spreiding van het profiel die schaalt als $\tau^{1/4}$, wat overeenkomt met een exponent $\gamma = 3/4$ in de frequentieafhankelijkheid van de correlatiefunctie ($C_l(\omega) \propto \omega^{-3/4}$). Numerieke resultaten bevestigen deze exponent.

4. Bijdragen en Significantie

1. Uniek Transportregime: Het artikel toont aan dat subdiffusie kan optreden in zuivere, translatie-invariante kwantumsystemen zonder wanorde of expliciete dipoolbehoudswetten, puur als gevolg van Hilbertruimte-fragmentatie en de koppeling tussen lading en spin.
2. Mechanisme-ontleding: Het biedt een fysiek inzichtelijk mechanisme voor subdiffusie via de "poreuze medium" dynamiek, waarbij de diffusie wordt vertraagd door de lokale magnetisatie die zelf weer diffuus is.
3. Verband met XXZ-modellen: Het bevestigt en verduidelijkt eerdere voorspellingen over subdiffusie in het "folded XXZ"-model, maar biedt een meer intuïtief beeld vanuit het Hubbard-model (oneindige afstoting).
4. Rol van Randvoorwaarden: Het onderzoek benadrukt het cruciale verschil tussen canonieke sectoren (waar spinstroom bij $m=0$ strikt nul is) en grand canonische gemiddelden (waar diffusie optreedt), en hoe open systemen met specifieke randdissipatie dit gedrag kunnen bevestigen of veranderen.

Conclusie:
De studie concludeert dat in het verstoord $t$-model de spintransportdynamica wordt gedicteerd door de fragmentatie van de Hilbertruimte. Dit leidt tot een unieke vorm van subdiffusie in de thermodynamische limiet, gedreven door de niet-lineaire koppeling tussen de diffusie van lading en de lokale spinpolarisatie. Dit onderscheidt het fundamenteel van andere bekende mechanismen voor anomal transport.