Infinite matrix product states for $(1+1)$-dimensional gauge theories

Gli autori presentano una costruzione di operatori prodotto di matrice (LEMPO) che permette di rappresentare in modo locale e traslazionalmente invariante le Hamiltoniane di teorie di gauge su reticolo infinito, applicando con successo il metodo allo studio del modello di Schwinger e della QCD$_2$ adiabatica.

L'essenziale

Immagina di dover risolvere un enigma cosmico: come funzionano le forze che tengono insieme l'universo, come quella che lega i quark dentro un protone? I fisici chiamano queste forze "teorie di gauge". Il problema è che sono incredibilmente complicate, come un labirinto infinito dove ogni volta che provi a calcolare qualcosa, i numeri esplodono e il computer si blocca.

Questo articolo, scritto da un gruppo di ricercatori della Princeton University, presenta una nuova chiave magica per aprire questo labirinto. Si chiama LEMPO (Link-Enhanced Matrix Product Operator), ma per semplicità, pensiamoci come a un "Super-Teletrasporto".

Ecco la spiegazione semplice, passo dopo passo:

1. Il Problema: Il Labirinto Infinito

Immagina di avere una fila infinita di scatole (i "nodi" del reticolo). In ogni scatola c'è una particella di materia (come un elettrone) e tra una scatola e l'altra c'è un "cavo" (il campo di gauge) che le collega.
Per capire come si comportano queste particelle, devi calcolare l'energia di tutto il sistema.

