Restoring detailed balance in non-Hermitian Markov processes

Résumé Technique : Restauration du Bilan Détaillé dans les Processus de Markov Non Hermitiens

Énoncé du Problème
Les processus stochastiques hors équilibre, particulièrement ceux modélisés par des processus de Markov à états discrets, impliquent fréquemment des taux de transition asymétriques qui brisent le bilan détaillé. Ces générateurs stochastiques non hermitiens ($H$) entraînent une production d'entropie non monotone, où l'entropie peut diminuer après des augmentations transitoires. Ce comportement compromet les interprétations thermodynamiques standards et les techniques d'inférence qui reposent sur la maximisation de l'entropie ou sur l'hypothèse d'une incertitude toujours croissante. De plus, l'absence de réversibilité temporelle dans ces systèmes complique l'interprétation physique des états d'équilibre, qui peuvent encoder des configurations de faible entropie. Bien que les méthodes de réseaux de tenseurs permettent de traiter la « malédiction de la dimensionnalité » dans les systèmes interactifs de haute dimension, elles ne résolvent pas le problème fondamental de la dynamique non hermitienne et de la production d'entropie non monotone qui en résulte.

Méthodologie
Les auteurs proposent l'utilisation de cartes de Dyson ($\eta$) pour transformer les générateurs stochastiques non hermitiens en opérateurs hermitiens tout en préservant le spectre des valeurs propres réelles. La méthodologie centrale implique :

1. Transformation de Dyson : Une transformation réversible est appliquée à l'équation de maîtrise $\partial_t |P(t)\rangle = -H|P(t)\rangle$, produisant une équation transformée $\partial_t |\phi(t)\rangle = -\tilde{H}|\phi(t)\rangle$, où $|\phi(t)\rangle = \eta|P(t)\rangle$ et $\tilde{H} = \eta H \eta^{-1}$.
2. Opérateur Métrique : La transformation introduit un opérateur métrique $\Omega = \eta^\dagger \eta$, qui repondère les contributions dans le système d'origine. Cet opérateur relie les mesures dans les coordonnées transformées au vecteur de probabilité d'origine, garantissant que les moyennes statistiques $\langle O \rangle$ peuvent être récupérées via $\langle \phi | \tilde{\Omega} \Xi' O' | \phi \rangle$.
3. Restauration de l'Entropie : Dans le système transformé, l'entropie de Rényi $S' = -\ln |\phi(t)|^2$ est prouvée comme augmentant de manière monotone au cours du temps ($\frac{1}{2}\frac{d}{dt}S'(t) \geq 0$). Cela contraste avec le système d'origine, où la norme au carré $|P(t)|^2$ (et donc l'entropie) peut présenter des minima locaux.
4. Construction Computationnelle :
   • Approche Analytique : Les auteurs cherchent un générateur hermitien $\Lambda = \ln \eta$ au sein d'une algèbre d'opérateurs fermée. En développant la condition $H = H^\dagger$ en termes de commutateurs, ils dérivent un ensemble d'équations algébriques (impliquant des coefficients $G_{ijk}$) pour résoudre $\Lambda$. Cette approche suppose une structure fermée sous le commutateur, ce qui peut ne pas être le cas pour tous les problèmes interactifs non triviaux.
   • Approche Numérique : Pour les cas généraux où la clôture analytique est difficile, un algorithme numérique itératif est proposé. Cette méthode traite $\Lambda$ comme un axe de rotation dans l'algèbre, effectuant des rotations itératives pour supprimer la composante anti-hermitienne ($\Delta H$) du générateur jusqu'à ce qu'une tolérance soit atteinte.

Résultats Clés
La méthodologie a été appliquée à un modèle épidémique Susceptible-Infecté-Susceptible (SIS) sur un graphe de contact statique :

• Croissance Monotone de l'Entropie : Dans le système transformé, l'entropie de Rényi $S'$ augmente de manière monotone jusqu'à l'atteinte de l'équilibre, alors que l'entropie du système d'origine présente un maximum transitoire suivi d'une diminution.
• Équivalence Statistique : Malgré la transformation, les moyennes statistiques (moyenne et écart-type des sommets infectés) calculées dans le système transformé correspondent aux résultats du système d'origine avec une grande précision (distance de Kolmogorov-Smirnov $< 10^{-3}$), à condition que l'opérateur métrique soit correctement appliqué.
• Détection de Transition de Phase : Le maximum de l'entropie de Rényi dans les deux systèmes (l'original et le transformé) coïncide au même ratio critique entre les taux d'infection et de guérison ($\beta/\gamma$). Cela confirme que le système transformé, bien qu'ayant restauré la réversibilité temporelle, conserve la signature de la transition de phase trouvée dans le système non hermitien d'origine.
• Préservation de la Structure : Les distributions à l'état stationnaire dans le système transformé ($|\phi\rangle$) préservent les indices structurels de la distribution d'origine ($|P\rangle$), suggérant que les contributions diagonales de la carte de Dyson régissent la transformation tandis que les composantes hors-diagonales servent à mélanger les composantes pour augmenter l'entropie.

Signification et Revendications
L'article affirme fournir une méthode générale et calculable pour restaurer le bilan détaillé dans une large classe de systèmes markoviens. En convertissant les générateurs non hermitiens en générateurs hermitiens, l'approche dote ces systèmes d'une « interprétation physique familière » caractérisée par une production d'entropie strictement positive. Cela permet l'application des outils standards de la physique statistique, tels que la maximisation de l'entropie, pour déterminer les états d'équilibre.

Les auteurs soulignent que ce formalisme établit une connexion entre les processus stochastiques et la dynamique quantique (via la rotation de Wick $t \to i\tau$) pour les systèmes possédant des spectres réels et non négatifs. Cette analogie suggère que des simulateurs quantiques pourraient potentiellement être utilisés pour résoudre des problèmes stochastiques classiques difficiles dominés par le bruit et les fluctuations. Cependant, les auteurs notent modestement certaines limites : l'approche actuelle nécessite des spectres réels et ne peut pas traiter directement les systèmes avec des valeurs propres complexes (comportement oscillatoire), et la construction analytique repose sur des clôtures algébriques spécifiques qui peuvent ne pas exister pour tous les problèmes interactifs. La méthode échange efficacement la complexité de l'entropie non monotone contre la complexité du calcul des moyennes statistiques via l'opérateur métrique $\tilde{\Omega}$.