Semiclassical Langevin dynamics of long-range dissipative time crystals

Résumé technique : Dynamique de Langevin semi-classique des cristaux temporels dissipatifs à longue portée

Énoncé du problème
Les cristaux temporels de bordure (BTC) dissipatifs représentent une phase de la matière où la symétrie de translation temporelle continue ou discrète est spontanément brisée dans la limite thermodynamique, résultant en des oscillations collectives persistantes. Bien que les modèles idéalisés avec des interactions à portée infinie (totalement collectives) soutiennent un comportement de cycle limite stable décrit par des équations de champ moyen, les systèmes réels impliquent souvent des interactions à portée finie et une dissipation locale. Une question cruciale reste ouverte concernant la robustesse de la phase de cristal temporel lorsque la structure collective idéale est affaiblie. Plus précisément, il est nécessaire de comprendre comment les effets de taille finie influencent la durée de vie de ces oscillations et d'identifier les seuils de portée d'interaction (paramétrés par un exposant de loi de puissance $\alpha$) au-delà desquels la phase de cristal temporel se dissout dans un état stationnaire. Les simulations numériques exactes utilisant l'équation de maître de Lindblad sont limitées à de petites tailles de système, ce qui nécessite des méthodes alternatives pour sonder la mise à l'échelle de taille finie dans la limite thermodynamique.

Méthodologie
Les auteurs développent une approche de Langevin semi-classique ancrée dans les cadres de l'approximation de Wigner tronquée (TWA) et de la TWA discrète (DTWA). Au lieu de résoudre l'équation de maître de Lindblad complète, la dynamique est formulée via des équations de Langevin quantiques pour les opérateurs de spin locaux.

• Formalisme : La méthode remplace les produits d'opérateurs par des produits de leurs valeurs attendues (approximation de champ moyen) tout en conservant des termes de bruit stochastique pour rendre compte des fluctuations quantiques de premier ordre générées par les canaux dissipatifs.
• Modèle de spin-1/2 : Pour une chaîne de spin-1/2 avec des opérateurs de Lindblad à loi de puissance (la dissipation est à longue portée, le Hamiltonien est local), la dynamique est mappée sur des équations différentielles stochastiques pour les composantes de spin locales $s_i^\mu$. Le bruit est représenté par des bruits blancs gausiens réels.
• Modèle de spin-1 : Pour un système de spin-1 avec une dissipation locale et des interactions Hamiltoniennes à longue portée, la dynamique est formulée en utilisant une base de matrices de Gell-Mann (les génératrices de $SU(3)$). Les variables locales sont les valeurs attendues de ces matrices, $\lambda_i^\mu$.
• Implémentation numérique : Les équations stochastiques sont résolues en utilisant une interprétation d'Itô avec un schéma de Runge-Kutta pour la dérive déterministe et des incréments de Wiener gausiens pour le bruit. Crucialement, les auteurs constatent que pour les modèles étudiés, l'échantillonnage sur différentes conditions initiales (standard dans la DTWA) est inutile ; une seule condition initiale suffit car la mémoire de l'état initial est perdue rapidement par rapport à la période d'oscillation.
• Diagnostic : L'ordre du cristal temporel est diagnostiqué en analysant la mise à l'échelle de taille finie de la stabilité des oscillations. Pour le spin-1/2, l'accent est mis sur le taux de décroissance de l'enveloppe d'oscillation. Pour le spin-1, deux sondes sont utilisées : le « temps de déviation » (le temps pendant lequel la trajectoire de taille finie reste proche de la limite du cycle limite de champ moyen) et la mise à l'échelle de la puissance du pic de Fourier dominant.

Résultats clés

