La recherche de phases magnétiques exotiques de la matière qui fondent même à T=0 uniquement sous l'action des fluctuations quantiques a été long et ardu, à la fois théoriquement et expérimentalement. La percée est venue récemment avec la découverte de l'Herbertsmithite, un composé formant un réseau kagome parfait avec des moments magnétiques de spin-1/2. Des expériences pionnières, mêlant des mesures de NMR, MuSR et de diffusion de neutrons, ont montré une absence totale de gel ou d'ordre des moments magnétiques de spin, fournissant ainsi une forte signature d'une phase paramgnétique quantique. Théoriquement, l'Herbertsmithite est extrêmement bien modélisé par le modèle de Heisenberg quantique antiferromagnétique pour des spins-1/2 sur le réseau kagome, problème qui n'a pas été résolu jusqu'à présent. Plusieurs méthodes approximatives numériques et analytiques ont donné différents états fondamentaux, allant des liquides de spins Z2 gappés et un liquide de spins exotique algébrique U(1) de Dirac aux liquides de spins chiraux et les cristaux à liaisons de valence. Dans cette thèse, le problème est traité dans le cadre d'une approche particule-esclave fermionique, à savoir le formalisme des fermions de Schwinger SU(2). Il est conclu qu'un liquide de spins sans gap algébrique de Dirac a l'énergie variationnelle la plus basse et peut en fait constituer un vrai état fondamental physique de liquide de spins. Une implémentation sophistiquée de méthodes numériques de pointes comme le Monte-Carlo variationnel, le Monte-Carlo fonctions de Green et l'application de pas Lanczos dans un schéma variationnel ont été utilisés. Il est montré que contrairement à la croyance habituelle, le liquide de spins de Dirac U(1) projeté en "2+1" dimensions est remarquablement robuste par rapport à une large classe de perturbations, incluant les liquides de spins topologiques Z2 et les cristaux à liaisons de valence. De plus, l'application de deux pas Lanczos sur la fonction d'onde du liquide de spins de Dirac U(1) montre que son énergie est compétitive avec celles proposées pour les liquides de spins topologiques Z2. Ce résultat, combiné avec les indications expérimentales qui pointent vers un liquide de spins sans gap pour l'Herbertsmithite, appuie l'affirmation que le vrai état fondamental de ce modèle est en fait un liquide de spins algébrique de Dirac.
La recherche de phases magnétiques exotiques de la matière qui fondent même à T=0 uniquement sous l'action des fluctuations quantiques a été long et ardu, à la fois théoriquement et expérimentalement. La percée est venue récemment avec la découverte de l'Herbertsmithite, un composé formant un réseau kagome parfait avec des moments magnétiques de spin-1/2. Des expériences pionnières, mêlant des mesures de NMR, MuSR et de diffusion de neutrons, ont montré une absence totale de gel ou d'ordre des moments magnétiques de spin, fournissant ainsi une forte signature d'une phase paramgnétique quantique. Théoriquement, l'Herbertsmithite est extrêmement bien modélisé par le modèle de Heisenberg quantique antiferromagnétique pour des spins-1/2 sur le réseau kagome, problème qui n'a pas été résolu jusqu'à présent. Plusieurs méthodes approximatives numériques et analytiques ont donné différents états fondamentaux, allant des liquides de spins Z2 gappés et un liquide de spins exotique algébrique U(1) de Dirac aux liquides de spins chiraux et les cristaux à liaisons de valence. Dans cette thèse, le problème est traité dans le cadre d'une approche particule-esclave fermionique, à savoir le formalisme des fermions de Schwinger SU(2). Il est conclu qu'un liquide de spins sans gap algébrique de Dirac a l'énergie variationnelle la plus basse et peut en fait constituer un vrai état fondamental physique de liquide de spins. Une implémentation sophistiquée de méthodes numériques de pointes comme le Monte-Carlo variationnel, le Monte-Carlo fonctions de Green et l'application de pas Lanczos dans un schéma variationnel ont été utilisés. Il est montré que contrairement à la croyance habituelle, le liquide de spins de Dirac U(1) projeté en "2+1" dimensions est remarquablement robuste par rapport à une large classe de perturbations, incluant les liquides de spins topologiques Z2 et les cristaux à liaisons de valence. De plus, l'application de deux pas Lanczos sur la fonction d'onde du liquide de spins de Dirac U(1) montre que son énergie est compétitive avec celles proposées pour les liquides de spins topologiques Z2. Ce résultat, combiné avec les indications expérimentales qui pointent vers un liquide de spins sans gap pour l'Herbertsmithite, appuie l'affirmation que le vrai état fondamental de ce modèle est en fait un liquide de spins algébrique de Dirac.
Etudes Monte Carlo de liquides de spins de Dirac et liquides de spins chiraux Les liquides de spins quantiques (QSL) sont des états exotiques de la matière caractérisés par leur intrication quantique à N corps et leurs excitations fractionnaires émergentes, dont la physique des basses énergies est décrite à l'aide des théories de jauge. La réalisation expérimentale de QSL dans des matériaux cristallins et leur existence en tant qu'états fondamentaux stables dans des modèles de spins ont été des questions fondamentales et largement étudiées en physique de la matière condensée ces dernières années. Dans cette thèse, nous utilisons des simulations quantiques de Monte Carlo pour explorer deux questions sur les QSL en deux dimensions, l'une sur les liquides de spins avec et sans bande interdite. La première question concerne la question de savoir si les états de paires intriqués projetés (PEPS, un type de réseau tensoriel) peuvent représenter avec précision les liquides de spins chiraux (CSL). Des questions ont été soulevées en raison d'un théorème interdisant les représentations PEPS à fermions libres d'états chiraux sans interaction. Nous construisons des approximants PEPS fermioniques (fPEPS) des isolants Chern projetés par Gutzwiller (un type de CSL) et démontrons que les fPEPS (de dimension de liaison finie) peuvent capturer la dégénérescence topologique correcte de l'état fondamental du CSL. De plus, des fPEPS plus généraux sont optimisés pour décrire la phase CSL d'un antiferromagnétique Heisenberg frustré, et il est démontré que les modes de bord chiraux dans le spectre d'intrication suivent les prédictions de la théorie conforme des champs. Par conséquent, notre travail fournit des preuves supplémentaires démontrant que le PEPS peut représenter efficacement les phases CSL. La deuxième question concerne la nature des excitations à faible énergie des liquides de spin de Dirac sans bande interdite. Contrairement aux liquides de spin espacés, tels que certains liquides de spin Z2, dont les excitations à faible énergie (par exemple, les spinons et les visons) peuvent être comprises qualitativement grâce à la théorie du champ moyen, les liquides de spin de Dirac présentent des défis plus subtils. Il reste à déterminer si les descriptions en champ moyen des DSL peuvent résister aux effets des fluctuations de jauge. De plus, la stabilité des DSL en tant qu'états fondamentaux dans les antiferromagnétiques de Heisenberg frustrés reste une question ouverte. Dans ce travail, nous construisons des analyses variationnelles pour les excitations à faible énergie des DSL. Nos résultats indiquent que les modèles Heisenberg sur les réseaux triangulaires et Kagome prennent en charge des excitations monopôles et spinons sans interruption au-delà du niveau de champ moyen (après application de la projection de Gutzwiller). Nous démontrons en outre que le liquide de spin de Dirac est énergétiquement stable contre la formation d'états liquides de spin chiraux générés par de multiples insertions de monopôles. De plus, les impulsions monopôles calculées sont cohérentes avec les récentes prédictions de la théorie des champs. Enfin, grâce à une diagonalisation exacte, nous constatons que les fonctions d'onde monopôle et spinon projetées par Gutzwiller présentent un recouvrement significatif dans les amas finis, ce qui suggère que ces quasi-particules pourraient ne pas être entièrement indépendantes dans les petits systèmes.
Le liquide de luttinger constitue le meilleur exemple de metal non-liquide de fermi connu. Il designe une classe d'universalite pour des systemes unidimensionnels caracterises par une dynamique purement collective sans quasi-particules de landau. Des realisations experimentales du liquide de luttinger existent, par exemple sur les bords d'un gaz d'electrons bidimensionnel sous champ magnetique tel que le fluide de hall fractionnaire. De nombreux modeles integrables exactement solubles font aussi partie de cette classe : la chaine d'heisenberg, le modele de hubbard de facon tres remarquable, les solutions exactes font apparaitre dans le spectre d'excitations des etats fractionnaires, qui ne portent plus les nombres quantiques de l'electron : par exemple le holon du modele de hubbard porte une charge unite mais n'a pas de spin. Cependant, aucune analyse detaillee de la theorie effective decrivant le liquide de luttinger qui ferait clairement apparaitre des etats fractionnaires dans le spectre n'existait. Dans cette these, nous avons montre explicitement que cette theorie effective admettait des excitations fractionnaires et avons decrit en detail leur construction. Nous avons pu ainsi obtenir toute une serie de resultats nouveaux sur les excitations fractionnaires, ce qui nous a conduit a faire un certain nombre de predictions. Nous avons trouve de facon tres remarquable qu'il n'y avait non plus une mais deux excitations elementaires correspondant au spinon des chaines de spin et une seconde qui est isomorphe a une quasi-particule de laughlin de l'effet hall fractionnaire. Cette derniere serait susceptible d'etre observee experimentalement. Les fonctions d'onde du fondamental et des excitations ont aussi ete calculees. Ces resultats ont ete ensuite generalises au liquide de luttinger avec spin : nous introduisons un formalisme qui nous permet de decrire de facon exhaustive le spectre fractionnaire dans des situations avec separation spin-charge ou meme sans separation spin-charge (comme par exemple pour des systemes sous champ magnetique comme le modele de hubbard sous champ).
La frustration magnétique est une propriété physique faisant l'objet de très nombreuses études. Elle intervient lorsque la géométrie du système ne permet pas la minimisation simultanée des énergies locales, menant alors à une dégénérescence macroscopique. Présent dans la nature au travers de différents minéraux, dont la célèbre herbertsmithite, le réseau kagomé antiferromagnétique est un parfait exemple de système bidimensionnel frustré. Lorsqu'à cette propriété s'ajoute un champ magnétique extérieur tendant à polariser le matériau, la compétition de ces deux phénomènes donne naissance à des processus d'aimantation exotiques, tels que les plateaux d'aimantation pour lesquels la susceptibilité magnétique est nulle. Nous avons étudié le réseau kagomé avec le hamiltonien AKLT de spin 2 et le hamiltonien de Heisenberg pour les spins 1/2, 1, 3/2 et 2 sous champ magnétique. Ce travail a été effectué en utilisant des méthodes numériques appelées méthodes Tensor Network. Tout d'abord développées dans le cadre de la théorie quantique de l'information, ces méthodes, basées sur le principe d'intrication, ont très vite été utilisées en matière condensée. Elles ont permis l'étude de systèmes bidimensionnels frustrés, jusqu'alors limitée aux méthodes numériques déjà connus telles que les méthodes Monte Carlo ou encore DMRG. Nous avons été en mesure de produire les diagrammes de phases, d'une grande richesse, pour tous les modèles cités, en fonction du champ magnétique.
Cette thèse porte sur trois modèles d’électrons fortement corrélés : le modèle de Kitaev-Hubbard, le modèle de Hubbard sur le réseau en treillis et le modèle de Hubbard sur le réseau en nid d’abeille. La résolution de ces trois modèles est basée sur les méthodes d’amas quantiques telles que l’approximation de l’amas variationnel (VCA), le champ moyen dynamique sur amas (CDMFT) et l’approximation de l’impureté dynamique sur amas (CDIA). Le solutionneur utilisé avec ces méthodes est la diagonalisation exacte (ED) à température nulle. Le modèle de Kitaev-Hubbard est une variante du célèbre modèle de Hubbard, défini sur le réseau du graphène et comportant un terme de saut t′ dépendant du spin, qui joue un rôle analogue à celui d’une interaction spin-orbite. Au demi-remplissage et dans la limite à couplage fort (U → ∞), ce modèle est une interpolation entre le modèle de Heisenberg pour l’antiferromagnétisme et le modèle de Kitaev utilisé en information quantique. Contrairement au modèle idéalisé de Kitaev, le modèle de Kitaev-Hubbard peut être réalisé dans des réseaux optiques d’atomes froids et ainsi supporter, dans certaines limites (t′ ∼ t, U ≫ t), des calculs quantiques topologiques. À couplage modéré, notre étude conduit à un diagramme de phase comportant de l’antiferromagnétisme, de liquide de spin algébrique ainsi que deux types de semi-métaux. Ces différentes phases du modèle se rencontrent en un seul point sur le plan (U, t′). La transition entre la phase antiferromagnétique et les phases semi-métallique est du deuxième ordre, alors qu’elle est du premier ordre entre le liquide de spin algébrique et l’antiferromagnétisme. La méthode d’amas quantique utilisée dans ce projet est l’approximation de l’impureté dynamique sur amas (CDIA). Le réseau en treillis constitue une description simplifiée du matériau Sr[indice inférieur 14−x]Ca[indice inférieur x]Cu[indice inférieur 24]O[indice inférieur 41], qui a une structure en échelles couplées. C’est un matériau supraconducteur découvert en 1996, mais dont la température critique est très sensible à la pression. La présence et le mécanisme de la supraconductivité, ainsi que la symétrie du paramètre d’ordre restent des sujets de recherche. Le matériau présente aussi des propriétés similaires aux cuprates et reste le seul supraconducteur à base d’oxyde de cuivre n’ayant pas une structure plane. Pour cette raison, l’étude du Sr[indice inférieur 14−x]Ca[indice inférieur x]Cu[indice inférieur 24]O[indice inférieur 41] est d’intérêt dans le but de comprendre le mécanisme de la supraconductivité causée par les corrélations électroniques, comme celle des cuprates. Notre étude montre que la supraconductivité est présente dans une gamme de densités électroniques s’échelonnant entre le demi-remplissage et un dopage d’environ 20%. De plus, les résultats montrent que le paramètre d’ordre est de type singulet d + id dans une gamme de dopage. Ainsi, la supraconductivité est peut-être chirale, c’est-à-dire qu’elle brise l’invariance par inversion du temps. Ce système porterait donc du courant spontané à la périphérie de l’échantillon même en absence de champ magnétique. Pour étudier ce modèle, nous avons utilisé l’approximation de l’amas variationnel (VCA) et le champ moyen dynamique sur amas (CDMFT). Le modèle de Hubbard étendu sur le réseau en nid d’abeille décrit approximativement le graphène et les matériaux analogues. L’étude de ce modèle montre que la supraconductivité peut être présente dans ce système. De plus, en fonction du dopage, notre étude montre que le paramètre d’ordre est de type triplet et de symétrie p + ip. Ainsi, la supraconductivité serait chirale. Les résultats indiquent aussi que la répulsion électron-électron entre atomes voisins favorise la supraconductivité dans sa forme chirale. Ces résultats ont été obtenus à la fois par l’approximation de l’amas variationnel (VCA) et le champ moyen dynamique sur amas (CDMFT).
Méthodes de réseaux tensoriels pour la dynamique hors équilibre des liquides de spin quantique La présence d'interactions fortes et d'intrication dans un système quantique à plusieurs corps peut conduire à plusieurs états exotiques de la matière. Les liquides Quantum Spin sont l'un de ces ensembles d'états exotiques qui ont suscité un intérêt considérable en raison de leurs propriétés intéressantes et de leurs applications technologiques possibles. Les progrès rapides des montages d'expériences sur les atomes froids, en particulier les simulateurs d'atomes de Rydberg, ont ouvert de nouvelles voies pour étudier ces liquides de spin en une et deux dimensions. Cela nécessite des techniques théoriques avancées pour étudier les propriétés de non-équilibre d'états quantiques purs isolés. Cette thèse vise à développer des algorithmes numériques pour étudier les propriétés de non-équilibre des liquides de spin sur un réseau carré suite à une trempe quantique en utilisant le cadre des réseaux de tenseurs. Une forme variationnelle de l'ansatz Projected Entangled Pair State (PEPS) est utilisée pour décrire les liquides de spin. Une nouvelle forme symétrique du groupe de renormalisation de la matrice de transfert de coin (CTMRG), l'algorithme nécessaire pour construire la matrice de densité réduite de la fonction d'onde dans la limite thermodynamique pour une fonction d'onde PEPS à valeurs complexes, est décrite. Pour étudier la dynamique de trempe, deux algorithmes appelés update simple et update complète, ainsi que les modifications nécessaires pour les adapter à l'étude des fonctions d'onde de spin liquide, sont expliqués. Les outils sont ensuite utilisés pour étudier la dynamique de non-équilibre de l'état de la liaison de valence résonnante (RVB) avec les configurations de liaison de valence du plus proche voisin (NN), un état spin-liquide connu pour présenter des corrélations dimère-dimère critiques, dans le Collecteur PEPS.
D ans cette présente thèse de doctorat et dans le but de découvrir le comportement magnétique des matériaux en fonction de di érents paramètres extérieurs, nous avons proposé quelques modèles de spin sur réseaux pour l’étude. Pour ce fait, nous avons donné en détails, dans le premier chapitre, l’origine du magnétisme dans les di érents matériaux et leur comportement en fonction du changement des paramètres. Dans le second chapitre, nous avons mis l’accent sur un phénomène extraordinaire qui se manifeste au niveau des matériaux qui est la transition de phase. Cette dernière n’est qu’un changement d’état subi à cause des modifications au niveau des paramètres. Nous avons cité les di érentes classifications des transitions de phase notamment la théorie d’Ehrenfest et celle de Landau. Nous avons donné un peu d’intérêt à la deuxième classification avec plus de détails. Nous avons présenté aussi les phénomènes tricritiques qui mènent au calcul des di érents exposants critiques. Nous avons expliqué aussi la notion des classes d’universalité auxquelles appartiennent les di érents systèmes étudiés. Dans le troisième chapitre, nous avons présenté notre travail de fond. Nous avons défini le modèle à étudier à savoir le modèle de Blume-Emery-Gri ths. nous avons donné un aperçu sur les méthodes qu’on a utilisé à citer la méthode de Monte-Carlo et la méthode (finite size scaling) TMFSS basée sur la matrice de transfert. Ces méthodes ont montré leurs preuves dans plusieurs études antérieures. Le diagramme de phase obtenu, riche avec ses caractéristiques, exhibe une ligne de transition de second ordre. Nous avons justifié la présence des di érentes phases par étude numérique de l’évolution thermique de l’aimantation du système, du moment quadrupolaire et du moment octupolaire. Les résultats sont en accord avec ceux des états fondamentaux. En utilisant la méthode de la matrice de transfert, nous avons étudié le comportement critique du système. Nous avons calculé l’exposant critique pour quelques valeurs des constantes d’interaction. Le quatrième chapitre a été consacré à l’étude des états fondamentaux de deux modèles importants. Le premier est le modèle de Blume-Emery-Gri ths à spin-1. Les interactions sont antiferromagnétiques avec le premier plus proche voisin. Nous avons considéré un système triangulaire avec sa particularité de frustration. Nous avons fait une étude analytique qui nous a donné les énergies des di érentes configurations. Une étude numérique nous a montré les phases les plus stables qu’on a représenté sur un diagramme de phase. Le deuxième modèle était le modele de Blume-Capel à spin-3/2 sur un réseau carré. Nous avons suivi les mêmes étapes d’étude. Comme perspective, il serait intéressant de compléter les deux derniers travaux par des études numériques. Nous prévoyons utiliser les deux méthodes exactes à savoir la méthode de Monte-Carlo et celle de la matrice de transfert. Nous essayerons de calculer les di érents exposants critiques et de déterminer leur classe d’universalité.