Finite-temperature formation of magnetic plateaus and simplex liquid states on the frustrated ruby lattice

En utilisant des états de réseaux de tenseurs infinis optimisés par propagation de croyance, cette étude révèle que l'antiferromagnétique de Heisenberg de spin 1/2 frustré sur le réseau rubis forme des plateaux magnétiques stables hébergeant un nouvel « état de liquide simplex » avec une entropie résiduelle à basses températures, un processus qui se produit continûment sans transition de phase à mesure que le système se refroidit.

L'essentiel

Imaginez une piste de danse bondée où tout le monde essaie de trouver un partenaire, mais les règles de la danse sont délicates. C'est l'histoire d'un nouvel article sur un type spécial de matériau magnétique appelé le « réseau rubis ».

Voici le détail de ce que les scientifiques ont découvert, en utilisant des analogies simples :

La piste de danse délicate (Frustration géométrique)

Dans les aimants normaux, les petits spins atomiques agissent comme de petites aiguilles de boussole qui veulent toutes s'aligner dans la même direction (comme des soldats qui marchent au pas). Mais dans cette structure spécifique de « réseau rubis », la géométrie est si tordue que les spins ne peuvent pas tous être heureux en même temps. C'est comme un jeu de chaises musicales où il y a plus de chaises que de personnes, mais les chaises sont disposées de telle manière qu'il est impossible pour tout le monde de s'asseoir confortablement sans bousculer quelqu'un d'autre. C'est ce qu'on appelle la frustration géométrique.

D'habitude, lorsque l'on refroidit ces matériaux, ils se sentent frustrés et finissent par se « figer » dans un motif rigide et ordonné (comme un cristal) pour résoudre le problème. Mais les scientifiques voulaient voir ce qui se passe si on les refroidit de manière très lente et méticuleuse.

La magie de la « Propagation du message »

Pour comprendre cela, les chercheurs ont utilisé une méthode informatique puissante appelée Réseaux de tenseurs, plus précisément une technique de Propagation du message (Belief Propagation - BP).

Voyez la Propagation du message comme une rumeur qui se répand dans une foule. Au lieu de demander à chaque personne de la pièce ce qu'elle fait, vous posez la question à quelques personnes, qui le disent à leurs voisins, qui le disent à leurs propres voisins, et ainsi de suite. Finalement, tout le monde a une bonne idée de ce que fait l'ensemble du groupe sans avoir besoin de vérifier chaque personne individuellement. Les chercheurs ont utilisé cette mathématique de la « propagation de rumeurs » pour simuler le comportement de ces spins magnétiques à différentes températures, même lorsque le système est infiniment grand.

La surprise : Pas de « blocage », juste un « liquide »

Lorsqu'ils ont refroidi le système, ils s'attendaient à ce que les spins se figent soudainement en un cristal rigide et ordonné (une transition de phase). Au lieu de cela, ils ont trouvé quelque chose de beaucoup plus fluide.

À mesure que la température chutait, les spins ne se figeaient pas dans un motif unique. À la place, ils formaient un « État liquide de simplexes ».

• L'analogie : Imaginez un groupe de personnes à une fête. Au lieu que tout le monde se tienne selon une grille parfaite (un cristal), ils forment de petits groupes très soudés de trois (appelés « simplexes »). Ces groupes dansent ensemble, mais l'arrangement des groupes ne cesse de changer.
• Le résultat : Même à des températures très basses, le système reste désordonné. C'est un « liquide » de ces groupes de danseurs. Parce qu'il existe de nombreuses façons pour ces groupes de s'organiser, le système conserve beaucoup d'« entropie résiduelle » (une mesure du désordre). C'est comme avoir un jeu de cartes qui est parfaitement mélangé à chaque fois que vous le regardez, ne se fixant jamais sur un ordre spécifique.

Les plateaux magnétiques

Les chercheurs ont également activé un champ magnétique (comme un vent fort soufflant sur la piste de danse). À mesure qu'ils augmentaient la force du vent, les spins tentaient de s'aligner avec lui.

Au lieu de tourner de manière fluide, les spins restaient bloqués sur des « plateaux » spécifiques.

