Emergence of caustics in dynamics of the Kitaev model

Résumé technique : Émergence de caustiques dans la dynamique du modèle de Kitaev

Énoncé du problème
L'article étudie la dynamique hors équilibre des systèmes quantiques à plusieurs corps fermés, en se concentrant spécifiquement sur la propagation des quasi-particules dans les systèmes intégrables bidimensionnels (2D). Alors que les structures en cône de lumière, délimitées par la limite de Lieb-Robinson, sont bien établies dans les systèmes unidimensionnels (1D) et servent de signatures de caustiques quantiques (analogues aux singularités en optique géométrique), leur comportement dans les modèles intégrables 2D reste moins exploré. Les auteurs visent à identifier ces structures en cône de lumière et ces caustiques quantiques dans le modèle de Kitaev en nid d'abeilles 2D, à analyser leur anisotropie spatiale, et à déterminer comment des excitations périodiques externes (ingénierie de Floquet) modifient ces dynamiques.

Méthodologie
L'étude emploie une combinaison de cadres analytiques et de simulations numériques :

1. Cartographie du modèle : Le modèle de Kitaev en nid d'abeilles de spin-1/2 est mappé sur un modèle de fermions non interactifs en utilisant la transformation de Jordan-Wigner. Le système est analysé dans le secteur où la constante du mouvement $\hat{D}_n = 1$.
2. Dynamique des quasi-particules : La dynamique d'une seule quasi-particule initiée à un site spécifique est décrite à l'aide d'une fonction d'onde $\psi(r, t)$ exprimée comme une somme sur les modes d'impulsion dans la zone de Brillouin. La phase de cette fonction d'onde est définie par un fonctionnel d'action $S(k; r, t)$.
3. Identification des caustiques : Les caustiques quantiques sont identifiées comme les lieux de densité de rayons infinie où l'amplitude de la fonction d'onde est maximisée. Mathématiquement, elles correspondent à la coalescence de points critiques (cols) du fonctionnel d'action. Les auteurs résolvent les conditions $\nabla_k S = 0$ et $\det(\nabla_k \nabla_k S) = 0$ pour dériver les trajectoires exactes de ces caustiques.
4. Excitation périodique : Pour étudier les effets hors équilibre, les auteurs introduisent des protocoles dépendants du temps où le paramètre de couplage $J_3$ est modulé périodiquement. Deux protocoles d'excitation spécifiques sont analysés :
   • Des coups de type $\delta$ appliqués périodiquement.
   • Des impulsions en onde carrée alternant entre deux forces de couplage.
     L'évolution est traitée en utilisant la théorie de Floquet, en construisant un hamiltonien de Floquet effectif ($H_{eff}$) à partir de l'opérateur d'évolution temporelle sur une période.

Principales contributions et résultats

1. Anisotropie spatiale en 2D : Contrairement aux systèmes 1D, la dynamique des quasi-particules dans le modèle de Kitaev 2D présente une forte anisotropie spatiale. L'amplitude de la fonction d'onde se propage le long de directions spécifiques dans le plan du réseau, déterminées par le rapport des paramètres de couplage ($J_2/J_1$). À mesure que ce rapport varie, la direction de propagation maximale change, et l'amplitude décroît à mesure que l'on s'éloigne de cet axe privilégié.
2. Trajectoires exactes des caustiques et limite de Lieb-Robinson : En résolvant les conditions de coalescence pour les points critiques, les auteurs dérivent des expressions analytiques exactes pour les enveloppes des caustiques. Ces enveloppes définissent la vitesse maximale de propagation de l'information, restaurant efficacement la limite de Lieb-Robinson pour le système 2D.
   • Dans la phase gappée ($J_3 \gg J_1, J_2$), les frontières du cône de lumière sont linéaires : $x_1 = 2J_1t$ et $x_2 = 2J_2t$.
   • Dans la phase sans gap ($J_3 = J_1 = J_2$), les frontières forment également une structure en cône de lumière où l'impulsion $k$ reste réelle uniquement dans des limites spécifiques.
3. Effet de l'excitation périodique : L'introduction d'une excitation dépendante du temps modifie fondamentalement la structure des caustiques.
   • Sous les protocoles à coup $\delta$ et à onde carrée, l'étalement caractéristique en cône de lumière observé dans le système statique disparaît.
   • Au lieu de cela, l'amplitude de la fonction d'onde présente un étalement local qui décroît à chaque période d'excitation.
   • La quasi-particule devient confinée (localisée) le long de lignes spécifiques (par exemple, $x_1 = 0$) plutôt que de se propager librement. Cette localisation persiste plus longtemps à mesure que la fréquence d'excitation tend vers zéro ($\omega \to 0$).

Signification et affirmations
L'article prétend fournir le premier cadre analytique pour étudier la dynamique d'une seule quasi-particule dans le modèle de Kitaev 2D en utilisant des fonctions de catastrophes. La signification principale réside dans :

• La démonstration que les systèmes intégrables 2D supportent des structures de type cône de lumière qui servent de caustiques quantiques, mais avec une dépendance directionnelle (anisotropie) absente en 1D.
• La fourniture de solutions analytiques exactes pour les enveloppes des caustiques, qui correspondent à la limite de Lieb-Robinson en 2D.
• La démonstration que des excitations périodiques externes peuvent supprimer complètement la structure en cône de lumière, conduisant à la localisation des quasi-particules. Cela suggère que la fréquence d'excitation est un paramètre de contrôle critique pour manipuler la propagation de l'information quantique et les échelles de temps de mesure dans de tels systèmes.

Les auteurs notent que, bien que la réalisation du modèle de Kitaev dans les réseaux optiques fasse l'objet de recherches en cours, leurs résultats théoriques offrent une base pour tracer la propagation des quasi-particules et potentiellement contrôler le mouvement d'une seule quasi-particule via des champs externes dans des configurations expérimentales futures.