Restoring detailed balance in non-Hermitian Markov processes

Technische Samenvatting: Herstel van Gedetailleerd Evenwicht in Niet-Hermitische Markov-processen

Probleemstelling
Stochastische processen buiten evenwicht, met name die gemodelleerd worden door Markov-processen met discrete toestanden, omvatten vaak asymmetrische overgangssnelheden die gedetailleerd evenwicht doorbreken. Deze niet-hermitische stochastische generatoren ($H$) leiden tot een niet-monotone entropieproductie, waarbij de entropie na een transiënte toename kan afnemen. Dit gedrag ondermijnt standaard thermodynamische interpretaties en inferentietechnieken die steunen op entropiemaximalisatie of de aanname van een altijd toenemende onzekerheid. Bovendien bemoeilijkt het gebrek aan tijdsomkeerbaarheid in deze systemen de fysieke interpretatie van evenwichtstoestanden, die configuraties met lage entropie kunnen bevatten. Hoewel tensornetwerkmethoden de "vloek van dimensionaliteit" in hoogdimensionale interagerende systemen aanpakken, lossen zij het fundamentele probleem van niet-hermitische dynamica en de resulterende niet-monotone entropie niet op.

Methodologie
De auteurs stellen het gebruik van Dyson-kaarten ($\eta$) voor om niet-hermitische stochastische generatoren te transformeren naar hermitische operatoren, terwijl het reële eigenwaardenspectrum behouden blijft. De kernmethodologie omvat:

1. Dyson-transformatie: Een omkeerbare transformatie wordt toegepast op de meestervergelijking $\partial_t |P(t)\rangle = -H|P(t)\rangle$, wat resulteert in een getransformeerde vergelijking $\partial_t |\phi(t)\rangle = -\tilde{H}|\phi(t)\rangle$, waarbij $|\phi(t)\rangle = \eta|P(t)\rangle$ en $\tilde{H} = \eta H \eta^{-1}$.
2. Metrische Operator: De transformatie introduceert een metrische operator $\Omega = \eta^\dagger \eta$, die bijdragen in het oorspronkelijke systeem herweegt. Deze operator verbindt metingen in de getransformeerde coördinaten terug met de oorspronkelijke waarschijnlijkheidsvector, waardoor statistische gemiddelden $\langle O \rangle$ kunnen worden hersteld via $\langle \phi | \tilde{\Omega} \Xi' O' | \phi \rangle$.
3. Entropieherstel: In het getransformeerde systeem wordt bewezen dat de Rényi-entropie $S' = -\ln |\phi(t)|^2$ monotoon toeneemt over de tijd ($\frac{1}{2}\frac{d}{dt}S'(t) \geq 0$). Dit staat in contrast met het oorspronkelijke systeem, waar de gekwadrateerde norm $|P(t)|^2$ (en dus de entropie) lokale minima kan vertonen.
4. Computationele Constructie:
   • Analytische Benadering: De auteurs zoeken naar een hermitische generator $\Lambda = \ln \eta$ binnen een gesloten operatoralgebra. Door de conditie $H = H^\dagger$ uit te breiden in termen van commutatoren, leiden zij een reeks algebraïsche vergelijkingen af (met behulp van coëfficiënten $G_{ijk}$) om $\Lambda$ op te lossen. Deze benadering gaat ervan uit dat er een gesloten structuur bestaat onder de commutator, wat niet noodzakelijkerwijs het geval is voor alle niet-triviale interagerende problemen.
   • Numerieke Benadering: Voor algemene gevallen waar analytische afsluiting moeilijk is, wordt een iteratief numeriek algoritme voorgesteld. Deze methode behandelt $\Lambda$ als een rotatie-as in de algebra en voert iteratieve rotaties uit om de anti-hermitische component ($\Delta H$) van de generator te onderdrukken totdat een tolerantie is bereikt.

Belangrijkste Resultaten
De methodologie werd toegepast op een Susceptible-Infected-Susceptible (SIS) epidemie-model op een statische contactgrafiek:

• Monotone Entropiegroei: In het getransformeerde systeem neemt de Rényi-entropie $S'$ monotoon toe totdat evenwicht is bereikt, terwijl de entropie van het oorspronkelijke systeem een transiënte maximum vertoont gevolgd door een afname.
• Statistische Equivalentie: Ondanks de transformatie komen de statistische gemiddelden (gemiddelde en standaarddeviatie van infectieuze knooppunten) berekend in het getransformeerde systeem met hoge precisie overeen met de resultaten van het oorspronkelijke systeem (Kolmogorov-Smirnov afstand $< 10^{-3}$), mits de metrische operator correct wordt toegepast.
• Detectie van Faseovergangen: Het maximum van de Rényi-entropie in zowel het oorspronkelijke als het getransformeerde systeem valt samen op dezelfde kritieke ratio van infectie tot herstel ($\beta/\gamma$). Dit bevestigt dat het getransformeerde systeem, ondanks het herstelde tijdsomkeerbaarheid, de signatuur van de faseovergang uit het oorspronkelijke niet-hermitische systeem behoudt.
• Structuurbehoud: De evenwichtstoestandsverdelingen in het getransformeerde systeem ($|\phi\rangle$) behouden de structurele indices van de oorspronkelijke verdeling ($|P\rangle$), wat suggereert dat de diagonale bijdragen van de Dyson-kaart de transformatie beheersen, terwijl de off-diagonale componenten dienen om componenten te mengen om de entropie te verhogen.

Betekenis en Claims
Het artikel claimt een algemene en computationeel hantebare methode te bieden om gedetailleerd evenwicht te herstellen in een brede klasse van Markov-systemen. Door niet-hermitische generatoren om te zetten in hermitische, voorziet de benadering deze systemen van een "vertrouwde fysieke interpretatie" die gekenmerkt wordt door strikt positieve entropieproductie. Dit maakt de toepassing van standaard statistische natuurkunde-instrumenten mogelijk, zoals entropiemaximalisatie, om evenwichtstoestanden te bepalen.

De auteurs benadrukken dat dit formalisme een verbinding legt tussen stochastische processen en kwantumdynamica (via de Wick-rotatie $t \to i\tau$) voor systemen met reële, niet-negatieve spectra. Deze analogie suggereert dat kwantumsimulatoren potentieel ingezet kunnen worden om complexe klassieke stochastische problemen die gedomineerd worden door ruis en fluctuaties op te lossen. De auteurs merken echter bescheiden de beperkingen op: de huidige benadering vereist reële spectra en kan complexe eigenwaarden (oscillerend gedrag) niet direct adresseren; bovendien rust de analytische constructie op specifieke algebraïsche afsluitingen die niet voor alle interagerende problemen bestaan. De methode ruilt effectief de complexiteit van niet-monotone entropie in voor de complexiteit van het berekenen van statistische gemiddelden via de metrische operator $\tilde{\Omega}$.