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

Technische Zusammenfassung: Quantensimulation von PT-symmetrischen Zwei-Niveau-Systemen mittels Hermitescher Hamiltonoperatoren

Problemstellung
Die Simulation offener Quantensysteme (OQSs) auf unitären gatterbasierten Quantencomputern stellt eine fundamentale Herausforderung dar, da OQSs oft nicht-unitäre Zeit-Propagatoren aufweisen, die nicht direkt auf das unitäre Framework aktueller Quantenhardware abgebelt werden können. PT-symmetrische Systeme sind eine spezifische Klasse von OQSs, die durch nicht-Hermitesche Hamiltonoperatoren charakterisiert sind, welche sich in ihrer „ungebrochenen“ Phase reelle Energiespektren besitzen. Während diese Systeme einzigartige Eigenschaften für das Quantensensing und die Erzeugung von Verschränkung bieten – wie etwa schnellere Evolutionszeiten zwischen Zuständen im Vergleich zu Hermiteschen Systemen –, erschwert ihre nicht-unitäre Dynamik die direkte Simulation auf Quantenprozessoren. Bestehende Methoden für die Simulation offener Quantensysteme, einschließlich der Imaginärzeit-Evolution und Block-Kodierungen, erfordern effiziente Ansätze, die speziell auf die PT-symmetrische Dynamik zugeschnitten sind, um deren einzigartige Vorteile zu nutzen.

Methodik
Die Autoren schlagen eine Methode zur Simulation der PT-symmetrischen Dynamik vor, indem sie Ähnlichkeitstransformationen nutzen, um nicht-Hermitesche PT-symmetrische Hamiltonoperatoren auf äquivalente Hermitesche Hamiltonoperatoren abzubilden. Diese Abbildung ist im Bereich der ungebrochenen PT-Symmetrie gültig, in dem der Hamiltonoperator pseudo-hermitisch ist. Der Kern des Ansatzes ist eine Ähnlichkeitstransformationsmatrix, $\tau = \sqrt{\eta}$, die den nicht-Hermiteschen PT-symmetrischen Hamiltonoperator ($\hat{H}_{PT}$) mit einem Hermiteschen Äquivalent ($\hat{h}$) verbindet. Die Zeitentwicklung wird dann als:
$$e^{-i\hat{H}_{PT}t} = \tau^{-1} e^{-i\hat{h}t} \tau$$
ausgedrückt. Diese Formulierung ermöglicht es, die Simulation in drei Schritte zu zerlegen: Initialisierung via $\tau$, unitäre Evolution via $e^{-i\hat{h}t}$ und eine abschließende Transformation via $\tau^{-1}$.

Das Paper führt zwei unterschiedliche Algorithmen zur Implementierung dieses Frameworks ein:

1. Hybrider klassisch-quantenmechanischer Algorithmus: Dieser Ansatz führt die Initialisierung und die unitäre Zeitentwicklung auf einem einzelnen Qubit durch. Der letzte Schritt, das Anwenden von $\tau^{-1}$, wird über klassische Post-Processing-Verfahren gehandhabt. Dies erfordert eine vollständige Quantenzustandstomographie, um die Dichtematrix vor der klassischen Multiplikation mit $\tau^{-1}$ zu rekonstruieren.
2. Dilationsalgorithmus: Um die Skalierbarkeitsbeschränkungen der vollständigen Tomographie zu umgehen, implementiert diese Methode $\tau^{-1}$ als dilatierten unitären Operator unter Verwendung eines Ancilla-Qubits. Durch Nutzung der Spektralzerlegung der Hermiteschen Matrix $\tau$ wird der nicht-unitäre Operator in einen größeren unitären Raum eingebettet. Die Evolution ist darauf konditioniert, dass das Ancilla-Qubit im Zustand $|0\rangle$ gemessen wird. Dieser Ansatz eliminiert die Notwendigkeit einer vollständigen Dichtematrix-Tomographie.

