Quantum Simulation of Generalized Parton Distributions in the Schwinger Model

Sintesi Tecnica: Simulazione Quantistica delle Distribuzioni di Partoni Generalizzate nel Modello di Schwinger

Enunciato del Problema
Le Distribuzioni di Partoni Generalizzate (GPD) sono fondamentali per comprendere la struttura interna degli adroni, codificando le correlazioni tra i parametri di impatto trasversali e il momento longitudinale. Tuttavia, il calcolo delle GPD a partire dai primi principi nella Cromodinamica Quantistica (QCD) è ostacolato dal problema del segno nella QCD su reticolo in spazio euclideo, che impedisce il calcolo diretto delle funzioni di correlazione sulla luce in tempo reale. Mentre gli approcci classici come la Large Momentum Effective Theory (LMET) offrono soluzioni alternative tramite le quasi-PDF, essi richiedono un estremo incremento del momento e soffrono di correzioni di ordine superiore. Il calcolo quantistico offre un cambio di paradigma consentendo la simulazione diretta della dinamica hamiltoniana e dell'evoluzione in tempo reale, potenzialmente superando il problema del segno. Nonostante questo potenziale, la simulazione delle GPD rimane una sfida formidabile, particolarmente per quanto riguarda la preservazione delle simmetrie discrete (Carica $\mathcal{C}$, Parità $\mathcal{P}$, Tempo $\mathcal{T}$) e la preparazione efficiente di stati hadronici con momento non nullo.

Metodologia
Gli autori propongono un algoritmo quantistico completo per simulare le GPD nel modello di Schwinger (elettrodinamica) in (1+1) dimensioni, fungendo da proxy per la piena QCD. La metodologia si distingue per diverse scelte tecniche chiave:

1. Fermioni di Wilson: A differenza delle attuali simulazioni quantistiche che utilizzano prevalentemente fermioni staggered, questo lavoro impiega i fermioni di Wilson per la discretizzazione su reticolo. Questa scelta è critica per la rigorosa preservazione della simmetria esatta di carica ($\mathcal{C}$) sul reticolo. Gli autori sostengono che la simmetria $\mathcal{C}$ sia vitale per studiare i mesoni carichi neutri, poiché elimina contaminazioni non fisiche della parte reale nelle funzioni di correlazione e assicura che le GPD rimangano funzioni dispari rispetto alla frazione di momento $x$.
2. Formulazione Hamiltoniana: Il modello di Schwinger è formulato nel gauge assiale. L'Hamiltoniana del continuo è discretizzata utilizzando i fermioni di Wilson, aumentata da un termine di Wilson per aggirare il problema del raddoppio dei fermioni pur preservando le simmetrie $\mathcal{C}$, $\mathcal{P}$ e $\mathcal{T}$. Il campo di gauge temporale viene eliminato tramite la legge di Gauss, risultando in un termine di interazione non locale nell'Hamiltoniana.
3. Preparazione dello Stato: L'algoritmo utilizza un Quantum-Number-Resolving Variational Quantum Eigensolver (VQE). Per preparare stati hadronici con specifici momenti non nulli ($|h(k)\rangle$) e il vuoto ($|\Omega\rangle$), gli autori costruiscono un insieme di stati di riferimento ortogonali. Un circuito di iniezione di momento, $D(k)$, mappa uno stato di riferimento in uno stato con momento definito. L'ansatz VQE, $U(\theta)$, è progettato per essere preservante della simmetria, ottimizzando una funzione di perdita multi-stato per estrarre simultaneamente il vuoto e gli stati eccitati hadronici all'interno del settore privo di carica ($Q=0$).
4. Misurazione della Funzione di Correlazione: Il cuore dell'algoritmo valuta la funzione di correlazione sulla luce invariante per gauge, che include le linee di Wilson. Gli autori implementano una versione controllata dell'operatore di traslazione del momento e i link di gauge controllati ($W_0$) utilizzando un qubit ausiliario. Misurando i valori di aspettazione off-diagonal ($\sigma_x$ e $\sigma_y$) del qubit ausiliario dopo l'evoluzione in tempo reale, si estraggono le funzioni di correlazione spaziale $\tilde{H}_{q/h}(N, \xi)$.
5. Analisi della Complessità: I requisiti di risorse sono analizzati utilizzando la formula del prodotto di Suzuki-Trotter del secondo ordine per l'evoluzione temporale. La complessità totale scala polinomialmente con la dimensione del sistema $N$ e l'inversa precisione $\epsilon$, specificamente come $O(N^p + N^{5/2}/\epsilon^{1/2})$.

