Digital Quantum Simulation of Nonequilibrium Dynamics in the Schwinger Model under a Strong External Electric Field

Sintesi Tecnica: Simulazione Quantistica Digitale della Dinamica di Non-Equilibrio nel Modello di Schwinger sotto un Forte Campo Elettrico Esterno

Problematica
Il documento affronta la sfida di simulare la dinamica di non-equilibrio nella QED (elettrodinamica quantistica) a campo forte, specificamente l'effetto Schwinger, in cui un forte campo elettrico induce il decadimento del vuoto e la produzione di coppie particella-antiparticella. I metodi perturbativi tradizionali non riescono a descrivere in modo affidabile questo fenomeno non-perturbativo e in tempo reale. Inoltre, i metodi standard di Monte Carlo basati su integrali di cammino euclideo incontrano il "problema del segno" quando si trattano l'evoluzione in tempo reale o pesi complessi, portando a una rapida crescita degli errori statistici. Sebbene il formalismo hamiltoniano eviti il problema del segno, esso richiede la gestione di uno spazio di Hilbert esponenzialmente grande. Gli autori propongono l'uso della simulazione quantistica digitale per superare tali limitazioni, rappresentando direttamente gli stati quantistici a molti corpi e l'evoluzione unitaria temporale su qubit. L'attenzione specifica è rivolta al modello di Schwinger in (1+1) dimensioni sotto un campo elettrico esterno costante, investigando la dinamica di non-equilibrio a dimensione finita a seguito di un quench di un campo elettrico.

Metodologia
Lo studio impiega un protocollo di simulazione quantistica digitale basato sul modello di Schwinger su reticolo in (1+1) dimensioni con condizioni al contorno aperte. La metodologia consiste in tre fasi principali:

1. Formulazione del Modello e Mappatura dei Qubit:

   • La teoria nel continuo viene discretizzata utilizzando il formalismo degli fermioni scalati di Kogut-Susskind.
   • La legge di Gauss viene utilizzata per eliminare i gradi di libertà di gauge, esprimendo il campo elettrico su ogni legame come funzione della distribuzione di carica della materia e del campo al contorno. Ciò riduce il modello a un Hamiltoniano puramente fermionico con interazioni a lungo raggio.
   • Gli operatori fermionici sono mappati in operatori di Pauli tramite la trasformazione di Jordan-Wigner, risultando in un Hamiltoniano di qubit composto da termini di energia cinetica, massa e campo elettrico.

2. Preparazione dello Stato (VQE):

   • Lo stato iniziale per la dinamica è lo stato di vuoto a campo nullo ($\epsilon=0$).
   • Questo stato è preparato utilizzando il Variational Quantum Eigensolver (VQE) con un ansatz RealAmplitudes hardware-efficient e l'ottimizzatore classico SLSQP.
   • Per studiare le proprietà statiche del sistema sotto campi esterni variabili, gli autori utilizzano inoltre il VQE per trovare lo stato fondamentale e il Variational Quantum Deflation (VQD) per stimare l'energia del primo stato eccitato.

3. Evoluzione in Tempo Reale (Trotter-Suzuki):

   • Il sistema è sottoposto a un improvviso "quench" in cui un campo elettrico esterno viene attivato al tempo $t=0$.
   • L'evoluzione in tempo reale è simulata utilizzando una decomposizione di Trotter-Suzuki del secondo ordine per approssimare l'operatore di evoluzione temporale $e^{-iHt}$.
   • La simulazione è eseguita su computer classici utilizzando il framework Qiskit, senza includere rumore quantistico, per confrontare l'algoritmo con la Diagonalizzazione Esatta (ED).

Contributi Chiave e Risultati
Il documento presenta i risultati numerici per dimensioni di reticolo $N=8$ e $N=12$ (con unità di reticolo $a=m=g=1$), confrontando il protocollo di simulazione quantistica con i benchmark di Diagonalizzazione Esatta (ED).

