Wigner crystallization in Bernal bilayer graphene

Résumé technique : Cristallisation de Wigner dans le graphène bicouche de Bernal

Énoncé du problème
L'article étudie la nature et la stabilité de l'état de cristal de Wigner (CW) dans le graphène bicouche de Bernal (BBG) sous l'influence d'un champ de déplacement perpendiculaire ($U$). Alors que des expériences récentes sur des multicouches de graphène torsadées et non torsadées ont révélé des phases corrélées riches, l'interaction spécifique entre la structure de bande unique du BBG (spécifiquement la dispersion en « chapeau mexicain » induite par $U$), la courbure de Berry et la déformation trigonale (trigonal warping) sur la phase du CW reste une question théorique ouverte. Les auteurs visent à déterminer les frontières de phase du CW en fonction de la densité électronique ($n$), du champ de déplacement ($U$) et de la température ($T$), ainsi qu'à caractériser les propriétés de l'état fondamental, notamment le magnétisme orbital et l'ordre de vallée, en les distinguant des phases de CW conventionnelles.

Méthodologie
Les auteurs emploient une description de type oscillateur harmonique (OH) à particule unique pour modéliser l'état de CW. Dans ce cadre, chaque électron du réseau de WC triangulaire est traité comme une particule confinée dans un puits de potentiel parabolique créé par la répulsion coulombienne de ses voisins. La force de ce confinement est déterminée par la densité électronique.

L'étude procède selon deux régimes théoriques principaux :

1. Modèle à symétrie rotationnelle (sans déformation trigonale) : Les auteurs négligent d'abord la déformation trigonale pour traiter la dispersion à basse énergie comme étant rotationnellement symétrique autour des points $K$ et $K'$. Ils dérivent une équation de Schrödinger effective à particule unique dans l'espace des impulsions, incorporant la connexion de Berry ($A(p)$) comme un potentiel vecteur. Cela permet des solutions analytiques pour estimer les spectres d'énergie, les fonctions d'onde et le critère de Lindemann pour la fusion.
2. Modèle de déformation trigonale : Les auteurs incluent ensuite les effets de la déformation trigonale, qui brise la symétrie rotationnelle et divise la dispersion en « chapeau mexicain » en « mini-vallées » discrètes (une poche centrale et trois poches latérales). Pour ce cas, ils effectuent une diagonalisation numérique de l'équation de Schrödinger bidimensionnelle sur une grille de l'espace des impulsions pour calculer le rapport de Lindemann et cartographier le diagramme de phase.

La stabilité du CW est évaluée à l'aide du critère de Lindemann, où la fusion se produit lorsque la largeur de la racine carrée de la moyenne de la fonction d'onde de l'électron dépasse une fraction critique ($\eta_c \approx 0,23$) de la constante de réseau. Les auteurs calculent également les intégrales d'échange pour estimer les températures d'ordre et évaluent le magnétisme via la valeur attendue de l'opérateur de moment angulaire orbital.

Contributions clés et résultats

• Dépendance au champ de déplacement et écrantage diélectrique : L'article établit qu'à champ de déplacement nul, l'état de CW est empêché dans le BBG en raison de la masse de bande s'annulant à basse énergie et d'un fort écrantage diélectrique interbande. Une phase de CW n'émerge que lorsqu'un champ de déplacement suffisant est appliqué pour ouvrir une bande interdite et aplatir le fond de la bande.
• Magnétisme orbital spontané (sans déformation trigonale) : En l'absence de déformation trigonale, les auteurs démontrent que la courbure de Berry induit une transition dans le moment angulaire orbital ($\ell$) des électrons du WC. À mesure que le champ de déplacement $U$ augmente, le flux de Berry enfermé par la fonction d'onde augmente. Lorsque ce flux dépasse une valeur critique ($\tilde{\Phi} > 1/2$), l'état fondamental passe de $\ell=0$ à $\ell=1$ (ou $\ell=-1$ dans la vallée $K'$). Cette transition se manifeste par un saut discontinu dans le magnétisme orbital de l'état de WC.
• Fusion doublement réentrante (avec déformation trigonale) : Lorsque la déformation trigonale est incluse, le diagramme de phase présente un comportement unique de « double réentrance » en fonction de la densité électronique à faible $U$.
  1. À très basses densités, les électrons se cristallisent dans les « poches latérales » de plus basse énergie.
  2. À mesure que la densité augmente, le WC fond en un liquide de Fermi (FL) à une densité critique $n_{side}^c$.
  3. À une densité plus élevée $n^*$, il devient énergétiquement favorable pour les électrons d'occuper la « poche centrale » à masse plus lourde, provoquant le recristallisation du système en un état de WC.
  4. Enfin, à une densité encore plus élevée $n_{central}^c$, ce WC de poche centrale fond à nouveau.
