Trimers in the Extended Hubbard Model

Résumé Technique : Instabilité des Trimères dans le Modèle de Hubbard Étendu

Énoncé du Problème
Ce travail étudie la stabilité de la phase isolante ferrimagnétique de Lieb au sein du modèle de Hubbard étendu sur une chaîne de trimères couplés bipartites à demi-remplissage (un électron par site). Alors que le modèle de Hubbard standard sur les réseaux bipartites avec saut de plus proche voisin et répulsion sur site $U$ est connu pour présenter un état ferrimagnétique robuste prédit par le théorème de Lieb, l'introduction d'une interaction coulombienne de plus proche voisin $V$ est censée induire des instabilités de charge. Le problème central abordé est l'évolution de cet état ferrimagnétique à mesure que $V$ augmente, en identifiant spécifiquement les conditions dans lesquelles la phase isolante homogène devient instable face à une séparation de phases (PS) et en caractérisant la nature des phases coexistantes résultantes.

Méthodologie
Les auteurs utilisent la méthode de la Densité de Matrice de Renormalisation (DMRG) pour résoudre le hamiltonien de Hubbard étendu pour des chaînes de trimères couplées avec des conditions aux limites ouvertes. Le système consiste en $L$ trimères, où chaque trimère contient trois sites nommés $A$, $B_1$ et $B_2$. L'hamiltonien comprend :

• L'amplitude de saut intra-trimère $t$ (fixée à $t=1$).
• L'amplitude de saut inter-trimère $t'$ (fixée à $t'=0,8$).
• La répulsion coulombienne sur site $U$.
• La répulsion coulombienne de plus proche voisin $V$.

Les calculs ont été effectués sur des chaînes de tailles allant jusqu'à $L=60$ (soit $N=180$ sites), utilisant des dimensions de liaison maximales entre 512 et 2048, avec une erreur de troncature de $\sim 1 \times 10^{-6}$. Pour sonder les propriétés de l'état fondamental, les auteurs ont analysé :

1. Les moyennes spatiales et les distributions locales de l'occupation particulaire ($\langle n_i \rangle$) et de l'occupation des doublons ($\langle d_i \rangle = \langle n_{i\uparrow} n_{i\downarrow} \rangle$).
2. Le spin total par trimère ($S/L$), calculé via les fonctions de corrélation de spin dans le secteur $S_z=0$.
3. La compressibilité ($\kappa$), dérivée des différences d'énergie de taille finie pour identifier les régimes incompressibles (gapés) versus compressibles (métalliques/séparés de phases).
4. Les métriques d'homogénéité, spécifiquement $\Delta_{max}$, définie comme l'écart maximal du nombre de particules dans n'importe quel trimère par rapport à la densité homogène de 3 électrons par trimère.

Contributions Clés et Résultats

• Analyse de la Limite Atomique ($t=t'=0$) : Les auteurs ont d'abord analysé les configurations statiques dans la limite atomique. Ils ont identifié qu'augmenter $V$ à partir de zéro favorise une transition d'un état standard à demi-remplissage (un électron par site) vers un état où les sites $A$ sont vides et les sites $B$ sont doublement occupés. Cette transition se produit à une valeur critique $V = U/4$, où l'énergie de l'état à demi-remplissage ($E/L = 4V$) est égale à l'énergie de l'état de transfert de charge ($E/L = U$).

• Le Cas $U=0$ : Pour une répulsion sur site nulle, toute valeur finie de $V$ induit une séparation de phases. Le système se sépare en deux phases distinctes : l'une riche en doublons sur le sous-réseau $B$ (isolante), et l'autre riche en doublons sur le sous-réseau $A$ (métallique). Le spin total par trimère est nul ($S/L = 0$) dans ce régime.

• Stabilité de la Phase Ferrimagnétique de Lieb ($U \neq 0$) :

  • Pour $V=0$, le système présente l'état isolant ferrimagnétique de Lieb attendu avec $S/L = 1/2$.
  • À mesure que $V$ augmente, le système reste dans une phase isolante homogène jusqu'à un point critique. Cependant, il y a un transfert de charge significatif des sites $A$ vers les sites $B$, augmentant la densité de doublons sur le sous-réseau $B$.
  • Les auteurs identifient une ligne de crossover près de $V \approx U/4$. En dessous de cette ligne, le système est homogène. Entre la ligne de crossover et la ligne de séparation de phases, le système reste homogène mais présente un caractère « riche en doublons » avec un transfert de charge accru, tout en maintenant l'ordre ferrimagnétique ($S/L = 1/2$).

• Séparation de Phases et Coexistence :

  • Au-dessus d'une ligne critique de séparation de phases (PS) (observée à $V \gtrsim U/4$), l'état homogène devient instable. Le système entre dans une région de coexistence où deux phases avec des occupancies de trimères distinctes existent spatialement.
  • Pour $U \gtrsim 1$ : La région de coexistence consiste en :
    1. Une phase isolante où les sites $B$ sont principalement occupés par des doublons (environ 4 électrons par trimère), correspondant à un état riche en doublons.
    2. Une phase ferromagnétique non saturée et métallique où la distribution de charge est plus uniforme mais ne respecte pas strictement la condition de demi-remplissage. Cette phase présente un moment magnétique net ($S/L > 0$), résultant de l'interaction entre l'itinérance de la charge et l'alignement ferrimagnétique des moments sur les sites $A$ et $B$.
  • Pour $U \lesssim 1$ : La région de coexistence ressemble au cas $U=0$, consistant en deux phases paramagnétiques ($S/L \approx 0$), l'une riche en doublons sur les sites $B$ et l'autre sur les sites $A$.

• Comportement à Haut $U$ : Les auteurs estiment que la ligne de crossover et la ligne de PS fusionnent vers $U \sim 16$. Pour $U > 16$, le système passe directement de la phase ferrimagnétique homogène à la région de séparation de phases, indiqué par un saut dans les moyennes spatiales plutôt que par un point d'inflexion.

Signification et Revendications
L'article affirme démontrer que la phase isolante ferrimagnétique de Lieb est remarquablement robuste contre la répulsion de plus proche voisin $V$, persistant bien au-delà du point où la densité de doublons augmente significativement. La principale contribution est la cartographie du diagramme de phase pour les chaînes de trimères couplés, distinguant une phase ferrimagnétique « riche en doublons » homogène d'une phase de séparation de phases ultérieure.

Les auteurs affirment que leurs résultats clarifient la nature de l'instabilité de l'état fondamental dans cette géométrie de réseau spécifique. Ils soulignent qu'au-dessus de la ligne de PS, le système supporte une coexistence unique d'une phase paramagnétique isolante et d'une phase ferromagnétique non saturée métallique pour des valeurs intermédiaires à élevées de $U$. Le travail fournit une estimation quantitative des lignes de PS et de crossover en utilisant la DMRG, confirmant que la transition est pilotée par la compétition entre le coût énergétique de la double occupation ($U$) et le gain d'énergie lié à l'évitement de la répulsion de plus proche voisin ($V$). Les auteurs concluent que ces découvertes fournissent une base pour de futures investigations sur les modèles dopés et les aspects topologiques de telles transitions de phase, sans proposer de nouvelles réalisations expérimentales spécifiques au-delà du cadre théorique.