Eigenstate Thermalization Hypothesis (ETH) for off-diagonal matrix elements in integrable spin chains

Lo studio utilizza l'Ansatz di Bethe Algebraico per dimostrare che, nelle catene di spin integrabili, gli elementi di matrice off-diagonali degli operatori locali decadono esponenzialmente con la dimensione del sistema sia per stati nella stessa macrostato termodinamico (con distribuzione di Gumbel) sia, in modo più rapido, per stati in macrostati diversi.

L'essenziale

Il Grande Mistero: Perché il mondo si "raffredda"?

Immagina di avere una stanza piena di persone (gli atomi o le particelle) che ballano tutte insieme seguendo una musica complessa. Se la musica è caotica e imprevedibile, dopo un po' le persone smettono di ballare in modo sincronizzato e si mescolano: la stanza raggiunge una temperatura uniforme. Questo è il thermalization (termalizzazione).

In fisica, c'è una regola d'oro chiamata ETH (Eigenstate Thermalization Hypothesis). Dice che, in sistemi "caotici" (come la maggior parte dei materiali reali), ogni singola "canzone" (stato energetico) contiene già in sé la statistica della stanza intera. Se guardi una sola nota della canzone, puoi capire come si comporta l'intera orchestra.

Ma cosa succede se la musica non è caotica, ma è perfettamente ordinata? Come una sinfonia dove ogni strumento sa esattamente cosa fare e non si mescola mai con gli altri? Questi sono i sistemi integrabili (come la catena di spin di Heisenberg studiata in questo articolo). Qui, la regola ETH classica dovrebbe rompersi.

Cosa hanno scoperto gli scienziati?

Rottoli e Alba hanno preso un sistema matematico molto specifico (una catena di spin, che puoi immaginare come una fila di calamite minuscole che possono puntare su o giù) e hanno chiesto: "Se prendiamo due note diverse di questa sinfonia ordinata e vediamo come interagiscono, succede qualcosa di strano?"

Hanno analizzato i matrici di elementi fuori diagonale. Tradotto in parole povere: hanno guardato quanto due stati energetici diversi "parlano" tra loro quando applichiamo un piccolo disturbo (un operatore locale).

Ecco i tre punti chiave della loro scoperta, spiegati con analogie:

1. La regola della "Soglia di Silenzio" (Decadimento Esponenziale)

In un sistema caotico, il "dialogo" tra due stati diversi diventa rapidamente impercettibile man mano che il sistema cresce. È come se due persone in una folla enorme iniziassero a sussurrare: più la folla è grande, meno si sentono.
Gli scienziati hanno scoperto che anche nei sistemi ordinati (integrabili), questo sussurro diventa inudibile molto velocemente.

• L'analogia: Immagina di cercare di sentire un sussurro in un stadio. Se lo stadio è piccolo, lo senti. Se lo stadio è enorme (come quelli studiati qui, con centinaia di posti), il sussurro sparisce quasi istantaneamente.
• La differenza: Nei sistemi caotici, il silenzio arriva con una certa velocità. Nei sistemi ordinati, se i due stati appartengono a "gruppi sociali" diversi (macrostati diversi), il silenzio arriva molto più velocemente (come un muro di gomma che blocca tutto). Se appartengono allo stesso gruppo, il silenzio arriva un po' più lentamente, ma comunque in modo esponenziale.

2. La forma della "Paura" (Distribuzione di Gumbel)

Qui arriva la parte più affascinante. Quando misuriamo quanto forte è questo sussurro (la probabilità che due stati interagiscano), ci aspettiamo che i risultati seguano una campana classica (Gaussiana), come quando lanci monete o misuri l'altezza delle persone.

• La sorpresa: Invece, in questi sistemi ordinati, i risultati non formano una campana. Formano una curva chiamata Distribuzione di Gumbel.
• L'analogia: Immagina di lanciare un dado. Se lanci 100 volte, ottieni una distribuzione normale. Ma se invece di lanciare un dado, guardi il record di qualcosa (es. l'onda più alta che ha colpito una spiaggia in un anno), i record estremi seguono la distribuzione di Gumbel.
  Gli scienziati hanno scoperto che i "sussurri" in questi sistemi integrabili si comportano come record estremi. Non è una media, è una statistica di eventi rari e speciali. È come se il sistema dicesse: "Di solito non parlo, ma quando parlo, lo faccio in modo molto specifico e non casuale".

3. Il problema delle "Stringhe" (Stati Legati)

