Chiral solitary waves in a nonlinear topological insulator model

1. Problématique et Contexte

Les isolants topologiques (TI) sont des matériaux qui conduisent le courant sur leurs bords tout en restant isolants dans leur volume, une propriété protégée par des invariants topologiques (comme le nombre de Chern). Bien que les TI linéaires soient bien compris et observés dans divers systèmes (photonique, acoustique, etc.), la réalisation de systèmes non linéaires supportant des ondes cohérentes en mouvement reste un défi majeur.

Le problème central identifié par les auteurs réside dans la gestion des pertes d'énergie dans les réseaux discrets fortement non linéaires. Les modèles traditionnels utilisant des non-linéarités de type Kerr (sur site) souffrent souvent d'effets liés à la barrière de Peierls-Nabarro (PN). Cette barrière, due au manque d'invariance par translation dans les réseaux discrets, crée une dépendance de l'énergie locale par rapport à la position du pic de l'onde (sur site ou entre sites). Cela conduit généralement à une radiation importante, à la perte de vitesse ou à la désintégration des solitons en mouvement, empêchant la formation d'états de bord cohérents et stables.

L'objectif de cette étude est de proposer et d'examiner un modèle d'isolant de Chern non linéaire capable de supporter des états de bord solitaires robustes et en mouvement, en s'inspirant de l'équation intégrable d'Ablowitz-Ladik.

2. Méthodologie

Les auteurs comparent deux modèles non linéaires basés sur le modèle linéaire de Haldane (un modèle de réseau en nid d'abeille avec brisure de symétrie d'inversion et de renversement du temps) :

• Modèle NH (Nonlinear Haldane) : Une extension standard avec une non-linéarité de type Kerr sur site (interaction locale $|\psi|^2\psi$).
• Modèle ALH (Ablowitz-Ladik-Haldane) : Un modèle inspiré de l'équation d'Ablowitz-Ladik, où la non-linéarité est de type « non-local » ou modulant l'opérateur de saut. L'équation de mouvement prend la forme $i\dot{\psi} = (1 + \sigma|\psi|^2)H_{lin}\psi$. Cette formulation préserve l'intégrabilité dans la limite 1D et modifie l'interaction des voisins les plus proches par l'intensité locale.

Approche numérique et analytique :

• Simulations dynamiques : Les auteurs utilisent une intégration Runge-Kutta d'ordre 4 pour faire évoluer des états initiaux (modes de bord linéaires) dans les deux modèles.
• Analyse topologique locale : Pour caractériser la topologie dans des systèmes finis et non linéaires (où le nombre de Chern global est difficile à définir), ils utilisent le localisateur spectral (spectral localizer). Cela permet de calculer un nombre de Chern local ($C_L$) en fonction de la position et de l'énergie.
• Étude des interactions : Ils simulent des collisions entre des ondes solitaires en mouvement (modes de bord) et des solitons de surface stationnaires.
• Recherche de solutions stationnaires : Utilisation d'une méthode de renormalisation spectrale pour trouver des solitons de surface stationnaires dans le modèle ALH.

3. Contributions Clés

1. Introduction du modèle ALH en 2D : C'est la première étude considérant la non-linéarité d'Ablowitz-Ladik dans un système d'isolant topologique 2D avec un mouvement non réciproque.
2. Suppression de la barrière de Peierls-Nabarro : Le modèle ALH démontre une capacité supérieure à supporter des courants de bord en mouvement par rapport au modèle NH, atténuant les pertes par radiation dues à la discrétisation du réseau.
3. Caractérisation topologique non linéaire : L'application réussie du localisateur spectral pour identifier la topologie non triviale dans des régimes non linéaires dynamiques, reliant le nombre de Chern local à la présence de modes de bord stables.
4. Découverte d'interactions inélastiques : La mise en évidence de dynamiques de collision complexes entre solitons mobiles et stationnaires dans un contexte topologique.

4. Résultats Principaux

A. Dynamique des ondes de bord

• Modèle NH (Kerr) : Les modes de bord fortement non linéaires se désintègrent rapidement. Une radiation massive est observée, et l'amplitude du pic diminue jusqu'à ce que l'état cohérent disparaisse, laissant place à une « mer » de radiation. Cela est attribué à l'effet de la barrière de Peierls-Nabarro et à la modification effective de la masse du système.
• Modèle ALH : En revanche, les courants de bord restent cohérents même à des puissances élevées. Au lieu de se désintégrer, le mode initial évolue vers un train d'ondes solitaires (wave train) qui se propage de manière robuste autour du réseau. Les solitons individuels conservent leur forme et leur vitesse sur de longues distances.

B. Topologie et Stabilité

• Le localisateur spectral confirme que la région intérieure du réseau possède un nombre de Chern local non nul ($C_L = 1$) dans la bande interdite, même en présence de non-linéarités fortes dans le modèle ALH.
• Une corrélation est établie : la présence d'un nombre de Chern local non nul est une condition nécessaire (mais non suffisante) pour l'existence de solitons de bord chirals stables. Lorsque le nombre de Chern local devient nul (par exemple en modifiant le terme de masse $M$), la radiation augmente considérablement.

C. Interactions Soliton-Soliton

Les auteurs étudient la collision entre un soliton de bord en mouvement (modèle ALH) et un soliton de surface stationnaire (trouvé dans la bande interdite supérieure) :

• Faible amplitude (A=1) : L'interaction est quasi-élastique. Le soliton mobile traverse le soliton stationnaire avec peu de changement, bien qu'une micro-mouvement cyclique soit observé sur le soliton stationnaire.
• Amplitude modérée (A=2) : L'interaction est plus forte. Le soliton stationnaire est repoussé de la surface vers le volume (bulk) du réseau.
• Forte amplitude (A=4) : L'interaction est inélastique et fatale. Le soliton mobile de haute puissance disloque complètement le soliton stationnaire, le radiant dans le volume du réseau.
• Conclusion sur les interactions : Ces résultats suggèrent un mécanisme de « nettoyage » de surface, où des ondes non linéaires fortes peuvent éliminer des défauts ou des modes localisés indésirables à l'interface.

5. Signification et Perspectives

Ce travail ouvre de nouvelles voies pour la réalisation et l'application des isolants de Chern non linéaires :

• Robustesse des états de bord : Il démontre que des couplages non linéaires légèrement non locaux (type Ablowitz-Ladik) sont essentiels pour surmonter les limitations des réseaux discrets et maintenir la cohérence des états de bord topologiques.
• Applications potentielles :
  • Communication optique solitonique : Utilisation de ces états pour le transport d'information robuste dans des guides d'ondes topologiques.
  • Nettoyage de surface : Le mécanisme d'interaction inélastique pourrait être exploité pour éliminer sélectivement des modes localisés ou des défauts dans des cristaux photoniques, en particulier dans des zones difficiles d'accès, en utilisant simplement la propagation d'ondes topologiques.
• Fondements théoriques : L'étude souligne la nécessité de développer de nouvelles théories de stabilité pour les isolants topologiques non linéaires, au-delà de la simple correspondance bulk-edge linéaire.

En résumé, l'article établit que le modèle Ablowitz-Ladik-Haldane est un candidat prometteur pour la réalisation expérimentale d'isolants topologiques non linéaires capables de supporter des solitons chirals en mouvement, offrant des perspectives concrètes pour le contrôle de la lumière et des ondes dans des milieux structurés.