• Proprietà Statiche dello Stato Fondamentale:

  • Lo stato fondamentale a campo nullo preparato tramite VQE raggiunge un'alta fedeltà ($>99,99\%$) rispetto ai risultati ED.
  • Le scansioni sull'intensità del campo esterno $\epsilon$ rivelano un comportamento a gradini nella condensazione chirale e un gap energetico che svanisce tra lo stato fondamentale e il primo stato eccitato a specifiche intensità di campo critiche. Ciò indica una ricostruzione del vuoto indotta dal campo, corrispondente a una transizione di fase del primo ordine in volumi finiti.
  • L'intensità del campo critico $\epsilon_c$ mostra una chiara dipendenza dalla dimensione finita da $1/N$. L'estrapolazione lineare produce un campo critico a volume infinito $\epsilon_c(\infty) \approx 0,469$, che gli autori attribuiscono all'influenza delle condizioni al contorno aperte che riducono il gap energetico.

• Dinamica di Non-Equilibrio:

  • Leggi di Conservazione: L'evoluzione di secondo ordine di Trotter preserva con successo la carica totale $U(1)$ e l'energia totale (entro la precisione numerica), validando la stabilità dell'algoritmo.
  • Separazione di Carica: All'attivazione del campo, la distribuzione spaziale della carica mostra la produzione e l'annichilazione periodica di coppie carica-anticarica vicino ai bordi aperti. L'ampiezza e il periodo di queste oscillazioni aumentano con l'intensità del campo esterno.
  • Ridistribuzione dell'Energia: L'energia del campo elettrico totale $H_E(t)$ mostra un comportamento quasiperiodico, riflettendo lo scambio di energia tra il termine del campo elettrico e il settore fermionico. Questo è distinto dalla periodicità più semplice della distribuzione di carica, suggerendo contributi complessi da parte di modi a energia più elevata.
  • Fedeltà del Vuoto: La fedeltà dello stato di vuoto $P_{vac}(t)$ decade man mano che il sistema viene allontanato dal vuoto iniziale. Nella finestra temporale iniziale ($0 \le t/a \le 1,0$), questo decadimento segue una legge esponenziale efficace, $P_{vac}(t) \sim e^{-\gamma_{eff} t}$. Il tasso di decadimento efficace $\gamma_{eff}$ aumenta monotonicamente con l'intensità del campo esterno.
  • Effetti di Dimensione Finita: Il confronto tra i sistemi $N=8$ e $N=12$ conferma che le caratteristiche dinamiche osservate (separazione di carica ai bordi, decadimento della fedeltà del vuoto) non sono artefatti numerici di piccole dimensioni di sistema ma persistono in reticoli più grandi, con sistemi più grandi che mostrano una deviazione più rapida dal vuoto iniziale.

Significatività e Rivendicazioni
Gli autori affermano che il loro lavoro dimostra la fattibilità della combinazione della preparazione dello stato basata su VQE con l'evoluzione digitale in tempo reale (tramite decomposizione di Trotter-Suzuki) per studiare la dinamica di non-equilibrio in teorie di gauge su reticolo a campo forte.

• Validazione: Il protocollo riproduce con successo le caratteristiche principali dei fenomeni indotti dal campo — specificamente la separazione di carica ai bordi, il decadimento della fedeltà dello stato di vuoto e la ridistribuzione quasiperiodica dell'energia — quando confrontato con la diagonalizzazione esatta.
• Utilità Metodologica: Lo studio funge da benchmark per la simulazione quantistica digitale in fisica delle alte energie, mostrando come l'evoluzione digitale esplicita di Trotter possa catturare efficacemente complesse dinamiche a campo forte in sistemi di reticolo finiti.
• Prospettive Future: Gli autori suggeriscono modestamente che questo approccio fornisce una base per future estensioni a dimensioni di sistema maggiori, ansatz variazionali migliorati, implementazione su hardware quantistico rumoroso e lo studio di campi esterni dipendenti dal tempo; non rivendicano un'immediata realizzazione sperimentale su l'attuale hardware rumoroso, ma enfatizzano la verifica algoritmica sul simulatore classico come un passo necessario.