Contributi Chiave

• Framework Unificato: Il documento stabilisce un framework algoritmico quantistico unificato capace di calcolare sia le Funzioni di Distribuzione dei Partoni (PDF) che le GPD, anziché trattarle come valutazioni disgiunte.
• Preservazione della Simmetria: Dimostra che l'implementazione dei fermioni di Wilson in un framework di calcolo quantistico non comporta una maggiore complessità computazionale rispetto ai fermioni staggered, ma offre proprietà di simmetria superiori. Nello specifico, preserva rigorosamente la simmetria $\mathcal{C}$ esatta, che si è dimostrato essenziale per ottenere risultati fisicamente significativi per i mesoni neutri.
• Costruzione Algoritmica: Il lavoro dettaglia la costruzione di circuiti quantistici specifici per l'iniezione di momento ($D(k)$), l'implementazione del link di gauge e l'ansatz VQE preservante della simmetria richiesto per il modello di Schwinger.

Risultati
Gli autori hanno testato il loro approccio utilizzando la Diagonalizzazione Esatta (ED) su un reticolo spaziale di dimensione $N=18$ (e fino a $N=22$ per controlli specifici):

• Spettro di Massa: Lo spettro di massa calcolato dello stato legato $q\bar{q}$ più leggero ha riprodotto il noto limite analitico del continuo ($m_h/g \approx 0.564$) tramite estrapolazione lineare, validando la formulazione hamiltoniana.
• Funzioni di Correlazione: Le funzioni di correlazione spaziali sulla luce $\tilde{H}_{q/h}(N, \xi)$ sono state estratte. I risultati hanno confermato che la parte reale della funzione di correlazione svanisce asintoticamente all'aumentare del volume del reticolo, in accordo con la simmetria $\mathcal{C}$.
• Confronto con i Fermioni Staggered: Un confronto diretto ha mostrato che il rapporto tra la parte reale e quella immaginaria della funzione di correlazione è significativamente soppresso per i fermioni di Wilson rispetto ai fermioni staggered. Ciò conferma che la formulazione di Wilson elimina efficacemente le contaminazioni non fisiche derivanti da effetti di volume finito e violazioni di simmetria.
• GPD e PDF: Le trasformate di Fourier delle funzioni di correlazione hanno prodotto le GPD nello spazio dei momenti. Le PDF risultanti (limite $\xi=0$) presentano un picco a $x=0.5$, riflettendo correttamente la natura a due quark di valenza dello stato legato in 1+1 dimensioni. Le GPD sono risultate essere funzioni dispari di $x$, soddisfacendo il requisito teorico per cui il fattore di forma della carica elastica per uno stato a carica neutra deve svanire.

Significatività e Rivendicazioni
Il documento sostiene che questo lavoro rappresenti un importante precursore per le simulazioni della QCD completa in (3+1) dimensioni. Dimostrando che i fermioni di Wilson possono essere implementati efficientemente su hardware quantistico senza aumentare la complessità, gli autori argomentano che questo approccio offre una strada robusta verso la simulazione di fenomeni non perturbativi rispettando rigorosamente le simmetrie fondamentali. La scalabilità polinomiale dell'algoritmo rispetto alla dimensione del sistema e alla precisione evidenzia il vantaggio teorico dei computer quantistici per la simulazione diretta delle correlazioni sulla luce. Gli autori pongono questo lavoro come un passo fondamentale verso la generalizzazione di questi algoritmi a teorie di gauge non abeliane (ad esempio, Yang-Mills SU(2) o SU(3)) e, infine, alla piena QCD su futuri hardware quantistici tolleranti ai guasti.