Quantum Simulation of Non-Hermitian Linear Response

Titel: Quantum Simulatie van Niet-Hermitiese Lineaire Respons

Auteur: Jeongbin Jo (Yonsei University, Seoel)
Datum: 19 maart 2026

---

1. Het Probleem

Lineaire respons-theorie en Green-functies vormen de hoeksteen voor het begrijpen van hoe macroscopische en sterk gecorreleerde systemen reageren op zwakke externe verstoringen. Hoewel de theoretische basis voor niet-Hermitiese lineaire respons recent is gelegd om open kwantumsystemen (met dissipatie of continue metingen) te beschrijven, blijft het vertalen van deze theorie naar praktische quantumcomputers een enorme uitdaging.

De kern van het probleem ligt in de niet-unitaire aard van de dynamica in open systemen. De tijdevolutie wordt bepaald door een superoperator $L$ (de Lindbladiaan), waarbij de evolutie-operator $e^{Lt}$ niet-unitair is.

• Traditionele quantumalgoritmen vereisen unitaire operatoren.
• Bestaande methoden om niet-unitaire dynamica te simuleren (zoals het gebruik van gedilateerde lineaire systeemoplossers of spectrale mapping) leiden tot exponentiële overhead in de diepte van de schakeling en de kosten voor het voorbereiden van de toestand. Dit maakt ze onpraktisch voor grootschalige simulaties.

2. Methodologie: Schrödingerisatie

De auteur presenteert een systematische algoritmische mapping die niet-unitaire multi-tijd correlatiefuncties transformeert in een unitaire vorm die uitvoerbaar is op quantumhardware. De centrale techniek is Schrödingerisatie (gebaseerd op werk van Jin et al.).

De methode verloopt als volgt:

1. Vectorisatie: De Lindblad-mastervergelijking voor de dichtheidsmatrix $\rho(t)$ wordt omgezet naar een lineaire differentiaalvergelijking in de Liouville-ruimte: $\frac{d}{dt}|\rho(t)\rangle = L |\rho(t)\rangle$.
2. Decompositie: De niet-Hermitiese generator $L$ wordt ontbonden in een Hermitiese en een anti-Hermitiese component: $L = H_1 - iH_2$.
3. Gedwongen Fase-transformatie: Er wordt een reële continue variabele $\xi > 0$ geïntroduceerd met de transformatie $w(t, \xi) = e^{-\xi} |\rho(t)\rangle$.
4. Fourier-transformatie: Door een Fourier-transformatie toe te passen op de variabele $\xi$ (naar de modus $\eta$), wordt het systeem omgezet in een reeks ongekoppelde Schrödinger-achtige vergelijkingen:
   $$i\partial_t \tilde{w} = (\eta H_1 + H_2) \tilde{w}$$
   De nieuwe effectieve Hamiltoniaan $H_{sch}(\eta) = \eta H_1 + H_2$ is strikt Hermities voor elke reële $\eta$. Dit stelt de quantumcomputer in staat om de evolutie unitair te simuleren via $U(t, \eta) = e^{-iH_{sch}(\eta)t}$.
5. Herconstructie: De fysieke dynamica wordt hersteld door een inverse Fourier-transformatie over de $\eta$-modi uit te voeren.

3. Belangrijkste Bijdragen

• Brug tussen Theorie en Implementatie: Het artikel sluit de kloof tussen de gevestigde fysica van niet-Hermitiese lineaire respons en kwantumcomputing. Het biedt het eerste schaalbare algoritmische pad om dissipatieve respons te simuleren zonder de exponentiële kosten van traditionele methoden.
• Generalisatie van Kubo-formule: De methode stelt in staat om de Kubo-formule voor Green-functies en responsfuncties $\chi(\tau)$ te evalueren in dissipatieve, niet-Hermitiese regimes.
• Optimale Kosten voor Toestandsvoorbereiding: In tegenstelling tot methoden zoals QSVT (Quantum Singular Value Transformation) of HHL die complexe blokken-codering vereisen voor niet-unitaire matrices, vereist Schrödingerisatie alleen de simulatie van de Hermitiese delen $H_1$ en $H_2$. Dit resulteert in optimale kosten voor het voorbereiden van de toestand.
• Analytische Validatie: De auteurs leiden de mapping expliciet af voor een fundamenteel open systeem: een enkele qubit met amplitude-demping (spontane emissie).

4. Resultaten

De auteurs hebben de methode gevalideerd via numerieke simulaties op een quantumcircuit (gebaseerd op Qiskit Sampler) voor het amplitude-dempingsmodel:

• Convergentie: De simulatie toont perfecte convergentie naar de exacte oplossing van de Lindblad-mastervergelijking wanneer de resolutie wordt verhoogd (hoge $N$ voor Fourier-grid en hoge $L$ voor truncatie).
• Foutanalyse:
  • De fout door truncatie van het Fourier-domein ($L$) schaalt exponentieel: $\epsilon \sim e^{-L}$.
  • De fout door discretisatie ($N$) neemt exponentieel af met het aantal qubits.
• Complexiteit: De schaalbaarheid is $O(\text{poly}(\log(1/\epsilon)))$ met betrekking tot de precisie $\epsilon$. Dit is een aanzienlijke verbetering ten opzichte van gedilateerde methoden die $O(\text{poly}(1/\epsilon))$ vereisen.
• Circuit-ontwerp: Er werd een aangepaste Hadamard-test-circuit ontworpen om het reële deel van de overlap van de dichtheidsmatrix te extraheren, waarbij een hulpqubit wordt gebruikt om de unitaire evolutie te controleren.

5. Betekenis en Toekomstperspectief

Dit werk is van fundamenteel belang voor de ontwikkeling van quantumcomputing in de studie van gedreven-dissipatieve systemen.

• Toepassingsgebied: Het opent de deur voor het simuleren van sterk gecorreleerde elektronensystemen en fysische chemie in open omgevingen, waar dissipatie een cruciale rol speelt.
• Efficiëntie: Door de noodzaak van complexe matrixinversie-algoritmen te omzeilen, biedt de methode een haalbare route voor het simuleren van hoge-dimensionale open systemen op toekomstige quantumhardware.
• Uitdagingen: Hoewel de theorie robuust is, wijzen de auteurs erop dat de implementatie op nabije-toekomstige hardware (NISQ) uitdagingen kent door de gevoeligheid van de inverse Fourier-transformatie voor ruis en de diepte van de schakeling.

Kortom, dit artikel introduceert een krachtige, unitaire "wrapper" voor niet-unitaire dynamica, waardoor het mogelijk wordt om lineaire respons in open kwantumsystemen efficiënt en nauwkeurig te simuleren op quantumcomputers.