• L'analogie : Imaginez un escalier. À mesure que vous poussez le système plus fort, l'aimantation (le degré d'alignement) fait des bonds, puis reste plate pendant un certain temps (un plateau), puis fait un autre bond, et ainsi de suite.
• Ils ont découvert des « marches » stables où l'aimantation était exactement de 1/3, 1/2 ou 2/3 du maximum possible.
• Le rebondissement : Même sur ces « marches » plates, le matériau ne devenait pas un cristal rigide. Il restait dans cet état « liquide », avec simplement un alignement moyen spécifique.

Le pic « Lambda » et le basculement

Il y a eu un moment très intéressant près du milieu de l'escalier (le plateau de 1/2).

• L'analogie : Imaginez que la piste de danse est divisée en deux. D'un côté, les groupes dansent d'une certaine façon ; de l'autre, ils dansent d'une autre. À une température et une vitesse de vent spécifiques, toute la piste bascule soudainement d'un style de danse à un autre.
• Ce basculement n'était pas progressif. Cela a créé un pic aigu dans la « capacité thermique » (la quantité d'énergie que le système absorbe), ayant la forme de la lettre grecque Lambda (λ). Cela suggère qu'au bord de ces plateaux, le système est au seuil d'un changement majeur, piloté par les fluctuations quantiques.

L'essentiel à retenir

La découverte la plus importante est que ce système ne se « fige » jamais dans un cristal traditionnel.

Même lorsqu'il devient incroyablement froid, il reste dans un état désordonné, semblable à un liquide, rempli de nombreuses configurations possibles. Les scientifiques ont prouvé cela en montrant que la « capacité thermique » (une mesure de la façon dont le système réagit aux changements de température) reste lisse et continue. Si le système s'était figé en un cristal, il y aurait eu un pic net et saccadé indiquant une transition de phase. Au lieu de cela, il a coulé de manière fluide vers ce nouvel état exotique.

En bref : Les chercheurs ont utilisé une astuce mathématique ingénieuse de « propagation de rumeurs » pour démontrer qu'un matériau magnétique frustré ne se fige pas en un cristal solide lorsqu'il est refroidi. Au lieu de cela, il se transforme en un « liquide » de groupes de spins danseurs qui reste désordonné et plein de possibilités, même à des températures proches du zéro absolu.

Résumé technique

Résumé technique : Formation à température finie de plateaux magnétiques et d'états de liquide de simplexes sur le réseau de rubis frustré

Énoncé du problème
La frustration géométrique dans les systèmes quantiques est connue pour stabiliser des phases non conventionnelles, telles que les liquides de spins quantiques, qui évitent l'ordre magnétique traditionnel à basses températures. Cependant, comprendre la formation de ces états à des températures finies et leurs signatures thermodynamiques reste un défi important. Plus précisément, l'état fondamental à température nulle ($T=0$) de l'antiferromagnétique de Heisenberg de spin-1/2 sur le réseau de rubis avec des interactions de second voisin a été précédemment montré comme présentant des plateaux magnétiques fractionnés ($m_z = 0, 1/3, 1/2, 2/3$). Ces plateaux ont été interprétés comme des phases cristallines discrètes brisant les symétries du réseau. Des questions restaient en suspens concernant la dégénérescence exacte de l'état fondamental, la stabilité de ces plateaux à température finie, et si le système subit des transitions de phase conventionnelles pour les atteindre. Résoudre ces problèmes nécessite des simulations capables de faire le pont entre les résultats computationnels et les réalisations expérimentales où les fluctuations thermiques et quantiques sont en compétition.

Méthodologie
Les auteurs utilisent un ansatz d'état de réseau de tenseurs infini purifié (iTNS) pour représenter la racine carrée non normalisée de la matrice densité d'équilibre, $\rho^{1/2}(\beta) = \exp(-\beta H/2)$, où $\beta = 1/T$. L'étude se concentre sur le cas isotrope ($J_t = J_h = J_d = J$) du Hamiltonien de Heisenberg sur le réseau de rubis.

• Représentation du réseau : Les sommets des unités triangulaires sont regroupés en triangles (simplexes), projetant le réseau de rubis vers une géométrie de type nid d'abeille macroscopique avec un nombre de coordination $z=3$. L'iTNS utilise une cellule unitaire à deux sites, où chaque tenseur possède six indices physiques (représentant les trois spins de chaque triangle) et trois indices virtuels.
• Algorithme : Le système est initialisé à température infinie ($\rho^{1/2}(0) = I$) et refroidi via une évolution dans le temps imaginaire en utilisant une Trotterisation de second ordre du propagateur.
• Observables : Les quantités thermodynamiques sont calculées selon deux approches :
  1. Propagation de croyance (BP) : Une contraction BP standard est appliquée au réseau de norme pour calculer les valeurs d'espérance.
  2. BP corrigée par les boucles : Pour améliorer la précision au-delà de l'approximation de type champ moyen de la BP, les auteurs appliquent des corrections de grappes de premier ordre (corrections de boucles) à la densité de l'énergie libre. Cela permet de prendre en compte les corrélations autour des hexagones du réseau.
