Lattice Quantum Chromodynamics for Quantum Simulations

Sintesi Tecnica: Cromodinamica Quantistica su Reticolo per Simulazioni Quantistiche

Enunciato del Problema
La Cromodinamica Quantistica (QCD) è la teoria fondamentale dell'interazione forte, eppure il suo regime non perturbativo — essenziale per comprendere la struttura degli adroni e la materia ad alta densità — rimane computazionalmente impegnativo. Gli approcci standard della QCD su reticolo si affidano a integrali di cammino euclidei e campionamento per importanza, che falliscono in presenza di un "problema del segno" (sign problem), come nel caso della dinamica in tempo reale o quando viene incluso un angolo $\theta$. Sebbene le formulazioni hamiltoniane offrano una via per la simulazione in tempo reale, la crescita esponenziale dello spazio di Hilbert con il volume del reticolo rende la simulazione classica intrattabile per sistemi di grandi dimensioni. Di conseguenza, esiste una critica necessità di framework che mappino le teorie di gauge su reticolo su computer quantistici, affrontando specificamente la codifica dei gradi di libertà di gauge e fermionici, la costruzione di stati fisici (invarianti per gauge) e lo sviluppo di algoriti efficienti di evoluzione temporale.

Metodologia
Gli autori sviluppano un framework completo per simulazioni quantistiche della teoria di gauge $SU(N_c)$ con quark, concentrandosi su $N_c=3$ (QCD) in due e tre dimensioni spaziali. La metodologia si basa sui seguenti pilastri:

1. Formulazione Hamiltoniana e Discretizzazione:
   Il lavoro utilizza la formulazione hamiltoniana di Kogut-Susskind. Gli autori considerano sia i fermioni staggered che quelli di Wilson per affrontare il problema del raddoppio dei fermioni. La teoria è discretizzata su un reticolo spaziale cubico con tempo continuo. L'Hamiltoniana è decomposta in termini elettrici ($H_E$), magnetici ($H_B$), $\theta$ ($H_\theta$), cinetici ($H_{kin}$) e di massa ($H_{mass}$). Il termine $\theta$ è incluso come accoppiamento alla densità di carica topologica, introducendo coefficienti complessi nell'Hamiltoniana.

2. Base di Rappresentazione e Spazio di Hilbert Fisico:
   Un contributo centrale è la formulazione dello spazio di Hilbert fisico in una base di rappresentazione. Invece di mappare direttamente gli elementi di gruppo su qubit, i gradi di libertà di gauge sui link e i gradi di libertà fermionici sui siti sono codificati utilizzando le rappresentazioni irriducibili (irreps) di $SU(N_c)$.

   • Spazio di Hilbert del Link: I link sono descritti da stati $|(R, r_R, r_L)\rangle$, dove $R$ è una irrep e $r_{R,L}$ sono vettori di base per le semilinks destra e sinistra. Questa base diagonalizza gli operatori di Casimir quadratico.
   • Spazio di Hilbert del Sito: I fermioni sono trattati tramite una base di numero di occupazione, che viene poi trasformata in una base di rappresentazione utilizzando i coefficienti di Clebsch-Gordan (CGC). Ciò consente di classificare gli stati particella e antiparticella in base alle loro proprietà di trasformazione sotto $SU(N_c)$.
   • Invarianza di Gauge: Gli stati fisici sono costruiti come singoletti (rappresentazioni triviali) formati dal prodotto tensoriale di tutte le irreps che si incontrano in un sito del reticolo (link e fermioni). Gli autori definiscono esplicitamente la convenzione di "ordine del sito" necessaria per costruire in modo univoco questi stati singoletti, garantendo che il vincolo della legge di Gauss sia soddisfatto.

3. Troncamento e Codifica:
   Poiché lo spazio di Hilbert di $SU(N_c)$ è a dimensione infinita, gli autori implementano schemi di troncamento. Questi includono cutoff sul peso massimo delle irreps ($n_1$) o sull'autovalore del Casimir quadratico ($C_2$) per i link di gauge, e limiti sui numeri di occupazione cromatica per i fermioni. Gli stati fisici sono codificati in qubit utilizzando uno schema di codifica locale, dove i qubit rappresentano specifiche irreps e indici di molteplicità del singoletto. Questo approccio è notato come più efficiente in termini di qubit rispetto alle codifiche dirette per grandi spazi di Hilbert.

4. Algoritmi di Simulazione Quantistica:
   Gli autori propongono un algoritmo di Trotterizzazione per l'evoluzione temporale. L'operatore di evoluzione temporale $e^{-iHt}$ è approssimato come un prodotto di unitari corrispondenti ai termini decomposti dell'Hamiltoniana.

   • Costruzione delle Porte: L'evoluzione degli operatori di plaquette e di link è mappata in rotazioni di Givens (unitari a due livelli) e spostamenti di fase. Gli autori dettagliano la costruzione dei circuiti quantistici per queste operazioni, incluso il trattamento delle stringhe di Jordan-Wigner (flip di fase) necessari per le relazioni di anticommutazione fermionica.
   • Elementi di Matrice: Il documento fornisce formule esplicite per gli elementi di matrice di tutti gli operatori dell'Hamiltoniana, espressi come prodotti di "fattori di sito" che coinvolgono i CGC. Queste formule sono essenziali per calcolare gli angoli di rotazione e le fasi richiesti per le porte quantistiche.

