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

Technische Samenvatting: Kwantumsimulatie van Twee-Niveau PT-Symmetrische Systemen met behulp van Hermitische Hamiltonia

Probleemstelling
Het simuleren van open kwantumsystemen (OQS'en) op unitair-gate-gebaseerde kwantumcomputers vormt een fundamentele uitdaging omdat OQS'en vaak niet-unitaire tijdspropagatoren vertonen, die niet direct naar het unitaire kader van de huidige kwantumhardware kunnen worden gemapt. Pariteit-Tijd (PT)-symmetrische systemen zijn een specifieke klasse van OQS'en die gekenmerkt worden door niet-Hermitische Hamiltonia die, in hun "ongebroken" fase, een reële energiespectra bezitten. Hoewel deze systemen unieke eigenschappen bieden voor kwantumsensing en verstrengelingsgeneratie—zoals snellere evolutietijden tussen toestanden vergeleken met Hermitische systemen—compliceren hun niet-unitaire dynamica de directe simulatie op kwantumprocessors. Bestaande methoden voor open kwantumsimulatie, waaronder imaginaire tijdevolutie en block-encodings, vereisen efficiënte benaderingen die specifiek zijn afgestemd op PT-symmetrische dynamica om hun unieke voordelen te benutten.

Methodologie
De auteurs stellen een methode voor om PT-symmetische dynamica te simuleren door gebruik te maken van gelijkenistransformaties om niet-Hermitische PT-symmetrische Hamiltonia te mappen op equivalente Hermitische Hamiltonia. Deze mapping is geldig in het ongebroken PT-symmetrische regime waar het Hamiltoniaan pseudo-Hermitisch is. De kern van de aanpak betreft een gelijkenistransformatie-matrix, $\tau = \sqrt{\eta}$, die de niet-Hermitische PT-symmetrische Hamiltoniaan ($\hat{H}_{PT}$) verbindt met een Hermitische equivalente ($\hat{h}$). De tijdevolutie wordt dan uitgedrukt als:
$$e^{-i\hat{H}_{PT}t} = \tau^{-1} e^{-i\hat{h}t} \tau$$
Deze formulering maakt het mogelijk om de simulatie te ontleden in drie stappen: initialisatie via $\tau$, unitaire evolutie via $e^{-i\hat{h}t}$, en een uiteindelijke transformatie via $\tau^{-1}$.

Het artikel introduceert twee onderscheidende algoritmen om dit kader te implementeren:

1. Hybride Klassiek-Kwantum Algoritme: Deze benadering voert de initialisatie en de unitaire tijdevolutie uit op een enkele qubit. De laatste stap, het toepassen van $\tau^{-1}$, wordt afgehandeld via klassieke post-verwerking. Dit vereist volledige kwantumtoestands-tomografie om de dichtheidsmatrix te reconstrueren voordat de klassieke vermenigvuldiging met $\tau^{-1}$ wordt uitgevoerd.
2. Dilatie Algoritme: Om de schaalbaarheidsbeperkingen van volledige tomografie te vermijden, implementeert deze methode $\tau^{-1}$ als een gedilateerde unitaire operator met behulp van een ancilla-qubit. Door gebruik te maken van de spectrale decompositie van de Hermitische matrix $\tau$, wordt de niet-unitaire operator ingebed in een grotere unitaire ruimte. De evolutie is geconditioneerd op het feit dat de ancilla-qubit wordt gemeten in de $|0\rangle$-toestand. Deze benadering elimineert de noodzaak voor volledige dichtheidsmatrix-tomografie.

Om deze methoden uit te breiden naar grotere systemen, passen de auteurs eerste-orde perturbatietheorie toe om twee zwak interagerende PT-symmetrische systemen te simuleren. Zij construeren een nieuwe gelijkenistransformatie gebaseerd op de geperturbeerde eigenfuncties van het gekoppelde systeem, waardoor de dynamica kan worden gemapt op een Hermitische evolutie die geschikt is voor het hybride algoritme.

Belangrijkste Resultaten
De algoritmen werden gevalideerd door middel van simulaties en experimenten op de IBM-Sherbrooke kwantumprocessor (QPU) en ruismodellen:

• Enkele Qubit Dynamica: Beide algoritmen reproduceerden succesvol de grondtoestand- en aangeslagen-toestand populatiedynamica van een enkel PT-symmetrisch systeem nabij een exceptional point. De experimentele data (cirkels) en ruismodel-simulaties (onderbroken lijnen) vertoonden een sterke overeenkomst met de exacte klassieke oplossingen (doorgetrokken lijnen).
• Twee-Qubit Verstrengeling: Met behulp van het hybride algoritme en perturbatietheorie simuleerden de auteurs twee zwak gekoppelde PT-symmetrische qubits. Zij demonstreerden succesvol de generatie van verstrengeling, gekwantificeerd door concurrence, die een maximum bereikte wanneer de populaties van de basis-toestanden ongeveer gelijk waren.
• Foutanalyse: De studie analyseerde de getrouwheid (fidelity) van beide algoritmen onder variërende meetfouten.
  • Het hybride algoritme bleek zeer gevoelig voor meetfouten, met name in het regime van snelle evolutie waarbij fouten worden versterkt door de klassieke post-verwerkingsstap. De getrouwheid daalde significant (tot ~20%) bij hoge foutensnelheden.
  • Het dilatie algoritme, hoewel het lijdt onder een verhoogd circuitvolume en twee-qubit gate-fouten, demonstreerde een grotere robuustheid tegen meetfouten, waarbij een hogere getrouwheid (~70%) werd behouden bij hoge foutensnelheden vergeleken met de hybride benadering. Het blijft echter gevoelig voor coherentiebeperkingen die inherent zijn aan twee-qubit gates.

Betekenis en Claims
Het artikel claimt een schaalbare benadering te bieden voor het simuleren van PT-symmetrische systemen en pseudo-Hermitische operatoren op unitaire gate-gebaseerde kwantumcomputers. Door equivalente Hermitische Hamiltonia te construeren, demonstreren de auteurs dat niet-unitaire PT-symmetrische dynamica gesimuleerd kan worden met behulp van real-time evolutie in plaats van individuele niet-unitaire gates.

Het werk benadrukt dat:

• PT-symmetrische systemen op huidige NISQ-apparaten (Noisy Intermediate-Scale Quantum) gesimuleerd kunnen worden met behulp van ofwel hybride ofwel volledig kwantum (dilatie) benaderingen.
• De methoden kunnen worden uitgebreid naar zwak gekoppelde multi-qubit systemen via perturbatietheorie, wat de studie van verstrengelingsgeneratie in niet-Hermitische systemen mogelijk maakt.
• De unieke foutkarakteristieken van PT-symmetrische simulaties—specifiek de amplificatie van fouten in snelle evolutieregimes—de unieke fysieke gedragingen van deze systemen nabij exceptional points weerspiegelen, wat wijst op potentieel nut voor kwantumsensing-toepassingen.

De auteurs concluderen dat dit werk een weg vrijmaakt voor verdere investigatie van niet-Hermitische Hamiltonia op kwantumcomputers, en een fundament legt voor toekomstig werk gericht op het verbeteren van algoritmische efficiëntie en het verkennen van praktische toepassingen van deze niet-Hermitische systemen.