Quantum Simulation of Two-Level $PT$-Symmetric Systems Using Hermitian Hamiltonians

Sintesi Tecnica: Simulazione Quantistica di Sistemi a Due Livelli con Simmetria PT Utilizzando Hamiltoniani Ermitici

Definizione del Problema
Simulare sistemi quantistici aperti (OQS) su computer quantistici basati su gate unitari presenta una sfida fondamentale, poiché gli OQS spesso esibiscono propagatori temporali non unitari, che non possono essere mappati direttamente nel framework unitario dell'attuale hardware quantistico. I sistemi con simmetria Parità-Tempo (PT) sono una specifica classe di OQS caratterizzati da Hamiltoniani non ermitici che, nel loro regime "non rotto", possiedono spettri di energia reali. Sebbene questi sistemi offrano proprietà uniche per il sensing quantistico e la generazione di entanglement — come tempi di evoluzione tra stati più rapidi rispetto ai sistemi ermitici — la loro dinamica non unitaria complica la simulazione diretta sui processori quantistici. I metodi esistenti per la simulazione di sistemi aperti, inclusi l'evoluzione nel tempo immaginario e il block encoding, richiedono approcci efficienti specificamente adattati alla dinamica PT-simmetrica per sfruttare i loro vanti unici.

Metodologia
Gli autori propongono un metodo per simulare la dinamica PT-simmetrica sfruttando le trasformazioni di similitudine per mappare gli Hamiltoniani non ermitici PT-simmetrici su equivalenti Hamiltoniani ermitici. Questa mappatura è valida nel regime di simmetria PT non rotto, dove l'Hamiltoniano è pseudo-ermitico. Il nucleo dell'approccio consiste in una matrice di trasformazione di similitudine, $\tau = \sqrt{\eta}$, che connette l'Hamiltoniano non ermitico PT-simmetrico ($\hat{H}_{PT}$) a un equivalente ermitico ($\hat{h}$). L'evoluzione temporale è quindi espressa come:
$$e^{-i\hat{H}_{PT}t} = \tau^{-1} e^{-i\hat{h}t} \tau$$
Questa formulazione permette alla simulazione di essere decomposta in tre fasi: inizializzazione tramite $\tau$, evoluzione unitaria tramite $e^{-i\hat{h}t}$, e una trasformazione finale tramite $\tau^{-1}$.

Il documento introduce due algoritmi distinti per implementare questo framework:

1. Algoritmo Ibrido Classico-Quantistico: Questo approccio esegue l'inizializzazione e l'evoluzione temporale unitaria su un singolo qubit. L'ultimo passaggio, l'applicazione di $\tau^{-1}$, è gestito tramite post-processing classico. Ciò richiede la tomografia completa dello stato quantistico per ricostruire la matrice di densità prima di eseguire la moltiplicazione classica per $\tau^{-1}$.
2. Algoritmo di Dilatazione: Per evitare i limiti di scalabilità della tomografia completa, questo metodo implementa $\tau^{-1}$ come un operatore unitario dilatato utilizzando un qubit ancilla. Utilizzando la decomposizione spettrale della matrice ermitica $\tau$, l'operatore non unitario viene incorporato in uno spazio unitario più ampio. L'evoluzione è condizionata dal fatto che il qubit ancilla sia misurato nello stato $|0\rangle$. Questo approccio elimina la necessità di una tomografia completa della matrice di densità.

Per estendere questi metodi a sistemi più grandi, gli autori applicano la teoria delle perturbazioni del primo ordine per simulare due sistemi PT-simmetrici debolmente interagenti. Costruiscono una nuova trasformazione di similitudine basata sugli autostati perturbati del sistema accoppiato, permettendo alla dinamica di essere mappata su un'evoluzione ermitica adatta all'algoritmo ibrido.

Risultati Chiave
Gli algoritmi sono stati validati attraverso simulazioni ed esperimenti sulla QPU IBM-Sherbrooke e su simulatori rumorosi:

• Dinamica di un Singolo Qubit: Entrambi gli algoritmi hanno riprodotto con successo la dinamica delle popolazioni dello stato fondamentale e dello stato eccitato di un sistema PT-simmetrico vicino a un punto eccezionale. I dati sperimentali (cerchi) e le simulazioni del modello di rumore (linee tratteggiate) mostrano un forte accordo con le soluzioni classiche esatte (linee continue).
• Entanglement a Due Qubit: Utilizzando l'algoritmo ibrido e la teoria delle perturbazioni, gli autori hanno simulato due qubit PT-simmetrici debolmente accoppiati. Hanno dimostrato con successo la generazione di entanglement, quantificato dalla concordanza, che ha raggiunto un massimo quando le popolazioni degli stati di base sono diventate approssimativamente uguali.
• Analisi dell'Errore: Lo studio ha analizzato la fedeltà di entrambi gli algoritmi sotto variabili errori di misura.
  • L'algoritmo ibrido è risultato altamente sensibile agli errori di misura, particolarmente nel regime di evoluzione rapida dove gli errori sono amplificati dal post-processing classico. La fedeltà è scesa significativamente (fino a ~20%) ad alti tassi di errore.
  • L'algoritmo di dilatazione, pur soffrendo dell'aumento del volume del circuito e degli errori dei gate a due qubit, ha dimostrato una maggiore robustezza contro gli errori di misura, mantenendo una fedeltà più elevata (~70%) ad alti tassi di errore rispetto all'approccio ibrido. Tuttavia, rimane sensibile ai limiti di coerenza inerenti ai gate a due qubit.

Significatività e Rivendicazioni
Il documento rivendica di offrire un approccio scalabile per simulare sistemi PT-simmetrici e operatori pseudo-ermitici su computer quantistici a gate unitari. Costruendo Hamiltoniani ermitici equivalenti, gli autori dimostrano che la dinamica PT-simmetrica non unitaria può essere simulata usando l'evoluzione in tempo reale piuttosto che singoli gate non unitari.

Il lavoro evidenzia che:

• I sistemi PT-simmetrici possono essere simulati su dispositivi quantistici NISQ (Noisy Intermediate-Scale Quantum) attuali utilizzando sia approcci ibridi che pienamente quantistici (dilatazione).
• I metodi possono essere estesi a sistemi multi-qubit debolmente accoppiati tramite la teoria delle perturbazioni, consentendo lo studio della generazione di entanglement in sistemi non ermitici.
• Le caratteristiche uniche dell'errore della simulazione PT-simmetrica — specificamente l'amplificazione degli errori nei regimi di evoluzione rapida — rispecchiano il comportamento fisico unico di questi sistemi vicino ai punti eccezionali, suggerendo un potenziale utilità per applicazioni di sensing quantistico.

Gli autori concludono che questo lavoro apre una strada per ulteriori investigazioni sugli Hamiltoniani non ermitici sui computer quantistici, fornendo una base per lavori futuri volti a migliorare l'efficienza algoritmica ed esplorare le applicazioni pratiche di questi sistemi non ermitici.