• Il vecchio modo: I fisici provavano a "tagliare" i cavi e a risolvere le equazioni per le scatole, ma questo rendeva il problema disordinato e non simmetrico. Era come cercare di descrivere una catena di montaggio smontando i robot: si perde la logica del flusso.
• Il limite: Se la fila è infinita (come nell'universo reale), i computer classici non possono gestire la quantità di informazioni. È come se dovessi contare ogni granello di sabbia di tutte le spiagge del mondo contemporaneamente.

2. La Soluzione: Le "Scatole Magiche" (MPS)

I ricercatori usano una tecnica chiamata MPS (Matrix Product State).
Immagina che invece di vedere l'intero universo come un blocco unico, lo guardi come una catena di anelli collegati. Ogni anello è una "scatola" che contiene un po' di informazione sulla particella e un po' di informazione sul cavo che la collega alla scatola successiva.
Queste scatole sono collegate da "fili invisibili" (spazi virtuali). Se i fili sono corti, il sistema è semplice; se sono lunghi, il sistema è complesso.

3. L'Innovazione: Il "Teletrasporto" sui Fili (LEMPO)

Qui arriva la genialità del nuovo metodo.
Nella fisica tradizionale, per calcolare l'energia del "cavo" tra due scatole, dovevi guardare tutte le scatole precedenti e sommare i loro effetti. Era come dover contare tutti i passi fatti da un amico prima di arrivare a te per sapere quanto è stanco.

I nuovi autori hanno inventato i LEMPO.

• L'analogia: Immagina che ogni scatola abbia un "pulsante magico" nascosto sul retro (lo spazio virtuale). Invece di dover guardare tutte le scatole precedenti per capire lo stato del cavo, il nuovo metodo ti permette di premere direttamente il pulsante sul retro della scatola corrente.
• Il risultato: Questo pulsante "sa" istantaneamente tutto lo stato del cavo, anche se non lo stai guardando direttamente. È come se avessi un teletrasporto che ti permette di toccare il cavo senza dover camminare lungo tutta la catena.

4. Perché è rivoluzionario?

Prima, per studiare queste teorie su computer, si doveva usare una "finestra" limitata (un numero finito di scatole) e sperare che i bordi non disturbassero il risultato.
Con i LEMPO, i ricercatori possono studiare una catena infinita direttamente.

• Analogia: Prima dovevi guardare un piccolo pezzo di un muro infinito e indovinare come era il resto. Ora puoi vedere il muro infinito intero, perfettamente liscio, senza bordi che disturbano.

5. Cosa hanno scoperto?

Hanno usato questo nuovo "Super-Teletrasporto" per studiare due modelli specifici:

1. Il Modello di Schwinger: Una versione semplificata dell'elettromagnetismo. Hanno calcolato con precisione estrema le masse delle particelle e come si comportano, ottenendo risultati che coincidono perfettamente con la teoria, ma con una precisione mai vista prima.
2. Adjoint QCD2: Un modello più complesso che assomiglia alla forza nucleare forte (quella che tiene insieme i nuclei atomici). Hanno scoperto dettagli sulle "tensioni" delle stringhe di energia e su come le particelle si comportano quando la massa cambia, confermando teorie sulla "supersimmetria" (un concetto dove ogni particella ha un "gemello" speculare).

In sintesi

Immagina di dover risolvere un puzzle di un milione di pezzi.

• Metodo vecchio: Prendi un pezzo alla volta, provi a indovinare come si collega agli altri, ma ti perdi spesso e il puzzle non è mai completo.
• Metodo nuovo (LEMPO): Hai una lente magica che ti permette di vedere come ogni pezzo si collega istantaneamente a tutto il resto, anche se il puzzle è infinito.

Questo metodo permette ai fisici di esplorare l'universo con una chiarezza senza precedenti, usando computer classici per risolvere problemi che sembravano impossibili, avvicinandoci a capire i segreti più profondi della materia.

Sommario tecnico

1. Il Problema

Lo studio delle teorie di gauge non perturbative, come il confinamento dei colori e la generazione del gap di massa, rappresenta una sfida fondamentale nella fisica moderna. Le teorie di gauge in $(1+1)$ dimensioni (spazio-tempo) fungono da modelli giocattolo gestibili per indagare queste dinamiche fortemente accoppiate.
Sebbene l'approccio Hamiltoniano ai reticoli (introdotto da Kogut e Susskind) offra vantaggi significativi rispetto all'approccio Euclideo di Wilson (assenza del problema del segno, accesso diretto allo spettro, dinamica in tempo reale), la dimensione esponenziale dello spazio di Hilbert rende la diagonalizzazione esatta impraticabile per reticoli di dimensioni moderate o superiori.
Le tecniche di rete tensoriale, in particolare gli stati prodotto matriciale (MPS), hanno avuto successo nella fisica della materia condensata, ma la loro applicazione alle teorie di gauge ha incontrato ostacoli:

• Rimozione dei gradi di libertà di gauge: Le strategie precedenti eliminavano i campi di gauge tramite la legge di Gauss, rendendo le interazioni non locali e rompendo l'invarianza traslazionale manifesta, il che impedisce l'uso di tecniche su reticoli infiniti.
• Troncamento manuale: L'aggiunta di siti fisici per rappresentare i campi di gauge richiede un troncamento manuale dello spazio di Hilbert infinito, introducendo complessità e limitazioni, specialmente per teorie non-abeliane.

2. Metodologia

Gli autori propongono una nuova costruzione basata su Stati Prodotto Matriciale Simmetrici e introducono un nuovo tipo di operatore: gli Operatori Prodotto Matriciale Potenziati dai Link (LEMPO - Link-Enhanced Matrix Product Operators).

• MPS Simmetrici: Utilizzano MPS che rispettano simmetrie globali. In questo contesto, i vincoli di simmetria locale sui tensori MPS sono identificati con la legge di Gauss della teoria di gauge. Di conseguenza, gli spazi virtuali (i legami tra i tensori) codificano naturalmente i gradi di libertà del campo di gauge, senza bisogno di gauge-fissare esplicitamente il sistema o rimuovere i gradi di libertà di gauge.
• LEMPO (Link-Enhanced MPO): A differenza dei tradizionali MPO che agiscono solo sugli spazi fisici, i LEMPO sono costruiti per agire sia sugli spazi fisici (materia) che sugli spazi virtuali (campi di gauge).
  • Gli operatori di campo elettrico (es. $L_n$) agiscono direttamente sugli indici virtuali dei tensori MPS.
  • Questo permette di esprimere l'Hamiltoniana della teoria di gauge in una forma locale e manifestamente invariante per traslazione, anche su reticoli infiniti.
• Implementazione Numerica: Il metodo è implementato utilizzando l'algoritmo VUMPS (Variational Uniform Matrix Product States) per trovare lo stato fondamentale su reticoli infiniti e l'ansatz di quasi-particelle per calcolare lo spettro delle eccitazioni. Il codice è stato sviluppato come estensione del pacchetto Julia MPSKit.jl.

3. Contributi Chiave

1. Costruzione LEMPO: Introduzione di una classe di operatori che agiscono sugli spazi virtuali, permettendo di trattare i campi di gauge in modo locale e simmetrico.
2. Generalità: Il metodo funziona sia per teorie di gauge abeliane (U(1)) che non-abeliane (SU(N)), superando le limitazioni delle strategie precedenti che faticavano a generalizzarsi al caso non-abeliano.
3. Assenza di Troncamento Esplicito: Non è necessario imporre un cutoff manuale sullo spazio di Hilbert del campo di gauge; la dimensione degli spazi virtuali viene ottimizzata dinamicamente dagli algoritmi variazionali.
4. Applicabilità a Retici Infiniti: Grazie alla preservazione dell'invarianza traslazionale, il metodo permette studi precisi nel limite termodinamico, evitando effetti di bordo indesiderati.

4. Risultati Numerici

Gli autori applicano il metodo a due modelli specifici:

A. Modello di Schwinger (QED in 1+1 dimensioni)

• Modello: Teoria di gauge U(1) accoppiata a fermioni di Dirac (massa $m$ e angolo $\theta$).
• Risultati:
  • Ottenuto un accordo eccellente con i risultati esatti nel caso di massa nulla ($m=0$).
  • Calcolo preciso dello spettro di massa e dei condensati chirali per $m \neq 0$ e vari $\theta$.
  • Confronto diretto con l'espansione ad accoppiamento forte su reticolo infinito: i risultati MPS mostrano una convergenza molto più rapida verso il limite continuo rispetto alle approssimazioni di Padé tradizionali.
  • Mappatura dettagliata dello spettro in funzione di $m$ e $\theta$, inclusa la transizione di fase a $\theta = \pi$.

B. Adjoint QCD$_2$ (SU(2) e SU(3))

• Modello: Teoria di gauge non-abeliana con fermioni di Majorana nella rappresentazione aggiunta.
• Risultati:
  • Calcolo ad alta precisione delle masse degli stati legati (bosoni e fermioni) e del condensato chirale.
  • Simmetria Supersimmetrica: Conferma numerica della supersimmetria $(1,1)$ al valore critico della massa $m_{SUSY} = g\sqrt{N_c/2\pi}$, evidenziando la degenerazione bosone-fermione nel settore di flusso banale e la presenza di un Goldstino senza massa nei settori non banali.
  • Tensione della Stringa: Calcolo della tensione della stringa fondamentale $\sigma$. Si osserva che $\sigma = 0$ per $m=0$ (comportamento di legge perimetrale) e cresce con $m$, avvicinandosi al valore della teoria di Yang-Mills pura per grandi masse.
  • Diagrammi di Fase: Mappatura delle fasi topologiche e dei punti critici per $N_c=2$ e $N_c=3$, inclusi i settori di flusso non banali.

5. Significato e Prospettive Future

Questo lavoro rappresenta un passo avanti significativo nella simulazione numerica delle teorie di gauge:

• Precisione: Fornisce stime di precisione per quantità fisiche che migliorano i risultati precedenti ottenuti con diagonalizzazione esatta o metodi di Monte Carlo.
• Scalabilità: La capacità di gestire teorie non-abeliane su reticoli infiniti apre la strada allo studio di fenomeni complessi come la dinamica in tempo reale, le funzioni di distribuzione dei partoni e le transizioni di fase in regimi non accessibili ad altri metodi.
• Estendibilità: Gli autori discutono come questo approccio possa essere generalizzato alle teorie di gauge in $(2+1)$ dimensioni utilizzando reti tensoriali PEPS (Projected Entangled Pair States), sebbene ciò presenti sfide computazionali aggiuntive legate alla contrazione dei tensori e alla gestione dei gradi di libertà di gauge dinamici.

In sintesi, gli autori dimostrano che l'uso combinato di MPS simmetrici e LEMPO offre un quadro potente, locale e invariante per traslazione per lo studio delle teorie di gauge in 1+1 dimensioni, superando le limitazioni storiche delle tecniche di rete tensoriale in questo dominio.