Hamiltonian dynamics from pure dissipation

Titel: Hamiltoniaanse dynamica uit pure dissipatie

Auteurs: Zhong-Xia Shang en Daniel Stilck França (Universiteit van Kopenhagen)

---

1. Probleemstelling

In de kwantummechanica wordt er fundamenteel onderscheid gemaakt tussen gesloten en open systemen:

• Gesloten systemen: Evolueren reversibel onder een unitaire Hamiltoniaanse dynamica ($d\rho/dt = -i[H, \rho]$). Er is geen interactie met de omgeving en geen informatieverlies.
• Open systemen: Koppelen continu aan een externe bad, wat leidt tot irreversibele dissipatie en decoherentie, beschreven door de Lindblad-vergelijking (GKSL-vorm). Deze bevat zowel een coherente term ($H$) als dissipatieve termen (sprongoperatoren $F_i$).

De centrale vraag van dit werk is: Kan coherente Hamiltoniaanse dynamica (reversibel) worden benaderd door uitsluitend dissipatieve interacties (irreversibel), en zo ja, tegen welke kosten?

Hoewel Stinespring-dilatatie toelaat om open systemen te modelleren als gesloten systemen door de Hilbertruimte uit te breiden, blijft de vraag of men binnen het systeem zelf, zonder een expliciete Hamiltoniaan-term in de generator, Hamiltoniaanse beweging kan "faken" puur via dissipatie. Eerdere resultaten toonden aan dat een exacte gelijkheid ($\sum (F_i \rho F_i^\dagger - \frac{1}{2}\{F_i^\dagger F_i, \rho\}) = -i[H, \rho]$) onmogelijk is voor $H \neq 0$. Dit werk onderzoekt echter de benadering.

2. Methodologie en Constructie

De auteurs presenteren een constructie waarbij een zuiver dissipatieve Lindblad-generator (zonder expliciete $H$-term) Hamiltoniaanse dynamica benadert.

De Kernconstructie:
Ze beschouwen een Lindblad-vergelijking met slechts één sprongoperator $F$:
$$\frac{d\rho}{dt} = F\rho F^\dagger - \frac{1}{2}\{F^\dagger F, \rho\}$$
Door de sprongoperator te kiezen als een kleine perturbatie van de eenheidsoperator:
$$F = I - i\delta H$$
waarbij $\delta$ een kleine parameter is en $H$ de doel-Hamiltoniaan.

Wanneer men dit in de Lindblad-vergelijking substitueert, ontstaat:
$$\frac{d\rho}{dt} = -i\delta[H, \rho] + \delta^2 \left( H\rho H - \frac{1}{2}\{H^2, \rho\} \right)$$

• De eerste term ($-i\delta[H, \rho]$) is de gewenste Hamiltoniaanse dynamica met sterkte $\delta$.
• De tweede term is een dissipatieve term met sterkte $\delta^2$.

Om de totale evolutie over een tijd $t$ te simuleren met een fout $\epsilon$, kiezen de auteurs $\delta = \epsilon/t$. Omdat de dissipatie ($\delta^2$) veel kleiner is dan de coherente term ($\delta$) voor kleine $\delta$, domineert de rotatie, maar de dissipatie is onvermijdelijk.

Verwachte Evolutietijd:
Om de totale evolutie $t$ te bereiken, moet het systeem evolueren gedurende een totale tijd $T = t/\delta$. Met $\delta = \epsilon/t$ wordt de vereiste tijd:
$$T = O\left(\frac{t^2}{\epsilon}\right)$$
Dit betekent dat de simulatie tijd kost die kwadratisch is in de doel-tijd $t$ en omgekeerd evenredig met de toegestane fout $\epsilon$.

3. Belangrijkste Resultaten

A. Bovenste Grens (Theorema 1)

De auteurs bewijzen dat een zuiver dissipatieve Lindblad-generator Hamiltoniaanse dynamica $e^{-iHt}$ kan benaderen binnen een diamantafstand (diamond norm) fout $\epsilon$ met een evolutietijd $T = O(\|H\|^2 t^2 / \epsilon)$.

• Dit geldt ook voor lokale systemen: als $H$ een som is van lokale termen $H = \sum H_i$, kan men lokale sprongoperatoren $F_i = I - i\delta H_i$ gebruiken. De schaling blijft polynomiëel.
• Dit toont aan dat coherente dynamica efficiënt kan worden gesimuleerd door dissipatie zonder exponentiële overhead.

B. Onderste Grens en Optimaliteit (Theorema 2)

Een cruciale bijdrage is het bewijs dat de $O(t^2/\epsilon)$ schaling noodzakelijk en optimaal is in het slechtste geval.

