Higher Berry curvature, second Chern numbers and magnetoelectric coupling in crystalline insulators

Résumé technique : Courbure de Berry supérieure, nombres de Chern de second degré et couplage magnétoélectrique dans les isolants cristallins

Énoncé du problème
Les invariants topologiques sont essentiels pour classifier les états fondamentaux quantiques à plusieurs corps. Alors que la courbure de Berry (une 2-forme) et le nombre de Chern associé caractérisent avec succès des familles d'états fondamentaux de dimension zéro (par exemple, les états de Bloch dans la zone de Brillouin), des développements théoriques récents suggèrent une généralisation aux systèmes étendus spatialement. Dans ce cadre, les familles d'états à intrication de courte portée dans $d$ dimensions spatiales sont caractérisées par des formes différentielles de degré supérieur (par exemple, une 3-forme pour les systèmes 1D). Spécifiquement, pour les systèmes unidimensionnels paramétrés par un espace tridimensionnel, l'invariant pertinent est le nombre de Dixmier–Douady–Kapustin–Spodyneiko (DDKS), $\nu$, dérivé d'une courbure de Berry supérieure.

Un écart important existe dans le calcul explicite et manifestement quantifié de ces invariants supérieurs pour des systèmes génériques de fermions libres ou de systèmes interagissants. Bien qu'il existe des expressions analytiques pour certains modèles de fermions libres (comme le modèle Dirac 4D), les méthodes numériques pour calculer le second nombre de Chern ($ch_2$) en quatre dimensions reposent généralement sur la théorie de jauge sur réseau ou sur les réponses de champs magnétiques dans l'espace réel. Ces méthodes standards manquent souvent de « quantification manifeste », ce qui signifie qu'elles ne convergent vers des valeurs entières que dans la limite du continuum ou nécessitent un raffinement extensif du maillage, plutôt que de retourner des entiers par construction. De plus, l'interprétation physique du nombre DDKS dans les trois dimensions spatiales, particulièrement sa relation avec le couplage magnétoélectrique orbital et l'angle axion de Chern-Simons ($\theta_{CS}$), demeure une question ouverte.

Méthodologie
Les auteurs abordent ces problèmes en réécrivant le modèle de réseau d'un isolant de Chern 4D comme une famille de chaînes unidimensionnelles (1D) invariantes par translation, infinies, paramétrées par la zone de Brillouin tridimensionnelle ($T^3$).

1. Reformulation du modèle : En partant de l'Hamiltonien de Dirac 4D, les auteurs effectuent une transformée de Fourier partielle sur une composante de quantité de mouvement ($k_4$), résultant en une famille d'Hamiltoniens 1D $\{H_\lambda\}$ dépendant de paramètres $\lambda = (\lambda_1, \lambda_2, \lambda_3)$. Ces Hamiltoniens décrivent des chaînes 1D avec deux orbitales par cellule unitaire.
2. Cadre numérique (iMPS) : Les états fondamentaux de ces chaînes 1D sont approximés par des états de matrice de produits (MPS) infinis (iMPS). Les auteurs utilisent l'algorithme proposé dans la Réf. [13], qui généralise la méthode de Fukui-Hatsugai-Suzuki pour calculer la courbure de Berry supérieure.
3. Courbure de Berry supérieure discrète :
   • Le paramètre $T^3$ est discrétisé en une grille et triangulé en tétraèdres.
   • Pour chaque sommet, un état fondamental iMPS est obtenu via l'iDMRG.
   • Les fonctions de transition entre les tenseurs iMPS voisins sont calculées en utilisant les vecteurs propres dominants des matrices de transfert mixtes.
   • Une phase invariante par jauge $\phi(\Delta_2)$ est définie pour les faces triangulaires, pondérée par le spectre de Schmidt pour gérer les variations de la dimension de liaison (bond dimension).
   • La courbure de Berry supérieure $F(\Delta_3)$ pour un tétraèdre est définie comme la somme orientée de ces phases sur son bord.
   • Le nombre DDKS $\nu$ est obtenu en sommant $F(\Delta_3)$ sur l'ensemble de la triangulation de $T^3$, produisant un entier manifestement quantifié.
4. Analyse comparative : Les auteurs calculent le second nombre de Chern $ch_2$ pour le même modèle en utilisant une approche standard de théorie de jauge sur réseau (Réf. [34]) et comparent les diagrammes de phases résultants.
5. Couplage magnétoélectrique : L'étude examine la relation entre le nombre DDKS et l'angle axion de Chern-Simons $\theta_{CS}$, qui régit le couplage magnétoélectrique orbital dans les 3D. Ils comparent la variation de $\theta_{CS}$ (dérivée analytiquement de la seconde forme de Chern) avec la variation d'une « phase de Berry supérieure » $\theta_{HBP}$ (dérivée de la construction iMPS) à travers un cycle de pompage.

