Site Basis Excitation Ansatz for Matrix Product States

Titolo: Site Basis Excitation Ansatz for Matrix Product States (SBEA)

Autore: Steven R. White (UC Irvine)
Oggetto: Un nuovo approccio efficiente per calcolare gli spettri di eccitazione elementare in sistemi quantistici unidimensionali utilizzando gli Stati a Prodotto di Matrice (MPS).

---

1. Il Problema

Lo spettro delle eccitazioni elementari è fondamentale per la comprensione teorica e il confronto con gli esperimenti, ma è computazionalmente più difficile da ottenere rispetto alle proprietà statiche.

• Limiti dei metodi esistenti: I metodi DMRG (Density Matrix Renormalization Group) convenzionali possono ottenere un numero limitato di stati eccitati, ma è difficile ottenere una relazione di dispersione completa. Gli algoritmi dinamici (come il DMRG dinamico o l'evoluzione temporale) sono alternativi valide ma computazionalmente costosi.
• L'Ansatz di Eccitazione (EA): L'approccio basato su MPS, noto come Excitation Ansatz (EA), rappresenta un'eccitazione come una sovrapposizione di onde piane di stati MPS, dove un tensore su un sito è diverso dallo stato fondamentale. Sebbene molto efficiente, l'EA standard richiede una diagonalizzazione completa separata per ogni momento $k$, rendendo il calcolo della dispersione oneroso. Inoltre, la convergenza del metodo VUMPS (Variational Uniform MPS) per lo stato fondamentale può essere instabile a seconda delle matrici iniziali.

2. Metodologia Proposta: SBEA

L'autore introduce il Site Basis Excitation Ansatz (SBEA), una variazione semplice ed efficiente dell'EA che risolve i problemi di convergenza e riduce drasticamente il costo computazionale.

A. Generazione dello Stato Fondamentale (iMPS)

Invece di utilizzare VUMPS, l'autore propone un metodo alternativo basato sul DMRG su sistemi finiti per ottenere uno stato fondamentale uniforme infinito (iMPS):

1. Si esegue un DMRG su una catena aperta lunga ($N \gg \xi$, dove $\xi$ è la lunghezza di correlazione).
2. Si inserisce un nuovo sito al centro della catena e si ottimizza il tensore centrale tramite un calcolo Lanczos (simile a un passo DMRG a sito singolo) per ottenere uno stato uniforme.
3. Si converte lo stato nella forma canonica simmetrica ($...AAA...$) utilizzando la procedura di Orus-Vidal.

• Vantaggio: Questo metodo è semplice, robusto e non richiede ricorsioni complesse, funzionando bene anche per catene a spin intero.

B. Costruzione della Base di Eccitazione

Invece di ottimizzare un tensore $B(k)$ per ogni momento $k$ separatamente, l'SBEA costruisce una base minima di tensori $\{B_\alpha\}$ che può rappresentare tutte le eccitazioni $B(k)$ simultaneamente:

1. Calcolo Lanczos a sito singolo: Si esegue un calcolo Lanczos su un singolo sito (inserito in un ambiente di stato fondamentale) per ottenere i primi $N_\alpha$ autostati dell'Hamiltoniana efficace.
2. Troncamento Pre-diagonalizzazione: Per limitare la dimensione della base, si tronca l'ambiente del tensore $B$ utilizzando i valori singolari $\lambda$ (riducendo la dimensione del legame a $\chi_{tr}$). Questo permette di ottenere una base compatta (es. $N_\alpha = 7$) che copre l'intera banda di dispersione.
3. Assenza di Gauge Sinistro: Un punto cruciale è non imporre la condizione di gauge sinistro (left-orthonormal gauge). Imporre tale gauge renderebbe gli stati locali ortogonali ma li sposterebbe a energie più alte, rendendo la base inefficiente. Mantenere la base non ortogonale è essenziale per la convergenza rapida con pochi stati.

C. Teoria delle Bande Non Ortogonale

Una volta ottenuti i tensori $B_\alpha$, si calcolano:

• Il nucleo di sovrapposizione: $O_{\alpha\alpha'}(j'-j) = \langle B_\alpha | B_{\alpha'} \rangle_{j'}$
• Il nucleo dell'Hamiltoniana: $H_{\alpha\alpha'}(j'-j) = \langle B_\alpha | H | B_{\alpha'} \rangle_{j'}$

Grazie alla lunghezza di correlazione finita, questi nuclei decadono esponenzialmente e possono essere troncati. Per ogni momento $k$, si risolve un problema agli autovalori generalizzato di dimensione ridotta ($N_\alpha \times N_\alpha$):
$$\tilde{H}(k)c = E(k)\tilde{O}(k)c$$
Questo è analogo alla teoria delle bande per un modello tight-binding con orbitali non ortogonali. La diagonalizzazione per tutti i $k$ diventa un'operazione banale e istantanea.

D. Funzioni di Wannier

L'autore dimostra come costruire eccitazioni di Wannier (localizzate nello spazio reale) combinando linearmente gli stati di momento. Queste funzioni sono ortonormali e localizzate, offrendo un punto di partenza per Hamiltoniani effettivi rinormalizzati per eccitazioni multi-magnone.

3. Risultati Chiave

Il metodo è stato testato sulla catena di Heisenberg antiferromagnetica a spin $S=1$ (catena di Haldane).

• Accuratezza: Il metodo SBEA produce la dispersione del magnone singolo con alta accuratezza. Con una base di soli $N_\alpha = 7$ stati (ottenuti con un troncamento $\chi_{tr}=8$), si ottiene un gap di Haldane $\Delta/J = 0.4107$, in ottimo accordo con il risultato esatto $0.410479...$.
• Efficienza:
  • Il calcolo dello spettro per 500 valori di $k$ richiede pochi secondi una volta calcolati i nuclei $O$ e $H$.
  • Il tempo di calcolo è dominato dalla costruzione dei nuclei (simile a un calcolo DMRG su una catena finita di lunghezza $L_c$), ma la fase di diagonalizzazione per ogni $k$ è trascurabile.
• Ruolo del Gauge: L'articolo dimostra sperimentalmente che imporre la condizione di gauge sinistro (come fatto in lavori precedenti) degrada severamente la convergenza, rendendo necessari molti più stati per ottenere la stessa accuratezza. Lasciare la base non ortogonale è fondamentale per l'efficienza.
• Limiti: Il metodo funziona bene per la banda a un magnone. Per $k$ piccoli, dove la linea a un magnone entra nel continuo a due magnoni, l'ansatz a singola particella non è più adeguato (come previsto dalla teoria).

4. Significato e Contributi

1. Semplificazione e Robustezza: SBEA riduce la complessità concettuale e computazionale rispetto all'EA standard e al VUMPS, rendendo il calcolo delle dispersioni accessibile e stabile.
2. Efficienza Computazionale: Trasforma un problema che richiedeva una diagonalizzazione completa per ogni $k$ in un problema di diagonalizzazione di matrici piccolissime, sfruttando la struttura della base locale.
3. Nuova Prospettiva sul Gauge: Dimostra che l'ortogonalità forzata (gauge sinistro) può essere controproducente per gli ansatz di eccitazione, suggerendo che la non-ortogonalità gestita correttamente è un vantaggio.
4. Versatilità: L'approccio è facilmente generalizzabile a celle unitarie più grandi e potrebbe essere applicato a sistemi bidimensionali più complessi, superando i limiti attuali della capacità computazionale.
5. Strumenti per la Fisica Teorica: La costruzione di funzioni di Wannier localizzate apre la strada alla definizione di Hamiltoniani effettivi per studiare le interazioni tra eccitazioni (es. scattering di magnoni).

In sintesi, il lavoro di White presenta un metodo "chiavi in mano" per ottenere spettri di eccitazione ad alta precisione in sistemi 1D, combinando la semplicità del DMRG finito con l'efficienza degli ansatz MPS, superando le limitazioni tecniche dei metodi precedenti.