• Échelle : Les simulations utilisent du matériel GPU pour atteindre des dimensions de liaison jusqu'à $\chi \sim 200$, permettant une mesure précise des observables thermodynamiques dans la limite thermodynamique sur une large gamme de températures.

Résultats clés

1. État de liquide de simplexes à champ nul :

   • Lors du refroidissement depuis la température infinie, le système ne subit pas de transition de phase vers un ordre cristallin. À la place, il forme un « état de liquide de simplexes » désordonné caractérisé par une entropie résiduelle non nulle ($s_0 \approx \frac{1}{3}\lambda$, où $\lambda$ est l'entropie résiduelle des recouvrements de dimères sur un réseau nid d'abeille).
   • La capacité calorifique ($C_v$) reste finie et continue à toutes les températures, sans pics divergents indiquant une transition de phase.
   • L'analyse de la matrice densité réduite révèle une valeur propre principale $\lambda_1 \approx 1/3$, correspondant à un état où deux simplexes sont fortement appariés. Les autres valeurs propres correspondent à des états faiblement corrélés. Le système est un mélange désordonné parfait d'un nombre exponentiellement grand de configurations cristallines $|\psi_d\rangle$ impliquant des paires de simplexes de spins.
   • La densité d'énergie de cet état liquide ($e/J \approx -0.498$) correspond à l'énergie du fondamental de configurations cristallines spécifiques trouvées dans des études antérieures à $T=0$, mais le système reste dans un état liquide invariante par translation plutôt que de sélectionner un cristal unique.

2. Plateaux magnétiques et comportement induit par le champ :

   • En présence d'un champ magnétique externe, le système présente des plateaux magnétiques stables aux densités de magnétisation $m_z = 0, 1/3, 1/2$ et $2/3$.
   • Signatures thermodynamiques : La formation de ces plateaux se produit de manière fluide sans transitions de phase conventionnelles. La capacité calorifique reste continue, bien que des pics peu profonds soient observés.
   • Entropie résiduelle : Les plateaux hébergent des « états de liquide de simplexes » dotés d'une dégénérescence extensive. Le plateau $m_z = 1/2$, en particulier, présente une grande entropie résiduelle (environ le double de celle du plateau $m_z=0$).
   • Stabilité : Le plateau $m_z = 0$ est le plus stable face aux fluctuations thermiques. Le plateau $m_z = 1/3$ est moins stable, n'apparaissant que dans une plage de champ étroite à très basse température ($T/J \le 10^{-2}$). Le plateau $m_z = 2/3$ est observé dans les calculs d'état fondamental mais nécessite des températures plus basses que celles actuellement accessibles dans la simulation à température finie pour être pleinement résolu.
   • Commutation de phase : À l'interface entre les plateaux $m_z = 1/2$ et $m_z = 1/3$, le système présente une commutation pilotée thermiquement entre les états liquides. Ceci est marqué par un pic « lambda » asymétrique dans la capacité calorifique et une augmentation brusque de la susceptibilité magnétique, reflétant probablement des fluctuations critiques proches d'une transition de phase quantique à température nulle.

Signification et affirmations
Cette publication démontre que les techniques de réseaux de tenseurs, spécifiquement les états de réseaux de tenseurs infinis optimisés avec la propagation de croyance et les corrections de boucles, constituent une voie puissante et scalable pour comprendre les aimants quantiques frustrés à des températures finies.

• Validation de la méthode : Ce travail valide l'approche de réseau de tenseurs basée sur la BP, capable de capturer avec précision les propriétés thermodynamiques sur de nombreux ordres de grandeur de température, incluant la détection des gaps d'énergie et des entropies résiduelles.
• Nature de l'état fondamental : L'étude clarifie que les plateaux magnétiques sur le réseau de rubis ne sont pas de simples ordres cristallins mais plutôt des « états de liquide de simplexes » — des mélanges désordonnés d'un nombre exponentiellement grand de configurations cristallines qui conservent une entropie résiduelle non nulle et ne brisent pas les symétries du réseau.
• Absence de transitions de phase : Un résultat clé est que le système atteint ces phases magnétiques stables sans subir de transitions de phase à température finie, comme en témoigne la capacité calorifique continue et finie.
• Directions futures : Les auteurs suggèrent que leur méthodologie peut être étendue pour étudier des systèmes frustrés en trois dimensions, tels que le réseau pyrochlore, et pour investiguer les perturbations qui pourraient briser la dégénérescence de l'état fondamental en faveur de superpositions à ordre topologique.