Contributions clés et résultats

• Calcul manifestement quantifié de $ch_2$ : Les auteurs démontrent que le nombre DDKS $\nu$, calculé via l'algorithme de courbure de Berry supérieure basé sur l'iMPS, génère un diagramme de phase topologique pour l'isolant de Chern 4D qui est exactement congruent au diagramme de phase du second nombre de Chern $ch_2$. Cet accord est maintenu à travers différents régimes de masse ($m$), y compris les régions où les calculs standards de la théorie de jauge sur réseau de $ch_2$ souffrent d'erreurs de grille finie significatives. La méthode iMPS produit des entiers exacts même sur des grilles grossières (par ex. $4^3$), alors que la méthode de jauge sur réseau nécessite des grilles beaucoup plus fines ($10^4$ points) pour converger.
• Lien théorique via la formule de Künneth : L'article établit un lien rigoureux entre le second nombre de Chern 4D et le nombre DDKS 3D en utilisant la formule de Künneth. En décomposant le tore 4D $T^4$ en $S^1 \times T^3$, les auteurs montrent que la seconde classe de Chern sur $T^4$ correspond au produit du générateur de $H^1(S^1, \mathbb{Z})$ et de la classe dans $H^3(T^3, \mathbb{Z})$ mesurée par $\nu$. Ainsi, pour cette décomposition spécifique, $ch_2 = \nu$.
• Couplage magnétoélectrique dans le modèle de Dirac 4D : Pour le modèle de Dirac 4D, les auteurs trouvent que le changement de l'angle axion de Chern-Simons $\Delta\theta_{CS}$ le long d'un cycle de pompage est reproduit précisément par le changement de la phase de Berry supérieure $\Delta\theta_{HBP}$. Cela suggère une correspondance directe entre la phase de Berry supérieure et la réponse magnétoélectrique dans les systèmes où les nombres de Chern premiers des tranches de la zone de Brillouin (BZ) sont nuls.
• Contre-exemple dans l'isolant de Hopf : Les auteurs testent cette correspondance sur l'isolant de Hopf, un modèle à deux bandes possédant un invariant de Hopf non trivial. Dans ce cas, bien que l'angle de Chern-Simons soit quantifié par l'invariant de Hopf, la phase de Berry supérieure calculée $\theta_{HBP}$ n'est ni quantifiée ni égale à $\theta_{CS}$. Les auteurs attribuent cette divergence à la nature « délicate » de l'invariant de Hopf, qui n'est strictement bien défini que pour les modèles à deux bandes dans l'espace des moments, une condition qui peut ne pas être préservée dans la description iMPS dans l'espace réel de la limite thermodynamique.

Signification et affirmations
L'article affirme fournir le premier calcul numérique, manifestement quantifié, du second nombre de Chern pour un modèle de réseau 4D en utilisant la courbure de Berry supérieure et l'iMPS. Cela valide le cadre théorique selon lequel les courbures de Berry supérieures peuvent servir de sondes pour les invariants topologiques de degré supérieur dans les familles d'états étendus spatialement.

Les auteurs positionnent leurs résultats dans le contexte de la géométrie quantique et des propriétés topologiques, notant que si la correspondance entre les phases de Berry supérieures et le couplage magnétoélectrique tient pour le modèle de Dirac 4D, elle n'est pas universelle, comme le montre le cas de l'isolant de Hopf. Ils soulignent que leur méthode est particulièrement adaptée aux systèmes interagissants où les descriptions de Bloch à un seul électron échouent, bien qu'ils se concentrent actuellement sur des tests de référence de fermions libres. Ce travail met en évidence le potentiel des dimensions synthétiques et de la commande quasi-périodique dans les plateformes expérimentales pour réaliser et sonder ces phénomènes topologiques de dimensions supérieures.

L'article conclut modestement, affirmant que la relation entre la courbure de Berry supérieure et la réponse physique dans les trois dimensions reste une question ouverte, particulièrement concernant l'universalité de la connexion avec $\theta_{CS}$. Il suggère que les phases de Berry supérieures peuvent capturer des informations géométriques distinctes des réponses de Chern-Simons standards, justifiant des recherches supplémentaires sur leur rôle dans le transport non linéaire et la physique des états fondamentaux.