Um diese Methoden auf größere Systeme zu erweitieren, wenden die Autoren die Störungstheorie erster Ordnung an, um zwei schwach wechselwirkende PT-sente Systeme zu simulieren. Sie konstruieren eine neue Ähnlichkeitstransformation basierend auf den gestörten Eigenzuständen des gekoppelten Systems, was die Dynamik auf eine Hermitesche Evolution ermöglicht, die für den hybriden Algorithmus geeignet ist.

Zentrale Ergebnisse
Die Algorithmen wurden durch Simulationen und Experimente auf der IBM-Sherbrooke Quantenprozessoreinheit (QPU) sowie auf verrauschten Simulatoren validiert:

• Einzel-Qubit-Dynamik: Beide Algorithmen reproduzierten erfolgreich die Grund- und angeregten Zustandspopulationen eines einzelnen PT-symmetrischen Systems nahe einem exzeptionellen Punkt. Die experimentellen Daten (Kreise) und die Rauschmodell-Simulationen (gestrichelte Linien) zeigten eine starke Übereinstimmung mit den exakten klassischen Lösungen (durchgezogene Linien).
• Zwei-Qubit-Verschränkung: Unter Verwendung des hybriden Algorithmus und der Störungstheorie simulierten die Autoren zwei schwach gekoppelte PT-symmetrische Qubits. Sie demonstrierten erfolgreich die Erzeugung von Verschränkung, quantifiziert durch die Konkurrenz (Concurrence), welche ein Maximum erreichte, wenn die Populationen der Basiszustände annähernd gleich groß wurden.
• Fehleranalyse: Die Studie analysierte die Fidelität beider Algorithmen unter variierenden Messfehlern.
  • Der hybride Algorithmus erwies sich als hochgradig sensitiv gegenüber Messfehlern, insbesondere im Regime schneller Evolution, in dem Fehler durch den klassischen Post-Processing-Schritt verstärkt werden. Die Fidelität sank signifikant (auf ~20%) bei hohen Fehlerraten.
  • Der Dilationsalgorithmus, obwohl er unter erhöhtem Schaltungsvolumen und Zwei-Qubit-Gatterfehlern leidet, zeigte eine größere Robustheit gegenüber Messfehlern und behielt bei hohen Fehlerraten eine höhere Fidelität (~70%) im Vergleich zum hybriden Ansatz bei. Er bleibt jedoch sensitiv gegenüber Kohärenzlimitierungen, die inhärent an Zwei-Qubit-Gattern sind.

Bedeutung und Ansprüche
Das Paper beansprucht, einen skalierbaren Ansatz für die Simulation von PT-symmetrischen Systemen und pseudo-hermiteschen Operatoren auf unitären gatterbasierten Quantencomputern anzubieten. Durch die Konstruktion äquivalenter Hermitescher Hamiltonoperatoren zeigen die Autoren, dass nicht-unitäre PT-symmetrische Dynamiken unter Verwendung von Echtzeit-Evolution statt einzelner nicht-unitärer Gatter simuliert werden können.

Die Arbeit hebt hervor, dass:

• PT-symmetrische Systeme auf aktuellen NISQ-Geräten (Noisy Intermediate-Scale Quantum) sowohl durch hybride als auch durch voll-quantenmechanische (Dilations-) Ansätze simuliert werden können.
• Die Methoden mittels Störungstheorie auf schwach gekoppelte Multi-Qubit-Systeme erweiterbar sind, was die Untersuchung der Verschränkungserzeugung in nicht-hermiteschen Systemen ermöglicht.
• Die einzigartigen Fehlercharakteristika von PT-Symmetrie-Simulationen – insbesondere die Verstärkung von Fehlern in schnellen Evolutionsregimen – die einzigartigen physikalischen Verhaltensweisen dieser Systeme nahe exzeptioneller Punkte widerspiegeln, was auf einen potenziellen Nutzen für Quantensensing-Anwendungen hindeutet.

Die Autoren kommen zu dem Schluss, dass diese Arbeit einen Weg für weitere Untersuchungen nicht-hermitescher Hamiltonoperatoren auf Quantencomputern ebnet und eine Grundlage für zukünftige Arbeiten zur Verbesserung der algorithmischen Effizienz und zur Erforschung praktischer Anwendungen dieser nicht-hermiteschen Systeme schafft.