• Ordre de mini-vallée : La présence de déformation trigonale brise la symétrie rotationnelle, rendant le moment angulaire $\ell$ mal défini. Au lieu de cela, l'état fondamental est caractérisé par un motif spatial non trivial d'« ordre de mini-vallée ». S'appuyant sur des travaux récents (Réf. [102]), les auteurs prédisent que les électrons se polarisent dans des mini-vallées spécifiques selon un motif de rayures alternées sur le réseau de WC.
• Frontières de phase : Les auteurs fournissent des estimations numériques des frontières de phase. La densité maximale pour l'état de WC dans le BBG est estimée de l'ordre de $1,2 \times 10^{10} \text{ cm}^{-2}$, se produisant à un champ de déplacement d'environ $0,7 \text{ V/nm}$. La température de fusion maximale est estimée de l'ordre de $10 \text{ K}$ (diminuant avec la densité), la région de réentrance de la poche centrale ayant une température maximale plus faible de $\sim 1 \text{ K}$.

Signification et affirmations
L'article prétend fournir un diagramme de phase théorique complet pour l'état de WC dans le BBG, soulignant comment la géométrie quantique non triviale (courbure de Berry) et la rupture de symétrie du réseau (déformation trigonale) modifient fondamentalement la nature du cristal de Wigner.

• Effets orbitaux : Le travail identifie un mécanisme pour le magnétisme orbital spontané dans un WC piloté purement par la courbure de Berry et les champs de déplacement, distinct de l'ordre basé sur le spin.
• Comportement de phase complexe : La prédiction de la fusion « doublement réentrante » est présentée comme une caractéristique novatrice découlant de la compétition entre différentes poches de bandes (centrale vs latérales) possédant des masses effectives et des énergies distinctes.
• Contexte expérimental : Les auteurs notent que la plage de densité prédite pour le WC dans le BBG ($n \lesssim 10^{10} \text{ cm}^{-2}$) est inférieure aux densités où des états isolants ont été observés lors d'expériences récentes, suggérant que les observations expérimentales actuelles dans le BBG pourraient être dominées par le désordre ou correspondre à des régimes différents. Ils suggèrent que le comportement de réentrance et l'ordre de mini-vallée pourraient être plus aisément observables dans le graphène à empilement rhomboédrique avec un nombre de couches ($N$) plus élevé, où la structure de bande est plus plate et le flux de Berry plus grand ($N\pi$).
• Complémentarité : Le papier est explicitement déclaré complémentaire à un travail compagnon (Réf. [1]) qui se concentre sur l'ordre de spin ; ce présent article se concentre sur les effets orbitaux et le diagramme de phase détaillé incluant la déformation trigonale.

Les auteurs maintiennent un ton modeste concernant les limites de l'approximation OH, reconnaissant qu'elle peut surestimer la stabilité de la phase de WC par rapport à des calculs de Monte Carlo quantique plus rigoureux, et que la limite de l'« Einstein phonon » utilisée ici ne peut reproduire les phases possédant des nombres de Chern non triviaux (tels que le cristal de l'effet Hall anomal quantique).