First-Quantized Relativistic Quantum Simulation with Periodic and Dirichlet Boundary Conditions

Sintesi Tecnica: Simulazione Quantistica Relativistica a Prima Quantizzazione con Condizioni al Contorno Periodiche e di Dirichlet

Definizione del Problema
La simulazione quantistica relativistica su domini finiti affronta una sfida fondamentale: l'operatore dell'energia cinetica relativistica è una funzione non lineare del momento, il che richiede la stima di momenti di ordine superiore del momento (specificamente $\langle \hat{p}^2 \rangle$ e $\langle \hat{p}^4 \rangle$ nel regime debolmente relativistico). Nelle codifiche su griglia a prima quantizzazione utilizzando registri qubit, l'operazione di shift naturale è ciclica. Sebbene questa traslazione ciclica si adatti perfettamente alle Condizioni al Contorno Periodiche (PBC), essa introduce accoppiamenti non fisici di "wrap-around" per le Condizioni al Contorno di Dirichlet (DBC), dove la funzione d'onda deve annullarsi ai bordi. I metodi esistenti spesso trattano le condizioni al contorno come un elemento secondario alla discretizzazione; tuttavia, questo lavoro sostiene che le condizioni al contorno siano intrinseche alla definizione dell'Hamiltoniana del dominio finito, rendendo necessari operatori di momento di reticolo e strategie di stima distinti per PBC e DBC.

Metodologia
Gli autori sviluppano un framework per simulare l'operatore a singola particella a energia positiva su un intervallo 1D di lunghezza $R$ utilizzando un registro di posizione a $L$ qubit ($N=2^L$ punti di griglia).

1. Formulazione dell'Hamiltoniana:

   • Espansione Debolmente Relativistica: L'operatore cinetico relativistico $\hat{T}_{rel}$ è espanso al quarto ordine nel momento: $\hat{T}_{rel} \approx \frac{\hat{p}^2}{2m} - \frac{\hat{p}^4}{8m^3c^2}$.
   • Implementazione PBC: L'operatore di momento di reticolo $\hat{P}_P$ è definito tramite un operatore di traslazione ciclica unitaria $\hat{A}$. I momenti $\langle \hat{P}_P^2 \rangle$ e $\langle \hat{P}_P^4 \rangle$ sono espressi interamente in termini delle parti reali dei momenti di traslazione $\langle \hat{A} \rangle$ e $\langle \hat{A}^2 \rangle$.
   • Implementazione DBC: L'obiettivo è un operatore di differenza finita a catena aperta. Gli autori costruiscono questo operatore prendendo l'operatore ciclico (valutato con spaziatura di griglia DBC) e aggiungendo termini di correzione locali ai bordi per cancellare il collegamento non fisico di wrap-around. L'operatore $\hat{P}_D$ è decomposto in una parte ciclica più termini che coinvolgono i proiettori di bordo ($\hat{E}_0, \hat{E}_1, \hat{E}_2$) che dipendono solo dalle probabilità agli estremi e dalle coerenze vicino ai bordi.

2. Stimatori Quantistici:

   • Momenti di Traslazione: Per le PBC e per il bulk delle DBC, le parti reali di $\langle \hat{A}^l \rangle$ ($l=1, 2$) sono stimate utilizzando un circuito di tipo Hadamard-test con un primitivo di controllo-unitario.
   • Sovrapposizioni di Bordo: Per le DBC, i termini di correzione sono stimati utilizzando un "identità di sovrapposizione di bordo". Questo metodo ricostruisce le coerenze $\langle |f\rangle\langle g| + |g\rangle\langle f| \rangle$ utilizzando solo le probabilità della base computazionale ($P_f, P_g$) e la probabilità dello stato di sovrapposizione $|s^+_{fg}\rangle = (|f\rangle + |g\rangle)/\sqrt{2}$. Ciò evita la necessità di una stima di fase complessa o di una tomografia completa dello stato.
   • Energia Potenziale: I potenziali diagonali sono stimati tramite campionamento diretto nella base di posizione e post-elaborazione classica.

Contributi Chiave

1. Hamiltoniane Risolte ai Bordi: La formulazione di hamiltoniane di reticolo debolmente relativistiche che distinguono esplicitamente tra PBC e DBC, garantendo che gli operatori di momento siano coerenti con i vincoli fisici del dominio.
2. Decomposizione degli Stimatori: La derivazione di stimatori osservabili che ricostruiscono i momenti richiesti $\langle \hat{P}^2 \rangle$ e $\langle \hat{P}^4 \rangle$. Fondamentalmente, lo stimatore DBC è mostrato essere una decomposizione dello stimatore ciclico più un insieme finito di termini locali al bordo, piuttosto che un primitivo di misura fondamentalmente diverso.
3. Framework di Validazione: Una strategia di validazione completa che separa tre fonti distinte di errore: discretizzazione della griglia finita, troncamento debolmente relativistico ed errore statistico di misura.

Risultati
Il framework è stato validato analiticamente e numericamente attraverso diversi benchmark:

• Particella Libera PBC: Ha dimostrato che la relazione di dispersione di reticolo converge all'energia relativistica del continuo all'aumentare della dimensione della griglia, con lo stimatore perturbativo che cattura correttamente l'errore di troncamento rispetto all'energia di radice quadrata esatta del reticolo.
• Pozzo Quadrato Infinito DBC: Ha verificato che lo stimatore ciclico corretto per il bordo riproduce gli autovalori esatti dell'operatore di differenza finita a catena aperta (e lo spettro sinusoidale analitico) con un accordo di precisione macchina (discrepanze relative $< 10^{-11}$).
• Potenziali Lisci: Testati su potenziali coseno periodici (PBC) e armonici (DBC). Il workflow completo dell'energia (cinetica + potenziale) ha mostrato un'indistinguibilità visiva tra la valutazione diretta della matrice e la ricostruzione dello stimatore, con differenze assolute massime dell'ordine di $10^{-7}$.
• Scaling Statistico: Le simulazioni a numero finito di shot hanno confermato che l'errore quadratico medio (RMSE) degli stimatori dell'energia cinetica scala come $O(M^{-1/2})$ con il numero di shot $M$, coerentemente con le normali aspettative statistiche. Lo stimatore DBC ha mostrato un prefattore leggermente maggiore a causa delle misurazioni di bordo aggiuntive, ma ha mantenuto lo stesso scaling asintotico.

Significato e Rivendicazioni
Il lavoro rivendica di aver stabilito una metodologia compatta e sistematica per la simulazione quantistica relativistica che rispetta le condizioni al contorno come componente fondamentale della definizione dell'Hamiltoniana. Gli autori sottolineano che il loro approccio non richiede nuovi primitivi di misura qualitativamente diversi per la principale correzione relativistica; esso estende semplicemente l'insieme degli osservabili misurati da $\langle \hat{A} \rangle$ a $\langle \hat{A} \rangle, \langle \hat{A}^2 \rangle$ e, per le DBC, un piccolo insieme di termini di sovrapposizione di bordo.

Il lavoro è esplicitamente inquadrato nel regime debolmente relativistico, dove il supporto del momento dello stato è piccolo rispetto a $mc$. Gli autori osservano che la validità della simulazione dipende dal fatto che lo stato fisico soddisfi questa condizione, non dalle identità degli stimatori stessi, che rimangono esatte per gli operatori di reticolo definiti. Il framework è presentato come un passo pratico verso simulazioni relativistiche affidabili di sistemi continui, con il potenziale di estensione a potenziali e geometrie più complessi.