Spin-imbalanced fermion on a dynamic lattice

Résumé Technique : Fermions à déséquilibre de spin sur un réseau dynamique

Énoncé du Problème
Cette étude examine les propriétés magnétiques de l'état fondamental d'un modèle de réseau dynamique de fermions de spin-1/2 en une dimension (1D). Alors que les systèmes 1D conventionnels (tels que le modèle de Fermi-Hubbard étendu) ont établi que les vecteurs d'ordre des ondes de densité de spin (SDW) suivent la relation $k = 2k_F$ dans les systèmes équilibrés, le comportement de ces systèmes sous un déséquilibre de spin des fermions ($N_\uparrow \neq N_\downarrow$) au sein d'un réseau dynamique reste largement inexploré. Plus précisément, l'article traite de la manière dont l'interaction entre les fermions itinérants et les spins localisés centrés sur les liaisons, médiée par un saut dépendant du spin, influence la modulation magnétique lorsque la population de fermions est partiellement polarisée.

Méthodologie
Les auteurs utilisent la méthode de la Densité de Matrice de Renormalisation (DMRG) pour simuler l'état fondamental du système. Le Hamiltonien du modèle décrit des fermions itinérants sautant entre les sites du réseau avec une amplitude $t$, laquelle est dynamiquement modulée par des spins localisés résidant sur les liaisons via un terme de couplage $\alpha$. Les spins localisés sont soumis à un champ longitudinal $\Delta$ et un champ transverse $\beta$. Les fermions itinérants subissent un champ de Zeeman externe $h$ (contrôlant le déséquilibre de spin) et une interaction de Hubbard répulsive sur site $U$.

Les simulations se concentrent sur le régime quasi-adiabatique ($\beta \ll t$), en fixant $\beta = 0,02$ et $t=1$. Des tailles de système allant de $L=24$ à $L=60$ sont analysées avec des conditions aux limites ouvertes. Les observables clés incluent le facteur de structure de spin localisé $S_z(k)$, la fonction de corrélation dans l'espace réel $S_z(r)$, et les fonctions de distribution de quantité de mouvement résolues en spin $n_\sigma(k)$. Une mise à l'échelle de taille finie est utilisée pour extrapoler les résultats vers la limite thermodynamique ($L \to \infty$).

Contributions Clés et Résultats

1. Identification des Phases Magnétiques :
   L'étude cartographie le diagramme de phase, identifiant deux phases primaires : une phase paramagnétique (PM) et une phase d'onde de densité de spin (SDW). La transition entre ces phases est pilotée par la compétition entre le couplage de saut dépendant du spin ($\alpha$) et le champ longitudinal ($\Delta$). Lorsque $\alpha$ et $\Delta$ sont comparables, des ordres SDW spatialement modulés émergent. Lorsqu'un paramètre domine, le système entre dans une phase PM caractérisée par un alignement de spin uniforme et l'absence de corrélations magnétiques à longue portée.

2. Phénomènes de Double Ordre SDW (SDW-I et SDW-II) :
   Dans le régime de spin partiellement polarisé, les auteurs découvrent deux motifs d'ordre SDW distincts, qui dévient fondamentalement de l'attente standard $k=2k_F$ pour les systèmes équilibrés :

   • SDW-I : Le vecteur d'onde d'ordre $k_{max}$ présente une dépendance linéaire vis-à-vis de la densité de spin majoritaire ($\rho_\uparrow$), suivant $k_{max} \approx 2\pi\rho_\uparrow$. Cela correspond au vecteur de mise en réseau (nesting) de la surface de Fermi des spins up.
   • SDW-II : Inversement, $k_{max}$ présente une dépendance linéaire vis-à-vis de la densité de spin minoritaire ($\rho_\downarrow$), suivant $k_{max} \approx 2\pi\rho_\downarrow$. Cela correspond au vecteur de mise en réseau de la surface de Fermi des spins down.
     La transition entre ces deux phénomènes est contrôlée par l'ajustement des paramètres de couplage ($\alpha, \Delta$) et le facteur de remplissage.

3. Origine Microscopique :
   L'analyse des fonctions de distribution de quantité de mouvement résolues en spin révèle que le pic SDW dominant est dicté non pas par la densité de remplissage totale, mais par les occupations relatives des composantes de spin. L'émergence de ces doubles motifs est interprétée comme un analogue magnétique d'une transition de type Peierls, pilotée par la mise en réseau de la surface de Fermi spécifique aux secteurs résolus en spin des fermions itinérants.

4. Robustesse face à l'Interaction de Hubbard :
   L'inclusion d'interactions de Hubbard répulsives sur site ($U$) est trouvée pour renforcer la tendance vers la polarisation de spin des fermions (abaissant le champ critique requis pour une polarisation complète) en supprimant la double occupation. Cependant, l'interaction répulsive ne modifie pas qualitativement la nature des phases SDW-I et SDW-II ; ces motifs d'ordre distincts restent robustes dans leurs régions de paramètres respectives.

Signification
L'article affirme révéler une nouvelle voie pour ajuster les modulations magnétiques dans les systèmes corrélés unidimensionnels. En démontrant que le déséquilibre de spin des fermions peut conduire le système vers des phases SDW distinctes caractérisées par différents vecteurs de mise en réseau ($2k_{F\uparrow}$ vs $2k_{F\downarrow}$), ce travail enrichit la compréhension microscopique du magnétisme de réseau dynamique. Les résultats suggèrent que l'ordre magnétique dans ces systèmes n'est pas simplement une propriété statique, mais est dynamiquement couplé à la configuration des spins localisés et à la population de spin spécifique des fermions itinérants. Cela fournit un cadre théorique pour comprendre comment le déséquilibre de spin façonne les phases quantiques dans les systèmes de faible dimension au-delà de la limite conventionnelle de l'équilibre.