Stochastic unravelings for Heisenberg picture and trace-nonpreserving dynamics

1. Problématique

Les désenchevêtrements stochastiques (stochastic unravelings) sont des outils puissants pour simuler la dynamique des systèmes quantiques ouverts. Traditionnellement, ces méthodes sont appliquées aux équations maîtresses (ME) dans le cadre de Schrödinger, où l'évolution de l'état du système est décrite par une application complètement positive et préservant la trace (TP). Dans ce cadre, la dynamique est interprétée comme une moyenne sur des trajectoires de vecteurs d'état purs.

Cependant, plusieurs limitations majeures existent :

• Cadre de Heisenberg : De nombreux problèmes physiques sont plus naturellement formulés dans le cadre de Heisenberg (évolution des opérateurs). Bien que des méthodes existent pour les systèmes gaussiens, il n'existait pas de méthode générale pour désenchevêtrer des équations maîtresses générales dans ce cadre.
• Non-préservation de la trace (TNP) : De nombreuses dynamiques physiques importantes ne préservent pas la trace. Cela inclut :
  • Les Hamiltoniens non hermitiens (points exceptionnels).
  • Les statistiques de comptage complet (full counting statistics) et les champs de comptage inclinés (tilted Lindbladians).
  • Les processus avec perte ou gain de particules.
    Les méthodes existantes (comme Monte-Carlo Wave Function - MCWF) échouent souvent ou nécessitent des hypothèses restrictives (comme la divisibilité CP stricte) pour traiter ces cas TNP.

L'objectif de ce travail est de développer un cadre général permettant d'étendre les désenchevêtrements stochastiques (notamment les processus déterministes par morceaux) au cadre de Heisenberg et à des équations maîtresses arbitraires non préservant la trace (TNP), tout en ne requérant que la positivité et l'hermiticité de la dynamique.

2. Méthodologie

Les auteurs proposent une généralisation des méthodes de désenchevêtrement existantes (MCWF et Rate Operator - RO) en introduisant deux mécanismes clés :

A. Gestion de la trace non préservée (Création et Disparition)

Dans les méthodes standards, la somme des probabilités de saut et d'évolution déterministe est égale à 1. Pour les équations TNP, cette somme peut être inférieure ou supérieure à 1. Les auteurs introduisent un processus stochastique supplémentaire pour compenser ce déséquilibre :

• Si la somme < 1 (Trace décroissante) : Une probabilité résiduelle correspond à la disparition (vanishing) de la réalisation stochastique. La trajectoire est supprimée de l'ensemble.
• Si la somme > 1 (Trace croissante) : Une probabilité résiduelle correspond à la création d'une copie (cloning) de la trajectoire actuelle. Les deux copies évoluent ensuite indépendamment.

Cela permet de maintenir la normalisation des vecteurs d'état tout en faisant varier le nombre total de réalisations $N(t)$ au cours du temps. La trace de l'opérateur moyen est alors donnée par le rapport $N(t)/N(0)$.

B. Extension au cadre de Heisenberg

Dans le cadre de Heisenberg, les opérateurs $X(t)$ ne sont pas nécessairement positifs ni hermitiens. Pour contourner cela :

1. Tout opérateur $X$ est décomposé en parties hermitienne et anti-hermitienne : $X = X_h + iX_a$.
2. Chaque partie est décomposée en différence de deux opérateurs positifs : $X_{h,a} = X^+_{h,a} - X^-_{h,a}$.
3. Chaque opérateur positif est proportionnel à une matrice densité. On applique alors le désenchevêtrement stochastique séparément à ces quatre composantes ($X^+_h, X^-_h, X^+_a, X^-_a$).
4. La solution finale est reconstruite comme une somme pondérée de ces trajectoires.

C. Adaptation des algorithmes (MCWF et RO)

• MCWF (Monte-Carlo Wave Function) : Les probabilités de saut et d'évolution sont recalculées en utilisant les générateurs de Heisenberg ($\tilde{H}, R_j, \Gamma_H$). Les probabilités de création/disparition sont dérivées de la différence entre l'opérateur qui rendrait l'équation TP et l'opérateur réel.
• RO (Rate Operator) : Cette méthode, qui utilise la décomposition spectrale d'un opérateur de taux, est adaptée de manière similaire. L'avantage du RO est qu'il peut gérer des dynamiques non divisibles (P-divisibles mais pas CP-divisibles) sans nécessiter de sauts inverses (reverse jumps) complexes, améliorant ainsi l'efficacité numérique.