Contributi Chiave

• Framework per $SU(3)$ con Angolo $\theta$: Il lavoro presenta la prima formulazione e simulazione della QCD su reticolo con un angolo $\theta$ in tre dimensioni spaziali utilizzando un approccio hamiltoniano su un computer quantistico.
• Formalismo della Base di Rappresentazione: Fornisce una costruzione rigorosa dello spazio di Hilbert fisico utilizzando le irreps e i singoletti CGC, offrendo un'alternativa scalabile alle basi basate sugli elementi di gruppo.
• Dettagli Algoritmici: Il lavoro dettaglia i specifici circuiti quantistici per l'evoluzione Trotterizzata, inclusa la decomposizione di termini complessi dell'Hamiltoniana (come il termine $\theta$ e i termini cinetici di Wilson) in porte elementari.
• Calcoli degli Elementi di Matrice: Gli autori derivano e tabulano gli elementi di matrice per gli operatori di plaquette, $\theta$ e cinetici fermionici, fornendo i dati necessari per implementare i circuiti di simulazione.

Risultati
Gli autori hanno eseguito simulazioni senza rumore su piccoli reticoli (fino a un reticolo periodico $2\times2$ e un cubo a bordi aperti $2\times2\times2$) utilizzando fino a 32 qubit su simulatori classici (Qiskit Aer e CUDA-Q). Le simulazioni hanno dimostrato:

• Effetti di Troncamento dello Spazio di Hilbert: Analisi di come diversi troncamenti di gauge e fermioni influenzino gli osservabili come il campo elettrico medio ($\langle \Pi^2 \rangle$) e il numero di occupazione cromatica ($\langle N_q \rangle$). I risultati hanno indicato che il forte accoppiamento sopprime gli stati ad alta popolazione, rendendo efficaci i troncamenti modesti per i fermioni staggered, mentre i fermioni di Wilson hanno mostrato una maggiore sensibilità al troncamento.
• Dipendenza dai Parametri: L'evoluzione degli osservabili è stata studiata in funzione della massa fermionica e del parametro di Wilson $r$. Masse più pesanti sono risultate capaci di sopprimere le oscillazioni verso stati ricchi di fermioni.
• Dinamica dell'Angolo $\theta$: Le simulazioni variando l'angolo $\theta$ hanno rivelato un comportamento aperiodico nel regime di forte accoppiamento, con dinamiche distinte per $\theta$ positivo e negativo.
• Stati Adronici e Dinamica delle Stringhe: Gli autori hanno dimostrato la preparazione di stati di "pioni carichi" e l'evoluzione di stati di corda (sia stringhe di loop che di Fermi). Hanno osservato la dinamica di rottura della corda, dove i tubi di flusso decadono in coppie fermione-antifermione, e hanno tracciato la probabilità di configurazioni con o senza rottura della corda.
• Potenziale Chimico dei Barioni: Utilizzando la diagonalizzazione esatta su un singolo cubo, gli autori hanno calcolato il costo energetico per l'aggiunta di un barione, servendo come prova di concetto per lo studio della materia ad alta densità.

Significatività e Rivendicazioni
Il documento si pone come un contributo fondamentale allo sviluppo di simulazioni quantistiche pratiche per la fisica fondamentale. Gli autori affermano che il loro lavoro fornisce gli "ingredienti" necessari per progettare simulazioni quantistiche della QCD su reticolo nella formulazione hamiltoniana.

• Prova di Concetto: La principale significatività risiede nel dimostrare che la complessa macchina della teoria di gauge su reticolo $SU(3)$, inclusi l'angolo $\theta$ e la materia fermionica, può essere codificata ed evoluta su hardware quantistico (simulato qui).
• Percorso di Scalabilità: La codifica locale e la base di rappresentazione sono presentate come una via percorribile verso simulazioni più grandi, richiedendo potenzialmente meno qubit rispetto alle codifiche globali.
• Prospettive Future: Gli autori dichiarano con modestia che i risultati attuali sono esplorativi e non traggono ancora conclusioni fisiche serie sulla QCD in tempo reale a causa delle piccole dimensioni del reticolo, del regime di forte accoppiamento e dei troncamenti aggressivi utilizzati. Identificano la necessità di Hamiltoniane migliorate (controtermini), migliori algoritmi di preparazione dello stato e hardware tollerante ai guasti per raggiungere il limite continuo e simulare regimi fisicamente rilevanti. Il lavoro è inquadrato come un gradino verso la simulazione della materia densa (es. stelle di neutroni) e della dinamica axion-QCD in futuro.