Tensor-network variational diagonalization of quantum many-body spectra

Technische Samenvatting: Tensor-Netwerk Variatonele Diagonalisatie van Kwantum-Vele-Deeltjes Spectra

Probleemstelling
De volledige eigenstructuur van een vele-deeltjes Hamiltoniaan bevat essentiële fysieke eigenschappen, waaronder thermodynamica, dynamische respons en de handtekeningen van kwantumchaos of lokalisatie. Echter, de exponentiële groei van de Hilbertruimte ($2^N$ voor $N$ spins) maakt exacte diagonalisatie (ED) onhandelbaar voor systemen groter dan ongeveer 30 spins, zelfs met het benutten van symmetrieën. Bestaande tensor-netwerkmethoden, zoals de Density-Matrix Renormalization Group (DMRG) en geprikkelde Matrix Product State (MPS) algoritmen, richten zich succesvol op specifieke spectrale sectoren (bijv. lage energie of gelokaliseerde toestanden), maar behandelen het energiespectrum doorgaans als een impliciete lijst in plaats van een opgeslagen variationeel object. De centrale uitdaging die in dit werk wordt aangepakt, is of de gehele exponentieel grote lijst van eigenenergieën gerepresenteerd en gesampled kan worden zonder expliciete enumeratie, mits het spectrum een comprimeerbare structuur bezit onder een passende etikettering.

Methodologie: Tensor-Netwerk Variatonele Diagonalisatie (TNVD)
De auteurs introduceren Tensor-Netwerk Variatonele Diagonalisatie (TNVD), een raamwerk dat twee gekoppelde tensor-netwerkobjecten gezamenlijk optimaliseert om de eigenwaardedecompositie $\hat{H} = \sum_{\alpha} E_\alpha |\alpha\rangle\langle\alpha|$ te representeren:

1. Spectrum MPS: De eigenenergieën $\{E_\alpha\}$ worden behandeld als een $N$-de orde tensor $E_{r_1 \dots r_N}$, waarbij de index $\alpha$ wordt gerepresenteerd door een binaire string $r = (r_1, \dots, r_N)$. Deze tensor wordt gecodeerd als een Matrix Product State (MPS). De virtuele bindingsdimensie $\chi_a$ van deze MPS dient als een directe sonde naar de "virtuele verstrengeling" van het energiespectrum in de gekozen labelruimte.
2. Variational Quantum Circuit (VQC): Een eindige diepte kwantumcircuit $\hat{U}$ brengt binaire producttoestanden $|r_\alpha\rangle$ in kaart naar benaderende eigentoestanden $|\alpha\rangle = \hat{U}|r_\alpha\rangle$. Het circuit bestaat uit een brick-wall architectuur van twee-deeltjes poorten met diepte $N_L$.
3. Gezamenlijke Optimalisatie: De binaire labels worden niet vastgesteld door energievolgorde, maar worden geleerd. De lokale tensoren van de spectrum MPS en de parameters van de VQC-poorten worden simultaan geoptimaliseerd om de logaritmische Hilbert-Schmidt loss te minimaliseren:
   $$F = \log_2 \|\hat{H} - \tilde{H}\|_{HS}^2 - N$$
   waarbij $\tilde{H}$ de ansatz Hamiltoniaan is. Deze globale loss-functie beperkt de gehele operator in plaats van een specifieke energieruimte te bevoordelen.
4. Computationele Kosten: De evaluatie van de loss omvat het contracteren van de spectrum MPS met de VQC en diens conjugaat, gemedieerd door super-diagonale tensoren. De kosten schalen polynomiaal met de systeemgrootte $N$ wanneer de MPS-dimensie $\chi_a$, de circuitdiepte $N_L$ en de TEBD-truncatiedimensie $\chi_t$ constant blijven.

Belangrijkste Resultaten

• Volledige Spectrale Reconstructie in de Ising-keten:
  Voor de transversale-veld Ising-keten reproduceert TNVD de exacte diagonalisatie-resultaten voor $N \le 16$. Cruciaal is dat voor $N=100$, de methode de dichtheid van toestanden (density of states) samplet uit een MPS die $2^{100}$ niveaus codeert zonder deze expliciet te enumereren. De gesamplede distributie komt overeen met Gaussische fits afgeleid van ED bij kleinere formaten.

  • Label Organisatie: Een controle-experiment met een willekeurige permutatie van binaire labels op het exacte spectrum resulteerde in een hoge virtuele verstrengelingsentropie, terwijl de door TNVD geoptimaliseerde labels een laag-verstrengelde, compacte MPS-representatie opleverden. Dit demonstreert dat de comprimeerbaarheid van de energietensor kritisch afhankelijk is van de geleerde organisatie van de labels.

• Gedisordreerde Systemen: Ising vs. XXZ:
  De auteurs vergeleken de willekeurige-veld Ising-keten met de willekeurige-veld XXZ-keten. Beide systemen vertonen vergelijkbare crossovers in niveau-spacing statistieken (van Poisson naar Gaussian Orthogonal Ensemble en terug naar Poisson) naarmate de wanorde (disorder) toeneemt, wat wijst op vergelijkbare spectrale chaos-eigenschappen.

  • Divergerende Fouten: Ondanks vergelijkbare niveau-statistieken waren de TNVD-fouten (Hilbert-Schmidt loss en gemiddelde energie-fout) aanzienlijk groter voor de XXZ-keten dan voor de Ising-keten onder identieke variationele middelen.
  • Oorzaakanalyse: De auteurs identificeerden dat spectrale chaos alleen niet de TNVD-complexiteit bepaalt. De bottleneck voor de XXZ-keten werd herleid tot de eindige-rang truncatie van de diagonaliserende transformatie.
    • Bij nul wanorde ($W=0$) vertoont de XXZ-keten (een $c=1$ kritische theorie) een veel tragere afname van de Schmidt-coëfficiënten en een zwaardere grondtoestand-Schmidt-staart vergeleken met de Ising-keten ($c=1/2$).
    • Consequent dragen de laag-energie eigentoestanden van de XXZ-keten substantieel meer Schmidt-gewicht buiten de truncatie-cutoff $\chi_t$ dan de Ising-keten, wat leidt tot hogere fouten. De Ising-keten, die bij $W=0$ een vrije-fermion integrabele theorie is, is beter geschikt voor laag-rang truncatie.

Betekenis en Claims
Het artikel claimt dat TNVD een directe test biedt voor wanneer volledige vele-deeltjes spectra een hanteerbare tensor-netwerk representatie toelaten. De methode verschuift het paradigma van volledige spectrale diagonalisatie van een enumeratie-taak naar een test van leerbare tensorstructuur.

De auteurs identificeren twee onafhankelijke structurele vereisten voor succesvolle TNVD:

1. Virtuele Verstrengeling in de Labelruimte: De energietensor moet een lage virtuele verstrengeling bezitten onder een geleerde binaire etikettering (gedemonstreerd door de Ising-keten).
2. Eindige-Rang Schmidt Truncatie: De diagonaliserende transformatie (de door het circuit geëvolueerde MPO) moet binnen de mogelijkheden van truncatie blijven, bepaald door het verval van de Schmidt-spectra (gedemonstreerd door het verschil tussen Ising en XXZ).

Het werk concludeert dat hoewel spectrale chaos (niveau-repulsie) een kenmerk van het spectrum is, de moeilijkheid van het representeren van het volledige spectrum via TNVD wordt bepaald door de interactie tussen de comprimeerbaarheid van de energietensor en de verstrengelingsstructuur van de eigentoestanden zelf. De broncode voor de implementatie is publiekelijk beschikbaar gesteld.