D. Gestion des taux négatifs (Non-Markovianité)

Pour les régimes non-Markoviens où les taux peuvent devenir négatifs, le cadre intègre des sauts inverses (reverse jumps), permettant de corriger les probabilités négatives en réutilisant des trajectoires existantes, tout en conservant les mécanismes de création/disparition pour la trace.

3. Contributions Clés

1. Cadre Unifié TNP et Heisenberg : Première méthode générale permettant de désenchevêtrer des équations maîtresses TNP dans le cadre de Heisenberg, au-delà des systèmes gaussiens.
2. Efficacité Numérique via la Divisibilité : L'article démontre qu'il existe des dynamiques qui sont CP-divisibles dans le cadre de Heisenberg mais non divisibles dans le cadre de Schrödinger. En travaillant dans le cadre de Heisenberg, on peut utiliser des méthodes simples (sans sauts inverses) là où le cadre de Schrödinger exigerait des simulations lourdes avec dépendance mutuelle entre trajectoires.
3. Application aux Statistiques de Comptage : La méthode est appliquée avec succès aux "tilted Lindbladians" pour calculer les moments de la distribution de probabilité des photons comptés, une tâche difficile avec les méthodes directes en raison du bruit numérique lors du calcul des dérivées.
4. Flexibilité Algorithmique : La méthode est compatible avec différents schémas (MCWF, RO) et gère nativement les processus de création/disparition de trajectoires pour maintenir la précision de la trace.

4. Résultats et Validation

Les auteurs valident leur approche à travers plusieurs exemples numériques :

• Exemple 1 (Divisible en Heisenberg) : Une dynamique de spin où les taux sont positifs en Heisenberg mais négatifs en Schrödinger.
  • Résultat : La méthode MCWF en Heisenberg reproduit exactement la solution analytique sans nécessiter de sauts inverses. La trace de l'opérateur observé oscille correctement autour de sa valeur initiale, confirmant que le nombre de trajectoires varie conformément à la dynamique de la trace.
• Exemple 2 (Non-divisible en Heisenberg) : Une dynamique de phase covariante avec des taux négatifs (violation de la divisibilité CP).
  • Résultat : L'utilisation de la méthode Rate Operator (RO) permet de simuler la dynamique sans sauts inverses, même avec des taux négatifs. L'efficacité est accrue car seules quelques états propres sont nécessaires pour décrire la dynamique réduite. La trace augmente monotonement, comme prévu.
• Exemple 3 (Statistiques de comptage de photons) : Application aux équations pour les moments $\tau_k$ dérivés des statistiques de comptage.
  • Résultat : La méthode permet de calculer les premiers moments de la distribution de photons avec une précision élevée, évitant le bruit excessif lié à la dérivation numérique directe de la trace.

5. Signification et Impact

Ce travail élargit considérablement le domaine d'application des simulations stochastiques en physique quantique :

• Nouveaux Régimes Physiques : Il ouvre la voie à la simulation efficace de dynamiques non-hermitiennes, de systèmes avec perte/gain de particules et de statistiques de comptage complet, qui étaient auparavant difficiles à simuler avec des méthodes de trajectoires pures.
• Optimisation des Calculs : En exploitant les propriétés de divisibilité spécifiques au cadre de Heisenberg, les auteurs montrent que certaines simulations complexes peuvent être réalisées de manière beaucoup plus efficace (moins de ressources mémoire, pas de dépendance entre trajectoires) en changeant simplement de picture.
• Outil Versatile : La méthode est compatible avec les régimes non-Markoviens et s'adapte aux schémas de désenchevêtrement existants, offrant un outil robuste pour l'étude des systèmes quantiques ouverts généraux.

En résumé, cet article fournit une généralisation théorique et pratique essentielle qui rend les méthodes de trajectoires stochastiques applicables à une classe beaucoup plus large de problèmes physiques, en particulier ceux impliquant des opérateurs non-hermitiens et des dynamiques non conservatrices de la trace.