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

Riassunto Tecnico: Ottimizzazione Variazionale di Stati di Coppie Entangled Proiettate su Reticolo Triangolare

Enunciato del Problema
La simulazione di sistemi molti-corpo quantistici frustrati, in particolare su reticoli triangolari e kagome, rimane una sfida significativa nella fisica della materia condensata. Sebbene gli Stati di Coppie Entangled Proiettate (PEPS) siano diventati un framework standard per simulare sistemi quantistici bidimensionali, la loro ottimizzazione variazionale su geometrie triangolari è stata storicamente limitata. Gli approcci convenzionali tipicamente mappano il reticolo triangolare su un reticolo quadrato, introducendo interazioni di secondo vicino diagonali. Questa mappatura costringe le correlazioni lungo la direzione diagonale a essere instradate attraverso le restanti direzioni del reticolo, il che rappresenta una rappresentazione inefficiente della struttura intrinseca dell'entanglement del sistema. Inoltre, i metodi esistenti spesso si affidano a mappature indirette che non riescono a catturare pienamente la simmetria rotazionale e la specifica connettività del reticolo di Bravais triangolare, potenzialmente limitando l'accuratezza delle minimizzazioni dell'energia variazionale a dimensione di legame fissata.

Metodologia
Gli autori introducono un algoritmo generale di Gruppo di Renormalizzazione della Matrice di Trasferimento d'Angolo (CTMRG) specificamente adattato per PEPS infiniti su reticolo triangolare. I componenti metodologici chiave includono:

• Formulazione Triangolare Nativa: Invece di mappare su un reticolo quadrato, l'algoritmo opera direttamente sulla geometria triangolare. L'ansatz PEPS è definito utilizzando tensori a sette indici (un indice fisico e sei indici virtuali che collegano tutti i vicini prossimi), riflettendo il numero di coordinazione del reticolo triangolare.
• Differenziazione Automatica (AD): L'approccio integra l'AD per calcolare i gradienti di energia in modo efficiente, consentendo la minimizzazione diretta dell'energia variazionale senza la necessità di derivare manualmente regole di aggiornamento complesse.
• Procedura CTMRG Radiale: A differenza del CTMRG su reticolo quadrato, che utilizza direzioni ortogonali, la variante triangolare impiega uno schema di assorbimento radiale. L'ambiente consiste in sei tensori d'angolo ($C$) e dodici tensori di bordo ($T$). L'algoritmo esegue metodi di potenza iterativi per calcolare i tensori dell'ambiente a punto fisso.
• Assorbimento e Renormalizzazione: L'assorbimento dei tensori PEPS nell'ambiente è eseguito radialmente. A causa dei vettori del reticolo non ortogonali, gli aggiornamenti dei tensori d'angolo e di bordo sono interdipendenti. Il metodo utilizza un processo di rinormalizzazione a due fasi che coinvolge due set di proiettori ($P$ e $Q$) per troncare le dimensioni di legame ingrandite riportandole a una dimensione di legame dell'ambiente fissa $\chi_E$, garantendo la stabilità numerica.
• Gestione della Simmetria: L'algoritmo non assume la simmetria rotazionale dei tensori di input, permettendo di accomodare stati arbitrari, inclusi stati a rottura di simmetria e a spirale. Può anche essere esteso per incorporare simmetrie fisiche per migliorare l'efficienza.

Contributi Chiave

