Site Basis Excitation Ansatz for Matrix Product States

Cet article présente une méthode efficace, appelée ansatz d'excitation de base de site (SBEA), pour calculer les spectres d'excitation élémentaire des systèmes quantiques unidimensionnels via des états de produit matriciel infinis, en exploitant une petite base non orthogonale et une construction d'excitations de Wannier qui permet d'obtenir avec précision la dispersion des magnons sans imposer de condition de jauge contraignante.

L'essentiel

Imaginez que vous essayez de comprendre comment une foule de personnes (les atomes d'un matériau) bouge et danse ensemble. En physique quantique, ces "danseurs" sont si nombreux et si complexes qu'il est impossible de les suivre un par un. Les scientifiques utilisent donc des outils mathématiques puissants pour décrire leurs mouvements collectifs.

Ce papier, écrit par Steven R. White, propose une nouvelle façon de calculer comment ces foules réagissent lorsqu'on les secoue légèrement pour créer une "vague" d'énergie (une excitation). Voici l'explication simplifiée, avec quelques analogies pour rendre les choses claires.

1. Le Problème : Trouver les vagues dans la foule

Imaginez une chaîne de dominos ou une foule de gens qui se tiennent par la main. Si vous poussez le premier, une onde de choc traverse la chaîne. En physique, on appelle cela une excitation (ici, un "magnon").

• Le défi est de prédire exactement comment cette onde se déplace, quelle vitesse elle a, et quelle énergie elle nécessite, pour toutes les directions possibles.
• Les méthodes anciennes étaient soit trop lentes (comme essayer de calculer chaque pas de chaque danseur), soit trop complexes à mettre en place.

2. La Solution : L'Ansatz SBEA (Le "Kit de Base")

L'auteur propose une méthode appelée SBEA (Site Basis Excitation Ansatz). Voici comment cela fonctionne avec une analogie simple :

L'ancienne méthode (l'approche "Sur mesure") :
Imaginez que vous vouliez connaître la vitesse d'une vague pour chaque direction possible (Nord, Sud-Est, etc.). Avec les anciennes méthodes, vous deviez construire un modèle mathématique spécial, unique, pour chaque direction. C'est comme si vous deviez fabriquer un nouveau type de chaussure pour chaque angle de marche possible. C'est efficace, mais fastidieux et long.

