Digital Quantum Simulation of the quantum $\beta$-FPUT Lattice: Formulation and Resource Estimation

Sintesi Tecnica: Simulazione Quantistica Digitale del Reticolo Quantistico $\beta$-FPUT

Enunciato del Problema
Il documento affronta la sfida di simulare la conduzione termica in sistemi a bassa dimensionalità, in particolare il reticolo quantistico $\beta$-Fermi–Pasta–Ulam–Tsingou ($\beta$-FPUT). Sebbene la dinamica molecolare (MD) classica catturi la dinamica in tempo reale, trascura le statistiche quantistiche e il moto di punto zero. Al contrario, il Monte Carlo con Integrali di Percorso (PIMC) modella con precisione le proprietà quantistiche di equilibrio, ma non riesce a fornire coefficienti di trasporto in tempo reale affidabili a causa della natura mal posta dell'analisi di continuazione analitica. Gli autori sostengono che il modello quantistico $\beta$-FPUT, che descrive una catena unidimensionale di masse con interazioni armoniche e anarmoniche di quarto ordine, rappresenta un paradigma minimale per studiare il trasporto anomalo e i tempi di vita dei fononi nei materiali quantistici a bassa dimensionalità. Tuttavia, l'accesso alla sua dinamica quantistica in tempo reale richiede risorse computazionali oltre le attuali capacità classiche e l'hardware quantistico intermedio su scala rumorosa (NISQ), rendendo necessaria un'approccio di simulazione quantistica digitale tollerante ai guasti.

Metodologia
Gli autori propongono un framework di simulazione quantistica digitale in prima quantizzazione progettato per computer quantistici tolleranti ai guasti. La metodologia evita la codifica bosonica (troncamento dei numeri di occupazione dei fononi) tipicamente utilizzata negli approcci in seconda quantizzazione, che introduce errori sistematici e sovraccarichi. Invece, l'approccio discretizza direttamente gli spostamenti del reticolo su griglie finite.

1. Formulazione dell'Hamiltoniana: Il sistema è definito da un'Hamiltoniana nello spazio reale contenente energia cinetica e un termine di energia potenziale con componenti sia armoniche ($\kappa$) che anarmoniche di quarto ordine ($\beta$). Gli autori utilizzano una rappresentazione in prima quantizzazione in cui posizione ($q_j$) e momento ($p_j$) sono operatori su griglie discrete.
2. Evoluzione Temporale (Trotterizzazione): L'evoluzione in tempo reale è approssimata utilizzando una formula di prodotto simmetrica di secondo ordine di Trotter–Suzuki.
   • Operatore Cinetico: Implementato tramite una Trasformata di Fourier Quantistica (QFT) su ogni registro di sito per passare alla base del momento, applicando una rotazione di fase diagonale (rimbalzo di fase) basata su $p^2$, e tornando alla base della posizione tramite QFT inversa.
   • Operatore Potenziale: Implementato nella base della posizione. L'algoritmo calcola gli spostamenti relativi ($\Delta_j = q_{j+1} - q_j$) utilizzando una sottrazione reversibile in registri ancilla. Vengono quindi applicate rotazioni di fase controllate basate sui termini quadratici e di quarto ordine di $\Delta_j$. Il sistema è partizionato in strati di legami pari e dispari per consentire l'esecuzione parallela.
3. Estrazione degli Osservabili (Correlatori): Per misurare i correlatori di spostamento risolti per modo $\langle \hat{Q}_k(t) \hat{Q}_k(0) \rangle$, gli autori affrontano il fatto che l'operatore di modo di Fourier $\hat{Q}_k$ non è hermitiano. Decompongono $\hat{Q}_k$ in quadrature hermitiane commutanti (componenti coseno e seno). La funzione di correlazione viene ricostruita stimando le funzioni generatrici di questi operatori hermitiani utilizzando un circuito basato sul test di Hadamard. Le derivate miste sono approssimate tramite stimatori alle differenze finite per isolare i termini correlatori desiderati.

Contributi Chiave

• Framework in Prima Quantizzazione: Il documento stabilisce un protocollo di simulazione che opera direttamente sugli spostamenti discretizzati del reticolo, evitando il sovraccarico e gli errori di troncamento associati alle codifiche di modo bosonico.
• Decomposizione in Quadrature Hermitiane: Viene introdotto un protocollo innovativo per gestire gli operatori di modo di Fourier non hermitiani. Decomponendoli in quadrature hermitiane seno e coseno commutanti, gli autori permettono la costruzione di circuiti quantistici poco profondi (profondità due) per le necessarie rotazioni unitarie, facilitando una misurazione efficiente.
• Flusso di Lavoro End-to-End: Il lavoro fornisce un progetto algoritmico completo, dalla formulazione dell'Hamiltoniana alla costruzione e misurazione del circuito Trotterizzato, specificamente adattato per l'esecuzione tollerante ai guasti.
• Stima delle Risorse: Gli autori forniscono stime dettagliate per il numero di qubit, la complessità delle porte e la profondità del circuito in funzione della dimensione del sistema ($N$) e della risoluzione della griglia ($b$ bit per sito).

Risultati e Stime delle Risorse
Il documento presenta leggi di scalatura asintotiche e stime concrete delle risorse per il flusso di lavoro proposto:

• Numero di Qubit: Il requisito totale di qubit logici è stimato come $Q_{tot} = \frac{3}{2}Nb + O(1)$, dove $N$ è il numero di siti del reticolo e $b$ è il numero di qubit per sito. Questo tiene conto del registro di sistema e dei qubit ancilla richiesti per l'aritmetica reversibile e l'elaborazione parallela dei legami.
• Complessità delle Porte: Il conteggio delle porte per passo di Trotter scala come $O(Nb^2)$. Ciò deriva dalle operazioni QFT ($O(b^2)$ per sito) e dall'aritmetica reversibile richiesta per le fasi dell'energia potenziale.
• Profondità del Circuito: La profondità per passo di Trotter scala come $O(b)$, assumendo l'esecuzione parallela di interazioni di legami non sovrapposti e operazioni sui siti.
• Precisione vs. Costo: Utilizzando una formula di Trotter di secondo ordine, il numero di passi $n$ richiesto per raggiungere una precisione target $\epsilon$ scala come $n = O(t^{3/2}/\sqrt{\epsilon})$. Il costo computazionale totale scala linearmente con il tempo simulato e il numero di punti temporali campionati.

Significato e Affermazioni
Gli autori affermano che questo lavoro stabilisce un progetto algoritmico concreto per simulare la dinamica del trasporto quantistico in modelli di reticolo non lineari a bassa dimensionalità su futuri hardware quantistici tolleranti ai guasti. Sottolineano che l'approccio è esplicitamente progettato per ambienti tolleranti ai guasti, riconoscendo che i circuiti profondi richiesti per simulazioni accurate in tempo reale sono oltre le capacità dei dispositivi NISQ attuali.

Il significato del lavoro risiede nel colmare il divario tra fisica fondamentale (trasporto anomalo nel modello $\beta$-FPUT) e progettazione pratica di algoritmi quantistici. Fornendo una formulazione in prima quantizzazione efficiente in termini di risorse e un protocollo specifico per l'estrazione di correlatori risolti per modo, il documento offre una roadmap per investigare fenomeni di trasporto termico quantistico attualmente inaccessibili ai metodi classici. Gli autori notano che, sebbene il framework sia una "roadmap tollerante ai guasti", i passi futuri includono il confronto con la diagonalizzazione esatta per sistemi piccoli e il perfezionamento delle strategie di preparazione dello stato per stati iniziali termici.