Relativistic Quantum Simulation under Periodic and Dirichlet Boundary Conditions: A First-Quantised Framework for Near-Term Devices

Sintesi Tecnica: Simulazione Quantistica Relativistica sotto Condizioni al Contorno Periodiche e di Dirichlet

Enunciato del Problema
Simulare sistemi quantistici è computazionalmente oneroso a causa della crescita esponenziale della complessità con la dimensione del sistema, una sfida esacerbata quando devono essere incorporati effetti relativistici per descrivere accuratamente i fenomeni fisici. Sebbene la teoria delle perturbazioni sia un approccio standard per studiare gli effetti relativistici oltre il quadro non relativistico, l'implementazione di queste simulazioni su dispositivi quantistici a breve termine rimane difficile. I metodi esistenti spesso si affidano alla seconda quantizzazione, che è uno standard per la chimica quantistica, ma potrebbe mancare di una gestione intuitiva per sistemi con non linearità. Esiste la necessità di un quadro pratico basato sulla prima quantizzazione capace di stimare le energie dello stato fondamentale relativistiche per sistemi sotto diverse condizioni al contorno su hardware near-term.

Metodologia
Gli autori propongono un quadro di prima quantizzazione per la simulazione quantistica relativistica che discretizza la funzione d'onda attraverso una griglia finita rappresentata da qubit di sistema. La metodologia centrale prevede i seguenti componenti:

1. Discretizzazione e Operatori di Traslazione:
   La funzione d'onda spaziale continua 1D viene mappata su una griglia discreta di $2^L$ punti utilizzando $L$ qubit, dove la posizione $x$ è scalata nel dominio $0 \le x < 1$. L'operatore del momento al quadrato $\hat{p}^2$ è approssimato tramite il metodo delle differenze finite basato su operatori di traslazione quantistica. L'operatore di traslazione $\hat{A}_x$ (addizionatore quantistico) sposta lo stato della posizione $|j\rangle$ in $|j+\delta x\rangle$, mentre il suo coniugato $\hat{A}_x^\dagger$ (sottrattore quantistico) sposta nella direzione opposta. Il momento al quadrato è espresso come $\hat{p}^2 = -(mc L_m)^2 (\hat{A}^{(1)} - 2\hat{1})$, dove $L_m = \lambda_m / \delta x$ è un parametro adimensionale che coinvolge la lunghezza d'onda di Compton ridotta.

2. Espansione Perturbativa dell'Energia Cinetica:
   L'Hamiltoniana dell'energia cinetica relativistica $\hat{H}_{re}^K$ è derivata dalla forma radicale della relazione energia-momento relativistica. Viene espansa in una serie perturbativa di termini di momento di ordine superiore ($\hat{p}^2, \hat{p}^4, \hat{p}^6, \dots$) valida nel regime in cui $\langle \hat{p}^2 \rangle \ll m^2 c^2$. L'energia cinetica relativistica totale è approssimata troncando questa serie, incorporando coefficienti derivati dall'equazione di Klein-Gordon.

3. Gestione delle Condizioni al Contorno:
   Il quadro affronta due distinte condizioni al contorno:

   • Condizioni al Contorno Periodiche (PBC): Gli operatori di traslazione preservano l'unitarietà avvolgendo la griglia ($\hat{A}_x |2^L-1\rangle = |0\rangle$). I valori di aspettazione delle potenze pari del momento sono espressi come combinazioni di $\langle \hat{A}^{(l)} \rangle$.
   • Condizioni al Contorno di Dirichlet (DBC): In questo caso, la funzione d'onda non si avvolge ($\psi(0) \neq \psi(1)$). Gli autori introducono un operatore di correzione $\hat{E}_0 = |2^L-1\rangle\langle 0| + |0\rangle\langle 2^L-1|$ per tenere conto del contorno. I valori di aspettazione per le DBC sono formulati come i termini PBC più termini aggiuntivi che coinvolgono $\hat{E}_0$ e operatori correlati ($\hat{E}_1, \hat{E}_2$), che dipendono dalle ampiezze di probabilità alle estremità della griglia.