1. Modèle de spin-1/2 (Dissipation en loi de puissance) :

   • Dans la limite totalement collective ($\alpha=0$), le système présente des oscillations robustes où l'amortissement de taille finie est supprimé à mesure que la taille du système $L$ augmente.
   • Pour des exposants de loi de puissance $\alpha$ finis, l'enveloppe d'oscillation décroît exponentiellement en temps, mais le taux de décroissance $D$ suit une loi algébrique avec la taille du système : $D(L) \sim L^{-\beta}$.
   • Un exposant de mise à l'échelle $\beta$ positif indique que le taux de décroissance s'annule dans la limite thermodynamique, signifiant un cristal temporel stable.
   • L'exposant $\beta$ diminue de manière monotone à mesure que $\alpha$ augmente. Le régime de cristal temporel persiste pour $\alpha \lesssim 1,2$. Notamment, ce seuil est légèrement supérieur à $\alpha=1$. Alors que la théorie du champ moyen prédit une relaxation vers un point fixe pour $\alpha > 1$ (où le champ moyen n'est plus exact dans la limite thermodynamique), l'approche de Langevin semi-classique prédit des oscillations persistantes dans la plage $1 < \alpha \lesssim 1,2$. Cela suggère que les fluctuations quantiques, capturées par la méthode, stabilisent le cristal temporel au-delà du régime de champ moyen.

2. Modèle de spin-1 (Dissipation locale, Hamiltonien à longue portée) :

   • Le système présente un comportement de cristal temporel lorsque l'interaction cohérente est suffisamment étendue à longue portée.
   • Temps de déviation : Pour $\alpha \leq 1$, le temps $t^*$ auquel la trajectoire de taille finie dévie de la solution de champ moyen suit une loi de type $t^*(L) \sim L^b$. L'exposant $b$ est positif pour de faibles valeurs de $\alpha$ (environ 0,37 pour $\alpha=0$) et diminue à mesure que $\alpha$ approche 1, s'annulant près de $\alpha \approx 1,2$. Cette analyse n'est strictement valide que pour $\alpha \leq 1$, où la théorie du champ moyen décrit la limite thermodynamique.
   • Mise à l'échelle du pic de Fourier : La puissance du pic de Fourier dominant suit une loi $P(L) \sim L^s$. L'exposant $s$ est positif pour $\alpha \lesssim 1,2$, indiquant un comportement de cristal temporel. Contrairement à l'analyse du temps de déviation, cette sonde ne repose pas sur l'hypothèse que la limite thermodynamique est de type champ moyen, ce qui la rend robuste pour $\alpha > 1$.
   • Les deux sondes identifient un seuil critique pour la disparition du cristal temporel dans l'intervalle $\alpha \in [1, 1,5]$, ce qui est cohérent avec les résultats précédents de l'expansion des cumulants.

Signification et affirmations
L'article affirme que la dynamique de Langevin semi-classique fournit une sonde de taille finie utile pour les cristaux temporels dissipatifs dans les régimes où les simulations de Lindblad exactes sont numériquement intraitables.

• Au-delà du champ moyen : L'approche capture avec succès la robustesse de l'ordre de cristal temporel dans les régimes où la théorie du champ moyen échoue (spécifiquement pour $\alpha > 1$ dans le modèle de spin-1/2). Elle démontre que la phase de cristal temporel s'étend au-delà de la limite stricte du champ moyen, à condition que la dissipation ou l'interaction reste suffisamment étendue à longue portée.
• Outils de diagnostic : Ce travail établit les exposants de mise à l'échelle de taille finie (taux de décroissance $\beta$, exposant du temps de déviation $b$, et exposant du pic de Fourier $s$) comme des diagnostics directs de la robustesse de la phase de cristal temporel.
• Validation méthodologique : Les résultats pour le modèle de spin-1 s'alignent sur les découvertes précédentes issues des expansions de cumulants, validant l'approche de Langevin semi-classique comme un outil fiable pour l'étude des systèmes quantiques ouverts avec des interactions à longue portée.
• Intuition physique : L'étude souligne que bien que la limite thermodynamique pour $\alpha \leq 1$ soit décrite par la dynamique de champ moyen, la mise à l'échelle de taille finie vers cette limite est non triviale et dépend de $\alpha$. De plus, la mise à l'échelle algébrique du temps de déviation (plutôt que logarithmique) suggère une dynamique régulière et non chaotique dans ces phases de cristal temporel.