La nouvelle méthode SBEA (l'approche "Boîte à outils") :
Au lieu de fabriquer une chaussure pour chaque angle, l'auteur dit : "Et si on créait une petite boîte à outils avec seulement 7 paires de chaussures de base ?"

1. Le Kit de Base : Il utilise une astuce mathématique (un calcul rapide appelé "Lanczos") pour trouver les 7 mouvements de base les plus probables d'un seul danseur dans la foule.
2. Le Mélange : Pour n'importe quelle direction (momentum), il suffit de mélanger ces 7 paires de chaussures avec des proportions différentes. C'est comme un mélangeur de couleurs : avec seulement 7 couleurs de base, vous pouvez créer n'importe quelle teinte.
3. Le Résultat : Au lieu de faire un gros calcul pour chaque direction, vous faites un petit calcul rapide (comme résoudre une équation simple) pour trouver le bon mélange. C'est extrêmement rapide et précis.

3. L'astuce secrète : Ne pas être trop "propre"

En mathématiques, on aime souvent que les choses soient "orthogonales" (c'est-à-dire parfaitement indépendantes, comme des axes X et Y qui ne se touchent pas). Les méthodes précédentes forçaient ces mouvements de base à être parfaitement indépendants.

L'auteur a découvert une astuce géniale : il ne faut pas les forcer à être indépendants !

• L'analogie : Imaginez que vous essayez de décrire un paysage. Si vous forcez vos pinceaux à ne jamais se toucher, vous aurez du mal à peindre un ciel nuageux. Mais si vous laissez vos pinceaux se chevaucher un peu (ce qui crée un mélange), vous pouvez peindre le ciel beaucoup plus vite et avec plus de naturel.
• En laissant les mouvements de base se chevaucher (non-orthogonaux), la méthode devient beaucoup plus puissante et converge plus vite. C'est contre-intuitif, mais cela fonctionne parfaitement ici.

4. La "Fonction de Wannier" : Les briques locales

Le papier introduit aussi un concept cool appelé les excitations de Wannier.

• L'analogie : Imaginez que vous avez une onde qui traverse toute la chaîne (comme une vague à la plage). C'est difficile à étudier car elle est partout.
• La méthode SBEA permet de reconstruire cette vague en utilisant des "briques" locales. Imaginez que vous prenez une petite pierre (une excitation locale) et que vous la déplacez d'un site à l'autre.
• L'auteur montre qu'en combinant intelligemment ces "pierres" locales, on peut reconstruire exactement n'importe quelle vague, peu importe sa direction. C'est comme dire : "Pour comprendre la marée, il suffit de bien comprendre comment une seule goutte d'eau bouge, puis de répéter ce mouvement partout."

En résumé

Steven White a inventé une méthode pour étudier les mouvements des atomes dans les matériaux magnétiques qui est :

1. Plus simple : Pas besoin de calculs complexes pour chaque direction.
2. Plus rapide : On calcule une fois une petite "boîte à outils" de mouvements, et on l'utilise pour tout.
3. Plus précise : En acceptant que les mouvements se mélangent un peu (au lieu de les séparer strictement), on obtient des résultats plus justes.

C'est comme passer d'une méthode où l'on construit une maison brique par brique pour chaque pièce, à une méthode où l'on utilise un kit de préfabriqués intelligents pour construire n'importe quelle maison en quelques minutes. C'est une avancée majeure pour comprendre les matériaux quantiques, surtout ceux qui sont difficiles à modéliser.

Résumé technique

1. Problématique

Le calcul des spectres d'excitation élémentaires (dispersion) des systèmes quantiques à un dimension (1D) est crucial pour la compréhension théorique et la comparaison avec l'expérience. Bien que la méthode DMRG (Density Matrix Renormalization Group) soit très efficace pour les états fondamentaux, l'obtention d'une relation de dispersion complète pour les états excités reste difficile.

Les approches existantes basées sur les états produits matriciels (MPS), telles que l'Ansatz d'Excitation (EA) standard, présentent plusieurs limitations :

• Coût computationnel : Chaque vecteur d'onde $k$ nécessite une diagonalisation complète et séparée, ce qui est inefficace pour obtenir une courbe de dispersion continue.
• Complexité et convergence : La méthode VUMPS (Variational Uniform MPS) pour obtenir l'état fondamental infini peut échouer à converger selon les matrices initiales.
• Choix de jauge : Les formulations précédentes imposent souvent une condition de jauge (gauche-orthonormale) qui, bien que simplifiant l'orthogonalité, peut nuire à la convergence et à l'efficacité de la base d'excitation.

2. Méthodologie : L'Ansatz d'Excitation en Base de Site (SBEA)

L'auteur propose une variation simple et efficace de l'EA, appelée Site Basis Excitation Ansatz (SBEA). La méthode se décompose en trois étapes principales :

A. Obtention de l'état fondamental infini (iMPS)

Au lieu d'utiliser des algorithmes complexes comme VUMPS, l'auteur propose une méthode simplifiée basée sur le DMRG sur système fini :

1. On effectue un DMRG sur une chaîne ouverte longue ($N \gg \xi$, où $\xi$ est la longueur de corrélation).
2. On insère un nouveau site au centre de la chaîne et on optimise ce site via une diagonalisation de Lanczos (similaire à une étape de DMRG à un site).
3. Cela permet de construire un MPS uniforme infini de la forme $\dots AAA \dots$ (ou $\dots \lambda \Gamma \lambda \Gamma \dots$) avec une grande précision, en évitant les problèmes de convergence initiaux.

B. Construction de la base d'excitation minimale

Au lieu de déterminer un tenseur d'excitation $B(k)$ pour chaque moment $k$, le SBEA introduit une base minimale de tenseurs $\{B_\alpha\}$ :

1. Ansatz de la base : On suppose que tout tenseur d'excitation $B(k)$ peut être écrit comme une combinaison linéaire de $N_\alpha$ tenseurs de base : $B(k) = \sum c_\alpha^k B_\alpha$.
2. Génération de la base : Les tenseurs $B_\alpha$ sont obtenus en effectuant une seule diagonalisation de Lanczos sur un site unique, en utilisant un Hamiltonien effectif. On conserve les $N_\alpha$ états propres de plus basse énergie.
3. Pré-troncature : Pour limiter le coût, on tronque l'environnement des tenseurs avant la diagonalisation, en ne gardant que les plus grandes valeurs singulières (jusqu'à une dimension de liaison réduite $\chi_{tr}$).
4. Absence de condition de jauge : Contrairement aux travaux antérieurs, aucune condition de jauge gauche n'est imposée. Cela permet de conserver l'orthogonalité naturelle (ou plutôt la non-orthogonalité) des états, ce qui s'avère crucial pour obtenir une base compacte et efficace. Imposer une jauge gauche rendrait les états locaux trop énergétiques et moins utiles.

C. Théorie des bandes non orthogonale

Une fois la base $\{B_\alpha\}$ définie, on calcule les matrices de recouvrement $O_{\alpha\alpha'}(d)$ et de Hamiltonien $H_{\alpha\alpha'}(d)$ dans l'espace réel (où $d$ est la distance entre sites).

• Grâce à la nature gappée du système, ces matrices décroissent exponentiellement et peuvent être tronquées.
• Pour chaque moment $k$, on effectue une transformation de Fourier pour obtenir $\tilde{O}(k)$ et $\tilde{H}(k)$.
• La dispersion $E(k)$ est obtenue en résolvant un problème aux valeurs propres généralisé de petite taille ($N_\alpha \times N_\alpha$) :
  $$\tilde{H}(k) \mathbf{c} = E(k) \tilde{O}(k) \mathbf{c}$$
  Cela équivaut à une théorie des bandes où les orbitales sur différents sites ne sont pas orthogonales.

D. Fonctions de Wannier

L'article montre également comment construire des excitations de Wannier (états localisés) à partir des états de moment $k$. Ces états, orthonormés et localisés, permettent de reconstruire exactement la bande de magnons et pourraient servir de base pour des Hamiltoniens effectifs multi-magnons.

3. Contributions Clés

• Simplification algorithmique : Remplacement des méthodes d'optimisation d'état fondamental infini complexes par une procédure simple basée sur l'insertion de site dans un système fini.
• Efficacité computationnelle : La méthode SBEA réduit le problème de la dispersion à la diagonalisation d'une petite matrice pour chaque $k$, après un pré-calcul coûteux mais unique de la base.
• Innovation sur la jauge : Démonstration que l'imposition d'une condition de jauge gauche (standard dans l'EA) est préjudiciable. Le maintien d'une base non orthogonale permet une troncature drastique de la taille de la base ($N_\alpha$) tout en maintenant une haute précision.
• Construction de Wannier : Introduction d'une méthode pour générer des excitations localisées (Wannier) dans le contexte des MPS, offrant un pont vers les modèles de type Hubbard pour les interactions multi-particules.

4. Résultats

Les tests ont été effectués sur la chaîne de Heisenberg antiferromagnétique de spin $S=1$ (modèle de Haldane), un système gappé bien compris.

• Précision : La méthode reproduit la dispersion du magnon avec une excellente précision. Pour $N_\alpha = 7$ (avec une troncature $\chi_{tr}=8$), le gap de Haldane est calculé à $\Delta/J = 0.4107$, très proche de la valeur exacte $0.410479$.
• Convergence : La dispersion converge rapidement vers la solution exacte à mesure que l'énergie maximale incluse dans la base ($E_{max}$) augmente.
• Efficacité : Le calcul du spectre pour 500 valeurs de $k$ ne prend que quelques secondes une fois les matrices de recouvrement et de Hamiltonien calculées (moins d'une minute au total).
• Limites : Comme l'EA standard, la méthode a du mal à décrire les états à deux magnons ou les régions où la bande de magnon unique entre dans le continuum à deux magnons ($k < 0.72\pi$), nécessitant des ansatz plus complexes (ex: ansatz à deux particules).

5. Signification et Impact

Ce travail propose une refonte pratique et efficace de la méthode d'excitation par MPS.

• Accessibilité : En éliminant la complexité de VUMPS et les conditions de jauge restrictives, le SBEA abaisse le seuil d'adoption pour les chercheurs souhaitant étudier les spectres d'excitation.
• Extensibilité : La méthode est conçue pour s'étendre facilement à des cellules unitaires plus grandes et à des systèmes bidimensionnels (2D), où les méthodes actuelles atteignent souvent leurs limites computationnelles.
• Fondement théorique : La capacité à construire des fonctions de Wannier localisées ouvre la voie à la construction de Hamiltoniens effectifs renormalisés pour étudier les interactions complexes entre plusieurs quasiparticules, un défi majeur en physique de la matière condensée.

En résumé, le SBEA combine la précision des méthodes variationnelles modernes avec une simplicité algorithmique inédite, rendant le calcul des spectres d'excitation dans les systèmes 1D (et potentiellement 2D) plus rapide, plus robuste et plus facile à mettre en œuvre.