Exact diagonalization study of energy level statistics in harmonically confined interacting bosons

Cet article emploie la diagonalisation exacte pour démontrer que, tandis que les bosons en interaction confinés de manière harmonique présentent des statistiques régulières (Poisson) ou faiblement chaotiques dans des régimes d'interaction modérés, ils transitent vers un comportement chaotique fort caractérisé par des distributions GOE dans des régimes d'interaction forts, la rotation accentuant davantage ce chaos.

L'essentiel

Imaginez une piste de danse bondée à l'intérieur d'un bol géant et invisible. Sur cette piste, nous avons des groupes de danseurs identiques (des bosons) qui essaient tous de bouger sur le même rythme. Le bol a la forme d'une crêpe plate (un « plan quasi-2D »), et la musique est un bourdonnement rythmique et régulier (le « piège harmonique »).

Les danseurs ont deux éléments principaux qui influencent leurs mouvements :

1. La forme du bol : Les parois du bol les poussent vers le centre. C'est l'« énergie de piège ».
2. L'espace personnel des danseurs : Les danseurs n'aiment pas se cogner les uns aux autres. Il existe une force de répulsion, comme du papier bulle invisible, qui les repousse lorsqu'ils sont trop proches. C'est l'« énergie d'interaction ».

Les scientifiques de cet article voulaient savoir : quand ces danseurs bougent, leur motif est-il ordonné et prévisible, ou chaotique et aléatoire ?

Pour le découvrir, ils n'ont pas seulement observé la danse ; ils ont examiné les « niveaux d'énergie » (les pas ou les notes spécifiques que les danseurs peuvent effectuer). Ils ont utilisé une boîte à outils mathématique spéciale pour voir si les écarts entre ces pas étaient aléatoires ou s'ils suivaient une règle stricte.

Les deux scénarios principaux

Les chercheurs ont testé deux réglages d'« humeur » différents pour la piste de danse :

1. La danse « Calme » (Interaction modérée)

• L'installation : Les danseurs sont polis. La force qui les pousse à s'écarter est faible par rapport à la force du bol qui les maintient en place.
• Le résultat (Sans rotation) : Lorsque le bol ne tourne pas, les danseurs se déplacent de manière très ordonnée et prévisible. Leurs pas suivent une « distribution de Poisson ».
  • Analogie : Imaginez une file de personnes attendant un bus. Elles se tiennent à des intervalles aléatoires, mais elles ne se soucient pas les unes des autres. Parfois deux personnes sont proches, parfois loin. Il n'y a pas de « répulsion de niveau » (elles ne s'évitent pas activement). C'est un système régulier, non chaotique.
• Le résultat (En rotation) : Si vous faites tourner le bol lentement (créant un vortex unique), les danseurs deviennent un peu plus agités. Ils commencent à montrer des signes de « chaos faible ». Ils ne sont pas encore totalement aléatoires, mais ils ne sont plus parfaitement ordonnés non plus.

2. La danse « Sauvage » (Interaction forte)

• L'installation : Les danseurs sont très insistants. La force qui les pousse à s'écarter est maintenant aussi forte que les parois du bol.
• Le résultat (Sans rotation) : Soudain, la piste de danse devient chaotique. Les pas ne semblent plus aléatoires ; ils ressemblent à un système complexe et chaotique.
  • Analogie : Maintenant, les danseurs s'évitent activement. Si une personne fait un pas, les autres se déplacent immédiatement pour éviter de la heurter. C'est ce qu'on appelle la « répulsion de niveau ». Le motif des pas correspond désormais à la « distribution GOE » (Ensemble Orthogonal Gaussien), qui est l'empreinte mathématique du chaos.
• Le résultat (En rotation) : Lorsque vous faites tourner le bol pendant que les danseurs sont très insistants, le chaos s'emballe. Le système devient fortement chaotique.

Le rebondissement : Combien de danseurs ?

Les chercheurs ont également modifié le nombre de danseurs (12, 16 ou 20).