1. Sviluppo Algoritmico: Il documento presenta il primo metodo computazionalmente efficiente e accurato per contrarre e ottimizzare variazionalmente i PEPS direttamente su reticolo triangolare, bypassando la necessità di mappature su reticolo quadrato.
2. Rappresentabilità Migliorata: Utilizzando la geometria nativa, l'ansatz TL-PEPS cattura le correlazioni quantistiche lungo tre direzioni indipendenti del reticolo in modo più fedele rispetto ai corrispondenti SL-PEPS. L'analisi teorica mostra che, per una data dimensione di legame $\chi_B$, il limite superiore dell'entropia di entanglement per TL-PEPS è significativamente più alto (scalando come $(8L-2)\log\chi_B$) rispetto a SL-PEPS (che scala come $4L\log\chi_B$).
3. Aumento dei Parametri Variazionali: A una dimensione di legame fissata, l'ansatz TL-PEPS possiede un numero maggiore di parametri variazionali ($p\chi_B^6$) rispetto al SL-PEPS ($p\chi_B^4$), offrendo una maggiore capacità espressiva per le stesse risorse computazionali in termini di dimensione di legame.
4. Framework di Benchmarking: Gli autori forniscono un protocollo di benchmarking rigoroso che confronta il TL-PEPS nativo con il SL-PEPS mappato, garantendo confronti equi accoppiando i parametri variazionali o analizzando la convergenza rispetto alle dimensioni di legame dell'ambiente.

Risultati
Il metodo è stato applicato a due modelli paradigmatici di Heisenberg spin-1/2 antiferromagnetico:

• Modello di Heisenberg su Reticolo Triangolare:

  • L'approccio TL-PEPS ha prodotto energie variazionali inferiori rispetto al SL-PEPS alle stesse dimensioni di legame.
  • Con $\chi_B = 5$, il TL-PEPS ha raggiunto un'energia fondamentale di $E = -0.55136$, che è il valore di energia variazionale di rete tensoriale più basso riportato finora per questo modello, migliorando i precedenti risultati di estrapolazione DMRG e iPEPS.
  • L'estrapolazione al limite di dimensione di legame infinita ha fornito $E = -0.55164$, coerente con i risultati SL ma ottenuto con proprietà di convergenza migliorate.
  • Il parametro d'ordine magnetico $m$ ha mostrato una riduzione sistematica in TL-PEPS ($m \approx 0.157$) rispetto a SL-PEPS ($m \approx 0.182$), attribuita alla migliore sensibilità dell'ansatz alla struttura dell'entanglement.

• Modello di Heisenberg su Reticolo Kagome:

  • Il TL-PEPS ha costantemente fornito energie variazionali inferiori in tutte le dimensioni di legame accessibili rispetto alla formulazione SL.
  • L'energia fondamentale estrapolata è $E = -0.43929$, migliorando leggermente i precedenti risultati di DMRG ed estrazione esatta (ED).
  • I dati del TL-PEPS hanno mostrato un comportamento di scaling più pulito rispetto alla lunghezza di correlazione inversa ($1/\xi$), facilitando estrapolazioni più affidabili anche a dimensioni di legame inferiori dove i dati SL-PEPS mostravano imprecisioni.
  • I risultati supportano uno stato fondamentale non magnetico coerente con un liquido di spin quantistico (QSL), sebbene la natura gapless o gapped rimanga una questione aperta nella letteratura.

Significatività e Rivendicazioni
Il documento sostiene che la formulazione nativa del reticolo triangolare offre un vantaggio significativo nell'accuratezza variazionale per i magneti quantistici frustrati. I miglioramenti sono attribuiti a:

1. Rappresentazione Fedele dell'Entanglement: L'ansatz rispetta naturalmente la geometria del reticolo e le simmetrie rotazionali, portando a una rappresentazione più accurata della struttura dell'entanglement.
2. Efficienza dei Parametri: Il maggior numero di parametri variazionali a una dimensione di legame fissata consente all'ansatz di esplorare lo spazio di Hilbert in modo più efficace.
3. Fattibilità: Nonostante il costo computazionale più elevato per passo (che scala come $O(\chi_B^{13})$ rispetto al reticolo quadrato), il metodo rimane fattibile e competitivo, in particolare quando messo a confronto con simulazioni su reticolo quadrato con equivalenti conteggi di parametri variazionali.

Gli autori concludono che questo metodo fornisce una prospettiva imparziale sui problemi aperti nel magnetismo quantistico frustrato e apre la strada allo studio di altri reticoli basati sul triangolo (es. esagonale, dice, maple-leaf) e modelli con interazioni di secondo vicino senza complessi instradamenti. L'implementazione è disponibile nella libreria open-source variPEPS.