• Geometrisch Argument: De auteurs analyseren de dynamica op de Bloch-sfeer (voor een qubit). Hamiltoniaanse dynamica is een rotatie (tangentiële beweging), terwijl dissipatie een radiale contractie veroorzaakt (decoherentie).
• Er bestaat een fundamentele relatie tussen de rotatiesnelheid ($\omega$) en de radiale contractie ($S$): $|\omega|^2 \leq m C^2 (-\text{Tr}(S))$.
• Om de rotatie te maximaliseren zonder dat de staat te veel van de Bloch-sfeer afwijkt (binnen fout $\epsilon$), moet de contractie klein zijn. Dit dwingt de rotatiesnelheid om laag te zijn, wat leidt tot een lange evolutietijd $T$.
• De ondergrens is: $T = \Omega(\|H\|^2 t^2 / (m C^2 \epsilon))$.
• Interpretatie: Je kunt Hamiltoniaanse snelheid niet "kopen" zonder de prijs van decoherentie te betalen. De $t^2$-factor is de fundamentele "decoherentie-kost" om de snelheid van Hamiltoniaanse dynamica bij te houden.

C. Uniekheid (Propositie 1)

Ze bewijzen dat voor elke gladde familie van zuiver dissipatieve Lindblad-generatoren die Hamiltoniaanse dynamica benadert, de sprongoperator noodzakelijkerwijs de vorm $F = I - \delta A + O(\delta^2)$ moet hebben, waarbij het anti-Hermietse deel van $A$ gelijk is aan $H$. Dit bevestigt dat hun constructie de enige mogelijke pad is voor dergelijke benaderingen.

4. Implicaties en Toepassingen

1. BQP-Completteit: Zuiver dissipatieve dynamica is BQP-compleet (voor berekeningen die niet alleen het steady-state bereiken, maar de dynamica zelf). Dit betekent dat zuivere dissipatie universeel kwantumcomputen kan uitvoeren, zelfs zonder een expliciete Hamiltoniaan, zolang men de dynamica onderweg kan lezen. Dit verschilt van eerdere resultaten die zich richtten op het steady-state.
2. Nieuw Zeno-effect: De auteurs beschrijven een "Zeno-adjacent" bevriezingseffect. Door een sprongoperator $F = \delta^{-1/2}(I + i\delta H)$ in te voeren, kan men de interne evolutie $H$ exact compenseren via dissipatie, waardoor elke toestand wordt bevroren (in tegenstelling tot het traditionele Zeno-effect dat alleen toestanden in een Zeno-ruimte bevriest).
3. Geen Super-kwadratische "Fast-Forwarding": Voor deze klasse van dissipatieve Lindblad-systemen geldt dat ze niet sneller dan kwadratisch kunnen worden "fast-forwarded" (versneld). Als ze dat wel zouden kunnen, zou dit in strijd zijn met de query-ondergrenzen voor het parity-check probleem. Dit suggereert een nieuwe klasse van Lindblad-systemen die kwadratisch versneld kunnen worden, maar niet sneller.
4. Vermindering van Simulatiekosten via "Gauge Changing": De auteurs tonen aan dat de simulatiekosten op een kwantumcomputer afhangen van de normen van de Hamiltoniaan en sprongoperatoren. Door een "gauge-transformatie" toe te passen (het verschuiven van termen tussen $H$ en $F_i$), kan men de normen verkleinen. Een Lindblad-systeem dat lijkt duur om te simuleren, kan in een andere "gauge" (met een andere verdeling van $H$ en $F$) veel goedkoper worden gesimuleerd, zelfs als de fysieke dynamica identiek blijft.

5. Conclusie en Significatie

Dit werk legt een fundamentele brug tussen twee ogenschijnlijk tegenstrijdige concepten: reversibele Hamiltoniaanse dynamica en irreversibele dissipatie.

• Fysische Inzicht: Het toont aan dat coherente rotatie kan worden "geënsceneerd" door externe dissipatie, maar dat dit een fundamentele prijs heeft: de tijd die nodig is om de rotatie te bereiken, schaalt kwadratisch met de doel-tijd. Dit is een directe manifestatie van de onmogelijkheid om reversibiliteit te creëren zonder decoherentie in een open systeem.
• Technologische Impact: De resultaten zijn relevant voor het ontwerp van dissipatieve kwantumalgoritmen (zoals voor Gibbs-state voorbereiding en optimalisatie). Het suggereert dat men niet noodzakelijk een complexe Hamiltoniaan hoeft te implementeren om coherente effecten te bereiken; zuivere dissipatie kan volstaan, mits men rekening houdt met de tijdschaal.
• Complexiteitstheorie: Het plaatst dissipatieve dynamica in de complexiteitsklasse BQP, wat betekent dat het een krachtig model is voor kwantumcomputatie, zelfs zonder traditionele unitaire poorten.

Samenvattend bewijzen de auteurs dat "dissipatie Hamiltoniaanse dynamica kan faken", maar dat dit faken een fundamentele, kwadratische tijdsstraf vereist die voortkomt uit de geometrie van decoherentie.