4. Algoritmo di Simulazione Quantistica:
   Gli autori progettano circuiti quantistici per stimare i valori di aspettazione richiesti:

   • Energia Cinetica: Gate di traslazione controllata vengono applicati iterativamente a un qubit di controllo. Misurando il qubit di controllo nella base Z si ottiene la parte reale di $\langle (\hat{A}_x^\dagger)^l \rangle$, che corrisponde ai termini di momento necessari.
   • Energia Potenziale: Una matrice di densità $\hat{\rho}_V$ che rappresenta il potenziale viene utilizzata in un circuito per stimare $\langle \hat{V} \rangle$ tramite le probabilità dello stato di sistema.
   • Ottimizzazione Variazionale: L'energia totale è stimata minimizzando la somma dei valori di aspettazione dell'energia cinetica e potenziale su un insieme di parametri variazionali $\{\bar{\kappa}\}$ in uno stato ansatz $|\psi(\{\bar{\kappa}\}, x)\rangle$.

Contributi Chiave

• Ricetta di Prima Quantizzazione: Il documento fornisce una ricetta specifica per la simulazione relativistica utilizzando la prima quantizzazione, offrendo un'alternativa ai metodi ampiamente utilizzati di seconda quantizzazione, particolarmente vantaggiosi per sistemi con non linearità.
• Quadro Unificato per le Condizioni al Contorno: Presenta un approccio unificato per gestire le condizioni al contorno Periodiche e di Dirichlet. Fondamentalmente, dimostra che i dati ottenuti per le PBC possono essere riutilizzati per stimare le energie DBC aggiungendo specifici termini di correzione al contorno, riducendo l'overhead computazionale.
• Compatibilità con Dispositivi Near-Term: Il metodo utilizza primitivi di circuito quantistico standard (addizionatori controllati, gate di Hadamard e misurazioni) e l'ottimizzazione variazionale, rendendolo adatto all'implementazione su dispositivi quantistici near-term.
• Design Espliciti dei Circuiti: Gli autori forniscono design dettagliati dei circuiti per calcolare i valori di aspettazione degli operatori di traslazione e i termini di correzione al contorno necessari per l'espansione perturbativa.

Risultati e Rivendicazioni
Il documento sostiene di offrire una "via pratica" per simulare effetti relativistici su dispositivi quantistici near-term. I risultati teorici mostrano che:

• I termini dell'energia cinetica relativistica fino a un ordine specifico di perturbazione possono essere decomposti in combinazioni lineari di valori di aspettazione di operatori di traslazione ($\langle \hat{A}^{(l)} \rangle$) e operatori di contorno ($\langle \hat{E}_d \rangle$).
• I coefficienti per queste espansioni possono essere pre-calcolati in base ai parametri fisici ($m, c, \hbar, \delta x$).
• Il metodo consente l'ottimizzazione variazionale degli stati ansatz per stimare sia le energie dello stato fondamentale non relativistiche che quelle relativistiche.

Significato
Gli autori posizionano questo lavoro come un passo fondamentale verso il supporto dei futuri sviluppi nella fisica e nella chimica quantistica. Abilitando la simulazione di effetti relativistici su hardware near-term, il framework colma una lacuna nelle attuali capacità di simulazione quantistica. Il documento sottolinea che, sebbene l'aumento del numero di qubit migliori la risoluzione della griglia, è necessario prestare attenzione affinché il parametro adimensionale $L_m$ rimanga sufficientemente piccolo affinché l'espansione perturbativa rimanga valida. Il lavoro non pretende di risolvere l'equazione di Dirac completa o di gestire i limiti ultra-relativistici, ma si concentra specificamente sul regime in cui le correzioni relativistiche al quadro non relativistico sono rilevanti. Il lavoro futuro suggerisce di esplorare operatori Laplaciani di ordine superiore e metodi di costruzione dell'ansatz più accurati.