Circular Rosenzweig-Porter random matrix ensemble

Cet article propose et valide numériquement un analogue unitaire (circulaire) de l'ensemble de matrices aléatoires de Rosenzweig-Porter, défini via le mouvement brownien de Dyson, pour modéliser les statistiques de niveaux et la fractalité des états propres de la localisation à plusieurs corps dans les systèmes périodiquement entraînés (Floquet).

L'essentiel

Imaginez une immense piste de danse chaotique remplie de milliers de danseurs. Dans le monde de la physique quantique, ces danseurs représentent les états possibles d'un système complexe (comme un groupe d'atomes en interaction).

Ce papier traite de la compréhension du mouvement et des interactions de ces danseurs lorsque la musique change de deux manières différentes :

1. Systèmes statiques : La musique est un bourdonnement stable et inchangeant (comme une pièce normale et calme).
2. Systèmes de Floquet (pilotés périodiquement) : La musique est un rythme répétitif qui modifie les règles toutes les quelques secondes (comme un stroboscope ou un laser pulsé).

Pendant longtemps, les physiciens disposaient d'un excellent « livre de règles » (appelé l'ensemble de Rosenzweig-Porter) pour le premier scénario (la pièce statique). Ce livre de règles les aide à prédire si les danseurs se mélangeront librement (chaos) ou resteront coincés dans leurs propres petits coins (localisation).

Cependant, personne n'avait de bon livre de règles pour le deuxième scénario (la pièce pulsée et rythmée). Puisque les règles de la mécanique quantique changent lorsque les systèmes sont pilotés par un rythme, les anciennes mathématiques ne convenaient pas tout à fait.

La nouvelle idée : une piste de danse circulaire

Les auteurs de ce papier se sont demandé : « Peut-on construire une version de ce vieux livre de règles qui fonctionne pour les systèmes rythmiques et pulsés ? »

Ils ont créé un nouveau modèle qu'ils appellent l'ensemble circulaire de Rosenzweig-Porter.

Voici comment ils l'ont construit, en utilisant une analogie simple :

• L'ancienne méthode (Mouvement brownien) : Imaginez les danseurs se déplaçant au hasard sur une ligne droite et plate. Si vous les poussez au hasard au fil du temps, ils se dispersent de manière prévisible. C'est ainsi que fonctionnait l'ancien modèle.
• La nouvelle méthode (Mouvement circulaire) : Pour les systèmes rythmiques, les auteurs ont réalisé que les danseurs ne se déplacent pas sur une ligne droite ; ils se déplacent sur un cercle. Imaginez les danseurs courant sur une piste circulaire. Leurs positions sont mesurées par des angles (comme un cadran d'horloge) plutôt que par des distances linéaires.

Ils ont défini leur nouveau modèle comme le résultat d'une « marche aléatoire » qui se produit spécifiquement sur ce cercle. Ils n'ont pas simplement deviné les mathématiques ; ils ont simulé ce processus sur un ordinateur pour voir ce qui se passait.

Ce qu'ils ont découvert

Les auteurs ont effectué d'énormes simulations informatiques (impliquant jusqu'à 1 000 danseurs) pour voir si leur nouveau modèle « circulaire » se comportait comme l'ancien modèle « en ligne droite ». Ils ont vérifié deux points principaux :

1. L'espacement des danseurs (Niveaux d'énergie)
Ils ont examiné les écarts entre les danseurs.

• Dans la zone « Chaotique » : Les danseurs sont répartis uniformément, et les écarts entre eux suivent un motif spécifique et complexe (comme une foule où tout le monde se bouscule).
• Dans la zone « Localisée » : Les danseurs se regroupent ou restent très éloignés les uns des autres de manière très prévisible et simple (comme des personnes alignées).
• Le résultat : Leur nouveau modèle circulaire a montré exactement le même passage du « chaos » au « regroupement » que l'ancien modèle. Le « point de bascule » où le comportement change se situait au même endroit.

2. La forme des danseurs (États propres)
Ils ont examiné à quel point l'influence d'un seul danseur était « étendue ».

