Quantum instanton approach to metastable collective spins

1. Problématique et Contexte

Les systèmes de spins collectifs, composés d'ensembles de spins couplés à un réservoir commun (comme une cavité optique amortie), sont fondamentaux en physique atomique et en physique de la matière condensée. Ces systèmes peuvent être décrits efficacement par un seul « macrospin ».

• Multistabilité et Transitions de Phase : En raison de leur non-linéarité intrinsèque, ces systèmes peuvent présenter plusieurs états métastables à longue durée de vie. Selon les paramètres de contrôle, l'état stationnaire dominant peut changer brusquement, signalant une transition de phase dissipative du premier ordre.
• Le Défi du Calcul : Pour un grand nombre de spins ($J \to \infty$), l'analyse via l'équation maîtresse quantique (QME) devient numériquement impossible. Les approches classiques utilisent une échelle asymptotique de type Arrhenius pour les taux de transition entre états métastables : $\kappa_{i \to j} \sim \exp(-J A_{i \to j})$, où $A_{i \to j}$ est la barrière d'activation.
• Limites des Approches Semiclassiques : Les études précédentes sur les systèmes de spins reposaient souvent sur l'approximation de Wigner (SW), qui utilise une équation de Fokker-Planck tronquée (négligeant les dérivées d'ordre supérieur). L'article montre que cette approximation échoue à capturer correctement les barrières d'activation car elle ignore les fluctuations non gaussiennes, conduisant à des prédictions erronées sur la position des transitions de phase et les temps de relaxation.

2. Méthodologie : L'Approche Instanton Quantique

Les auteurs développent une nouvelle méthode basée sur la dynamique des quasi-probabilités quantiques exactes, sans troncature, pour calculer les barrières d'activation.

• Représentation des Quasi-probabilités : Au lieu de la fonction de Wigner, ils utilisent les représentations de Husimi ($p_H$) et de P ($p_P$) pour la matrice densité. Ces distributions satisfont des équations d'évolution différentielles exactes de la forme $\partial_t p_\alpha = \mathcal{L}_\alpha p_\alpha$.
• Limite de Grand Spin et Ansatz WKB : Dans la limite $J \to \infty$, le propagateur de ces équations est analysé via un ansatz WKB : $K_\alpha \sim \exp(-J S_\alpha)$. L'action $S_\alpha$ obéit à une équation de Hamilton-Jacobi avec un hamiltonien auxiliaire $H_\alpha(x, \pi)$, où $\pi$ joue le rôle de moment conjugué.
• Détermination des Instantons :
  • Les trajectoires instantons (chemins de transition les plus probables) sont des solutions des équations de Hamilton associées à $H_\alpha$.
  • Contrairement aux systèmes stochastiques classiques où le hamiltonien est convexe, ici $H_\alpha$ est non convexe en $\pi$. Cela permet l'existence de solutions avec une action négative (physiquement non pertinentes) et nécessite des critères de sélection rigoureux.
  • Les auteurs imposent des conditions aux limites spécifiques : le moment $\pi$ doit être nul aux points de départ et d'arrivée (les bassins d'attraction stables), et l'action doit converger vers une valeur positive constante pour $t \to \infty$.
  • Grâce à la symétrie du modèle, le problème est réduit à un plan bidimensionnel ($w, \pi_w$), permettant un calcul efficace par méthode de continuation numérique.

3. Résultats Clés

L'approche a été testée sur un modèle de spin collectif avec un forçage cohérent, un pompage linéaire et une dissipation non linéaire.

• Précision des Barrières d'Activation :
  • La méthode instanton quantique prédit correctement les barrières d'activation $A_{i \to j}$.
  • Le point de croisement des barrières (où $A_{\ell \to u} = A_{u \to \ell}$), qui marque la transition de phase, est trouvé à $\Gamma \approx 8.9\gamma$.
  • En comparaison, l'approche semiclassique de Wigner (SW) prédit un croisement à $\Gamma \approx 7.9\gamma$, ce qui est en désaccord avec les simulations QME exactes pour des $J$ finis. L'approche SW sous-estime les barrières car elle tronque les termes d'ordre supérieur (fluctuations non gaussiennes).
• Validation par le Gap de Liouvillien :
  • Le taux de relaxation le plus lent du système est caractérisé par le gap de Liouvillien $\lambda$.
  • Les auteurs montrent que $\lambda \sim \exp(-J A_{\min})$.
  • La décroissance exponentielle de $\lambda$ calculée via la QME correspond parfaitement aux prédictions de la méthode instanton quantique, mais est significativement sous-estimée par l'approche SW.
• Résolution du Paradoxe de Keizer : La méthode explique comment, pour un $J$ fini mais grand, les attracteurs de la dynamique de champ moyen (MF) se transforment en états métastables avec des durées de vie finies, et comment le système bascule d'un état à l'autre en fonction de la taille du système.

4. Contributions et Signification

• Preuve de Concept : C'est la première démonstration explicite d'une approche instanton basée sur les équations de mouvement exactes des quasi-probabilités pour les systèmes de spins, dépassant les limites de l'approximation semiclassique.
• Supériorité sur Wigner : L'article démontre de manière convaincante que l'approximation de Wigner, souvent utilisée par défaut, est insuffisante pour décrire la métastabilité dans les systèmes de spins quantiques non linéaires en raison de la négligence des fluctuations non gaussiennes.
• Généralité : Bien que appliquée aux spins collectifs, la méthode est générale pour tout système admettant une dynamique de quasi-probabilité de la forme étudiée. Elle pourrait être étendue aux systèmes bosoniques (offrant une alternative conceptuellement plus simple aux intégrales de chemin de Keldysh), aux complexes spin-boson, et aux systèmes avec dissipation locale ou contrôlée par rétroaction.
• Impact : Cette approche fournit un outil théorique robuste pour prédire les temps de relaxation et les points de transition de phase dans les systèmes quantiques ouverts complexes, cruciaux pour le développement de technologies quantiques comme les qubits de type "chat de Schrödinger" ou les simulateurs quantiques.

En résumé, cet article établit un cadre théorique rigoureux pour analyser la métastabilité quantique au-delà des approximations semiclassiques, en corrigeant les erreurs systématiques des méthodes précédentes et en offrant une description précise de la dynamique à long terme des systèmes de spins collectifs.