• Dans le scénario Calme, ajouter plus de danseurs rendait le système plus ordonné (plus proche de la file d'attente aléatoire du bus).
• Dans le scénario Sauvage, ajouter plus de danseurs faisait fluctuer le chaos. Parfois, il devenait plus chaotique, parfois il se calmait un peu, mais il restait généralement dans la zone chaotique.

Le facteur « Rotation »

L'article a découvert que la rotation est l'amplificateur de chaos ultime.

• Même lorsque les danseurs étaient seulement modérément insistants, faire tourner le bol les faisait agir de manière plus chaotique.
• Lorsque les danseurs étaient déjà très insistants, faire tourner le bol rendait le chaos encore plus fort.
• Ils ont même testé la rotation rapide du bol (créant 2 ou 3 vortex). Dans ces cas, le système était purement chaotique, quel que soit le nombre de danseurs sur la piste.

Les outils utilisés (Simplifiés)

Pour mesurer cela, les scientifiques ont utilisé quatre « règles » différentes :

1. Espacement entre voisins les plus proches (NNSD) : Mesurer la distance entre un pas et le suivant.
2. Rapport des espacements : Comparer la distance entre le pas A et B, à la distance entre B et C. (C'est une astuce ingénieuse pour éviter certaines erreurs mathématiques).
3. Règles à longue portée (Dyson-Mehta & Variance du nombre de niveaux) : Ces outils regardaient le motif sur une longue série de pas pour voir si toute la piste de danse était rigide ou flexible.

L'essentiel à retenir

L'article conclut que le comportement de ces atomes piégés est un bras de fer entre le piège (qui veut de l'ordre) et l'interaction (qui crée de la complexité).

• Interaction faible + Pas de rotation = Ordonné (Régulier).
• Interaction forte OU Rotation = Chaotique.
• Interaction forte + Rotation = Chaos maximum.

Essentiellement, l'étude montre qu'en changeant simplement la force avec laquelle les particules se poussent ou la vitesse à laquelle le système tourne, on peut faire passer un système quantique d'une machine d'horlogerie prévisible à une tempête chaotique et sauvage. Cela aide les scientifiques à comprendre comment le « chaos quantique » émerge dans le monde réel, spécifiquement dans les gaz ultra-froids comme ceux faits d'atomes de Rubidium.

Résumé technique

Résumé Technique : Étude par Diagonalisation Exacte de la Statistique des Niveaux d'Énergie de Bosons en Confinement Harmonique

Énoncé du Problème
L'article étudie les propriétés statistiques des spectres d'énergie dans les systèmes à plusieurs corps quantiques, en se concentrant spécifiquement sur les bosons interagissants confinés harmoniquement dans un plan quasi-2D. Bien que la Théorie des Matrices Aléatoires (RMT) fournisse un cadre pour distinguer les systèmes réguliers (intégrables) des systèmes chaotiques (non intégrables) via les conjectures de Bohigas-Giannoni-Schmit (BGS) et de Berry-Tabor, les effets spécifiques de la variation de la force d'interaction à deux corps, du nombre de particules ($N$) et de la rotation sur les bosons piégés et interagissants n'avaient pas été pleinement explorés. Les études précédentes utilisaient souvent des méthodes approximatives ou se concentraient sur des limites spécifiques. Ce travail vise à combler cette lacune en effectuant une étude de diagonalisation exacte pour comprendre comment l'interrelation entre l'énergie d'interaction et l'énergie de piégeage, combinée à la rotation, induit la transition d'un comportement spectral régulier vers un comportement chaotique.

Méthodologie
Les auteurs emploient une approche de diagonalisation exacte pour résoudre le hamiltonien à plusieurs corps pour $N = 12, 16$ et $20$ bosons de $^{87}\text{Rb}$.

• Modèle de Système : Les bosons sont confinés dans un piège harmonique quasi-2D avec un paramètre d'anisotropie $\lambda_z = \omega_z/\omega_\perp = 4$. L'interaction est modélisée par un potentiel gaussien répulsif avec une portée $\sigma = 0,1 a_\perp$. Le système est analysé à la fois dans des référentiels non rotatifs ($L_z = 0$) et rotatifs (état de vortex unique $L_z = N$, et moments angulaires plus élevés $L_z = 2N, 3N$).
• Régimes d'Interaction : Deux régimes distincts sont examinés :
  1. Interaction Modérée : L'énergie d'interaction est faible par rapport à l'énergie de piégeage ($g_2 = 0,3669$, correspondant à $a_{sc} = 1000 a_0$).
  2. Interaction Forte : L'énergie d'interaction est comparable à l'énergie de piégeage ($g_2 = 3,669$, correspondant à $a_{sc} = 10000 a_0$).
• Outils Statistiques : L'étude analyse les 100 premiers niveaux d'énergie en utilisant des mesures de corrélation à courte et longue portée :
  • Courte Portée : La distribution des espacements entre voisins les plus proches (NNSD) $P(s)$ et la distribution du rapport des espacements de niveaux consécutifs $P(r)$. Les auteurs soulignent l'importance de $P(r)$ car elle évite la procédure de dépliage (unfolding) requise pour $P(s)$, minimisant ainsi les erreurs systématiques. La distribution de Brody est utilisée pour ajuster les données, avec un paramètre $b$ qui interpole entre Poisson ($b=0$) et l'Ensemble Orthogonal de Wigner (GOE, $b=1$).
  • Longue Portée : La statistique de Dyson-Mehta $\Delta_3(L)$ et la variance du nombre de niveaux $\Sigma^2(L)$.

Contributions Clés et Résultats

1. Cas Non Rotatif ($L_z = 0$) :

   • Interaction Modérée : Le système présente un comportement régulier. Les distributions NNSD et $P(r)$ s'alignent sur la distribution de Poisson, indiquant une absence de répulsion des niveaux. Le paramètre de Brody $b$ diminue systématiquement avec l'augmentation de $N$ (passant de $0,33$ pour $N=12$ à $0,04$ pour $N=20$), suggérant que les corrélations spectrales diminuent à mesure que le nombre de bosons augmente dans ce régime. Les statistiques à longue portée ($\Delta_3$ et $\Sigma^2$) suivent également un comportement de Poisson pour de petits intervalles d'énergie.
   • Interaction Forte : Le système transite vers un comportement chaotique. Les distributions s'alignent sur le GOE, signifiant une forte répulsion des niveaux. Cependant, le degré de chaos est modulé par $N$. Pour $N=16$, le système présente les signatures de chaos les plus fortes ($b=0,85$, $\langle r \rangle \approx 0,538$), tandis que pour $N=12$ et $N=20$, le système revient à un régime faiblement chaotique avec une certaine régularité. La statistique $\Delta_3(L)$ sature à des valeurs significativement plus basses dans le régime fort par rapport au régime modéré, confirmant la transition vers un système non intégrable.

2. Cas Rotatif (État de Vortex Unique $L_z = N$) :

   • Interaction Modérée : La rotation induit une transition d'un comportement régulier vers un comportement faiblement chaotique. Le système présente des signatures de chaos avec un certain degré de régularité (paramètres de Brody $b \approx 0,52 - 0,64$). Le rapport de gap moyen $\langle r \rangle$ est proche de $0,5$, intermédiaire entre Poisson et GOE.
   • Interaction Forte : Le système affiche un comportement fortement chaotique cohérent avec les statistiques GOE ($b \approx 1,0$, $\langle r \rangle \approx 0,55$). Il a été constaté que la rotation renforce la nature chaotique du système par rapport au cas non rotatif dans le même régime d'interaction.

3. Moments Angulaires Supérieurs ($L_z = 2N, 3N$) :

   • Dans le régime d'interaction forte, les systèmes avec $L_z = 2N$ et $L_z = 3N$ présentent de fortes signatures de chaos quantique, avec les statistiques à courte portée ($P(r)$) et à longue portée ($\Delta_3$) s'alignant sur la distribution GOE.

Signification et Revendications
Les auteurs affirment que leur étude fournit une compréhension complète de la manière dont la force d'interaction, le nombre de particules et la rotation gouvernent les corrélations spectrales dans les bosons piégés et interagissants.

• Mécanisme de Transition : L'interrelation entre l'énergie d'interaction et l'énergie de piégeage régit la transition du comportement régulier (Poisson) vers le comportement chaotique (GOE).
• Rôle de la Rotation : La rotation est identifiée comme un facteur qui renforce le comportement chaotique dans les régimes d'interaction modérée et forte.
• Robustesse Méthodologique : L'étude valide l'utilisation du rapport des espacements de niveaux consécutifs $P(r)$ comme un outil fiable pour l'analyse des corrélations à courte portée, car elle évite les erreurs potentielles associées à la procédure de dépliage requise pour la NNSD.
• Universalité : Les résultats soutiennent l'applicabilité universelle de la Théorie des Matrices Aléatoires pour décrire les corrélations spectrales dans les systèmes quantiques à plusieurs corps, tels que les systèmes de Bose ultra-froids.

L'article conclut modestement, suggérant que ces découvertes font progresser la compréhension du chaos quantique dans les systèmes de Bose ultra-froids et notant que les travaux futurs pourraient étendre cette analyse aux facteurs de forme spectraux et aux modèles de Bose-Hubbard étendus avec des interactions à longue portée.