• Étendu : L'énergie d'un danseur est partagée parmi beaucoup d'autres.
• Fractale : Un danseur se trouve dans un état intermédiaire étrange — étendu, mais pas complètement. C'est comme un nuage ayant une forme floue et auto-similaire.
• Localisé : Un danseur est coincé à un endroit précis.
• Le résultat : Le modèle circulaire a reproduit exactement ces mêmes formes. Que les danseurs soient complètement mélangés, partiellement mélangés (fractals) ou coincés, le nouveau modèle correspondait parfaitement à l'ancien.

La conclusion

L'article affirme qu'ils ont réussi à construire une version unitaire (circulaire) du célèbre modèle de Rosenzweig-Porter.

En traitant le système comme un cercle plutôt que comme une ligne droite, ils ont créé un outil qui décrit avec précision le comportement des systèmes quantiques périodiquement pilotés (pulsés). Tout comme l'ancien modèle était un outil « phénoménologique » (descriptif) pour les systèmes statiques, ce nouveau modèle circulaire sert d'outil descriptif pour les systèmes qui sont secoués ou pilotés rythmiquement.

Ils l'ont prouvé en montrant que les « empreintes digitales » statistiques de leur nouveau modèle (la façon dont les niveaux sont espacés et dont les états sont façonnés) sont indiscernables de celles du modèle original, bien compris. Cela offre aux physiciens une nouvelle méthode fiable pour étudier les systèmes quantiques complexes et rythmiques sans avoir à résoudre des équations extrêmement difficiles à partir de zéro.

Résumé technique

Résumé technique : Ensemble de matrices aléatoires Rosenzweig-Porter circulaire

Énoncé du problème
L'ensemble de matrices aléatoires Rosenzweig-Porter (RP) est un modèle phénoménologique bien établi utilisé pour décrire les statistiques des niveaux et la fractalité des états propres à travers la transition de localisation à plusieurs corps (MBL) dans les systèmes quantiques interactifs désordonnés et statiques. Dans les systèmes statiques, la dynamique est régie par des Hamiltoniens hermitiens. Cependant, les systèmes périodiquement pilotés (Floquet) sont décrits par des opérateurs de Floquet unitaires, dont les valeurs propres se situent sur le cercle unité dans le plan complexe. Bien que la MBL ait été étendue aux systèmes de Floquet, un analogue unitaire (circulaire) direct de l'ensemble RP, capable de capturer la même phénoménologie pour les systèmes périodiquement pilotés, n'avait pas été construit. Cet article aborde la question de savoir si un tel analogue unitaire peut être défini et s'il partage les propriétés statistiques clés de l'ensemble RP original.

Méthodologie
Les auteurs proposent une construction de l'ensemble Rosenzweig-Porter circulaire (CRP) basée sur l'observation que l'ensemble RP standard peut être considéré comme le résultat d'un mouvement brownien de Dyson (DBM) à temps fini.

1. Construction : Les auteurs définissent l'ensemble CRP comme le résultat d'une évolution stochastique d'une matrice unitaire symétrique dépendante du temps $S(t)$ de dimension $N$. L'évolution suit un processus de mouvement brownien de Dyson adapté aux matrices unitaires :
   $$dS(t) = i\sqrt{dt} U^T(t) X U(t)$$
   où $X$ est une matrice de l'Ensemble Orthogonal Gaussien (GOE) échantillonnée à chaque étape de temps infinitésimale, et $U(t)$ est dérivée de la décomposition $S(t) = U^T(t)U(t)$.
2. Conditions initiales : Le processus commence avec une matrice diagonale $S(0) = \text{diag}(e^{i\theta_1(0)}, \dots, e^{i\theta_N(0)})$, où les phases initiales $\theta_n(0)$ sont échantillonnées indépendamment selon une distribution uniforme sur $[-\pi, \pi)$. Cela assure une densité d'états uniforme, cohérente avec les opérateurs de Floquet génériques.
3. Définition de l'ensemble : L'ensemble CRP est défini comme l'état de la matrice $S(t)$ à un temps spécifique $t = \epsilon^2 / N^\gamma$, où $\epsilon$ est une force de perturbation et $\gamma$ contrôle la force relative des termes diagonaux et hors diagonale.
4. Implémentation numérique : Les auteurs évaluent numériquement le processus stochastique en utilisant un algorithme qui met à jour la matrice via $S(t+dt) = S(t)e^{i\sqrt{dt}X}$. Cette approche préserve l'unitarité pour des pas de temps finis et évite les décompositions explicites de matrices. Les simulations ont été réalisées pour des dimensions de matrice $N \leq 1000$, avec une moyenne sur au moins $10^6$ valeurs propres ou états propres.