In questi sistemi ordinati, le particelle possono legarsi insieme formando "pacchetti" chiamati stringhe (come un gruppo di amici che camminano sempre tenendosi per mano).
Fino a poco tempo fa, calcolare come queste "stringhe" interagivano era un incubo matematico perché le formule si rompevano (diventavano infinite).

• Il trucco: Gli autori hanno usato una tecnica avanzata (l'Ansatz di Bethe Algebrico) che permette di "aggiustare" le formule, rimuovendo i buchi matematici. È come se avessero trovato un modo per contare le persone in una stanza anche se alcune di loro sono incollate l'una all'altra, senza contare due volte lo stesso gruppo.
• Risultato: Hanno dimostrato che anche con queste "stringhe" (legami complessi), la regola del silenzio esponenziale e la forma "Gumbel" restano vere.

In sintesi: Perché è importante?

Prima di questo lavoro, pensavamo che i sistemi ordinati (integrabili) fossero un mondo a parte, dove le regole del caos non valevano affatto.
Questo studio ci dice: "No, non è così!".
Anche in un sistema perfettamente ordinato, c'è una forma di termalizzazione locale. Le particelle smettono di parlarsi tra loro in modo significativo quando il sistema diventa grande, proprio come nei sistemi caotici.
La differenza sta solo nel come si comportano statisticamente: invece di essere casuali come un lancio di moneta (Gaussiano), seguono una logica più "estrema" (Gumbel), tipica dei sistemi che hanno molte regole di conservazione.

La morale della favola: Anche in un universo perfettamente ordinato e prevedibile, le cose locali tendono a "dimenticarsi" l'una dell'altra molto velocemente, ma lo fanno con un'impronta digitale statistica unica e diversa dal caos. Gli scienziati hanno finalmente trovato il modo di leggere questa impronta digitale.

Sommario tecnico

Titolo: Ipotesi di Termalizzazione degli Autostati (ETH) per gli elementi di matrice fuori diagonale in catene di spin integrabili

1. Il Problema e il Contesto

L'articolo affronta uno dei problemi fondamentali della fisica quantistica dei molti corpi fuori equilibrio: comprendere come l'evoluzione unitaria di sistemi isolati porti alla termalizzazione locale.

• Contesto: L'ipotesi di termalizzazione degli autostati (ETH) è il paradigma standard per i sistemi caotici (non integrabili). Essa afferma che gli elementi di matrice di un operatore locale $O$ tra autostati dell'energia $|E_i\rangle$ e $|E_j\rangle$ seguono una legge specifica: gli elementi diagonali sono funzioni lisce dell'energia, mentre gli elementi fuori diagonale sono soppressi esponenzialmente con la dimensione del sistema $L$ e distribuiti secondo una variabile casuale gaussiana.
• La Sfida: Nei sistemi integrabili, la presenza di un numero esteso di cariche conservate viola l'ETH standard. Sebbene sia noto che gli elementi diagonali decadono come $L^{-1}$ (invece che esponenzialmente), la statistica e il decadimento degli elementi di matrice fuori diagonale in modelli integrabili reticolari (come le catene di spin) rimanevano poco esplorati, specialmente alla luce di recenti risultati sul modello di campo continuo Lieb-Liniger (che non presenta stati legati).
• Obiettivo: Investigare il comportamento degli elementi di matrice fuori diagonale in una catena di spin integrabile con stati legati (stringhe), specificamente la catena di Heisenberg isotropa spin-1/2 (modello XXX), utilizzando operatori con supporto su uno o due siti.

2. Metodologia

Gli autori impiegano risultati all'avanguardia dell'Ansatz di Bethe Algebrico (ABA) per calcolare efficientemente gli elementi di matrice di operatori locali tra autostati generici, superando i limiti computazionali della diagonalizzazione esatta.

• Modello: Catena di Heisenberg isotropa (XXX) con Hamiltoniana $H = \frac{1}{4} \sum \vec{\sigma}_i \cdot \vec{\sigma}_{i+1}$.
• Strumenti Computazionali:
  • Utilizzo delle formule esatte per gli elementi di matrice di operatori locali (basate sulle formule di Slavnov e Gaudin).
  • Gestione delle singolarità fittizie introdotte dal limite omogeneo e dalla presenza di stringhe di Bethe (stati legati di quasiparticelle). Per gli operatori a un sito, le singolarità possono essere rimosse analiticamente; per operatori a due siti, si restringe l'analisi agli autostati privi di stringhe.
  • Campionamento degli autostati tramite un algoritmo Metropolis per generare ensemble termici (macrostati) a diverse temperature (inclusi $\beta=0$ e $\beta=\infty$) mantenendo fissa la densità di particelle $M/L$.
• Osservabili:
  • Operatori a un sito: $E^{(11)}_i = \frac{1}{2}(1 + \sigma^z_i)$.
  • Operatori a due siti: $E^{(11)}_i E^{(22)}_{i+1}$.
  • Grandezza studiata: $M_{ij} = \ln |\langle E_i | O | E_j \rangle|^2$.

3. Risultati Chiave

A. Autostati nello stesso macrostato termico

• Decadimento: Gli elementi di matrice fuori diagonale decadono esponenzialmente con la dimensione del sistema $L$, ma con una correzione logaritmica:
  $$|\langle E_i | O | E_j \rangle| \propto \exp\left( -c_O L \ln L - c_0 L \right)$$
  Questo conferma che la presenza di stati legati (stringhe) non altera qualitativamente il quadro osservato nel modello Lieb-Liniger.
• Distribuzione Statistica: La distribuzione di probabilità (PDF) della quantità $M_{ij}$ (o meglio, della sua parte logaritmica) non è gaussiana, come previsto dall'ETH standard per i sistemi caotici.
  • La distribuzione è ben descritta da una distribuzione di Gumbel.
  • I parametri della distribuzione di Gumbel sono indipendenti da $L$ nel limite termodinamico.
  • Questo risultato è in netto contrasto con il modello Lieb-Liniger (dove si osservava una distribuzione di Fréchet) e con i sistemi caotici (distribuzione gaussiana).
• Non-Gaussianità: Gli autori dimostrano che la non-gaussianità è rilevabile tramite il rapporto $\Gamma = \overline{|O_{ij}|^2} / |\overline{O_{ij}}|^2$. Mentre per variabili gaussiane $\Gamma = \pi/2$, nel modello XXX questo rapporto diverge esponenzialmente con $L$ (raggiungendo valori dell'ordine di $10^5$), confermando la natura non gaussiana delle fluttuazioni.

B. Autostati in macrostati diversi

• Decadimento: Quando si considerano elementi di matrice tra autostati appartenenti a macrostati termodinamici diversi (es. stato fondamentale vs stato a temperatura infinita), il decadimento è molto più rapido:
  $$|\langle E_i | O | E_j \rangle| \propto \exp(-d_O L^2)$$
• Statistica: Anche in questo caso, la distribuzione dei logaritmi degli elementi di matrice segue una legge di Gumbel.

C. Operatori a due siti

• Nonostante le limitazioni computazionali (dimensioni fino a $L \approx 60$ a causa della complessità esponenziale nel supporto dell'operatore), i dati numerici per operatori a due siti confermano lo stesso scenario qualitativo osservato per gli operatori a un sito: decadimento esponenziale con correzione logaritmica e distribuzione di Gumbel.

4. Contributi e Significato

1. Estensione ai Modelli Reticolari: Il lavoro estende la comprensione della statistica degli elementi di matrice fuori diagonale dai modelli di campo continuo (Lieb-Liniger) ai modelli reticolari integrabili, dimostrando che il quadro teorico è robusto anche in presenza di stati legati complessi.
2. Ruolo degli Stati Legati: Dimostra che la presenza di stringhe di Bethe (stati legati) non cambia la scala di decadimento esponenziale, ma influenza la statistica fine (distribuzione di Gumbel invece che Fréchet).
3. Violazione dell'ETH Standard: Fornisce una prova numerica solida che, sebbene la scala di decadimento sia simile a quella ETH ($e^{-L}$), la statistica sottostante è radicalmente diversa (Gumbel vs Gaussiana). Questo suggerisce che l'integrabilità modifica la natura delle fluttuazioni quantistiche in modo fondamentale.
4. Metodologia Computazionale: Il successo nel calcolare elementi di matrice per $L \sim 500$ (per operatori a un sito) utilizzando l'ABA, superando i limiti della diagonalizzazione esatta ($L \sim 20-30$), rappresenta un avanzamento tecnico significativo nello studio delle proprietà di termalizzazione in sistemi integrabili.
5. Implicazioni Future: I risultati aprono la strada a nuove indagini sulla relazione tra l'ETH, la teoria della probabilità libera (Free Probability Theory) e la dinamica di rilassamento in sistemi integrabili, suggerendo che la statistica di Gumbel potrebbe essere un segno distintivo universale per certi tipi di sistemi integrabili.

In sintesi, il paper stabilisce che, sebbene i sistemi integrabili mostrino un decadimento esponenziale degli elementi di matrice fuori diagonale simile ai sistemi caotici, la loro statistica interna è governata da leggi di distribuzione estreme (Gumbel) piuttosto che gaussiane, riflettendo la struttura peculiare dello spazio di Hilbert degli stati legati.