Contributions clés

• Définition du CRP : L'article introduit l'ensemble Rosenzweig-Porter circulaire comme un analogue unitaire de l'ensemble RP standard, défini explicitement par un processus de mouvement brownien de Dyson circulaire.
• Correspondance théorique : Les auteurs fournissent un argument heuristique selon lequel les statistiques spectrales de l'ensemble CRP au milieu du spectre devraient correspondre à celles de l'ensemble RP, à condition que la densité d'états soit mise à l'échelle de manière appropriée (spécifiquement, en mapant le paramètre RP $\epsilon$ vers $\epsilon/\sqrt{2\pi}$ dans le contexte CRP).
• Vérification numérique : L'étude fournit des preuves numériques complètes caractérisant les propriétés statistiques des valeurs propres et des états propres de l'ensemble CRP.

Résultats
L'enquête numérique confirme que l'ensemble CRP présente des comportements statistiques analogues à ceux de l'ensemble RP à travers différents régimes du paramètre $\gamma$ :

1. Statistiques des espacements de niveaux : Le rapport moyen des espacements de niveaux consécutifs, $r$, sert de diagnostic pour le chaos quantique.
   • Pour $\gamma < 2$, le système présente des statistiques de Wigner-Dyson ($r \approx 0,530$), caractéristiques des systèmes chaotiques.
   • Pour $\gamma \geq 2$, le système transitionne vers des statistiques de Poisson ($r \approx 0,386$), caractéristiques des systèmes intégrables.
   • Une transition de phase se produit au point critique $\gamma_c = 2$ avec un exposant critique $\nu = 1$. Les données montrent un effondrement d'échelle de $r$ en fonction de $(\gamma - \gamma_c) \ln(N)^{1/\nu}$, cohérent avec l'ensemble RP.
2. Fractalité des états propres : Le rapport d'inverse de participation (IPR) des états propres a été analysé dans la base où $S(0)$ est diagonale.
   • Pour $\gamma \leq 1$, les états propres sont étendus ($IPR_2 \sim N^{-1}$).
   • Pour $\gamma \geq 2$, les états propres sont localisés ($IPR_2 \sim O(1)$).
   • Dans le régime intermédiaire $1 < \gamma < 2$, les états propres sont fractals, évoluant selon $IPR_2 \sim N^{-(2-\gamma)}$.
3. Forme des états propres : La distribution des amplitudes des états propres $|\langle m | \psi_n \rangle|^2$ s'est avérée suivre des statistiques de Breit-Wigner. Les données s'alignent avec la prédiction théorique :
   $$|\langle m | \psi_n \rangle|^2 \sim \frac{1}{(\theta_n - \theta_m^{(0)})^2 + \Gamma^2}$$
   où $\Gamma$ est la largeur de diffusion. Bien que les effets de taille finie et les approximations de continuum entraînent de légères discordances quantitatives, l'accord qualitatif et les différences structurelles distinctes entre les régimes $\gamma = 0,5$ et $\gamma = 1,5$ sont clairement observés.

Signification et affirmations
Les auteurs affirment que l'ensemble Rosenzweig-Porter circulaire sert de modèle phénoménologique valide pour les statistiques des niveaux et la fractalité des états propres observées à travers la transition de localisation à plusieurs corps dans les systèmes périodiquement pilotés (Floquet). En établissant que l'ensemble CRP partage des propriétés statistiques clés avec l'ensemble RP, ce travail fournit un outil théorique maniable pour étudier les phases non ergodiques dans les systèmes quantiques pilotés.

L'article présente modestement ce travail comme une « première étape » vers la construction de généralisations plus vastes, telles que des modèles incorporant la multifractalité ou d'autres modèles de Floquet avec des états propres multifractals. Les auteurs suggèrent que de futures généralisations pourraient être obtenues en considérant des processus stochastiques avec des incréments corrélés, similaires aux propositions récentes pour l'ensemble RP standard. De plus, ils notent l'utilité potentielle de l'ensemble CRP comme bloc de construction non maximiquement aléatoire dans les études de circuits quantiques aléatoires.