Initialement conçue pour traiter les fortes corrélations électroniques dans des hamiltoniens modèles, la théorie du champ moyen dynamique (DMFT) est devenue, en combinaison avec la théorie de la densité fonctionnelle (DFT), un outil puissant pour la simulation de matériaux réels. Au cœur de la DMFT se trouve la solution d'un modèle d'impureté quantique. Seulement les solutionneurs d'impureté Monte Carlo en temps continu (CT-QMC) donnent des solutions exactes. En même temps, ces solutionneurs sont plutôt complexes et gourmands en temps de calcul. Le but de cette thèse est donc d'améliorer les solutionneurs d'impureté CT-QMC. Après une courte introduction à la DMFT et à ses extensions pour les amas, on commence par une revue de deux des solutionneurs CT-QMC, celui en développement d'interaction ou de "Rubtsov" (CT-INT) et celui en développement d'hybridation (CT-HYB). Mettant l'accent sur le dernier, qui est l'algorithme de choix dans le cadre des matériaux réels, on montre alors comment réduire un problème de signe, nous permettant ainsi de traiter la criticalité inhabituelle des organiques en couche. Avec la supraconductivité à haute température critique comme exemple, on discute ensuite comment assurer l'ergodicité du solutionneur CT-HYB dans le cadre des symétries brisées. Finalement, des optimisations algorithmiques sont présentées et combinées, amenant à des accélérations allant jusqu'à un facteur de 500 dans le contexte des matériaux réels.
Initialement conçue pour traiter les fortes corrélations électroniques dans des hamiltoniens modèles, la théorie du champ moyen dynamique (DMFT) est devenue, en combinaison avec la théorie de la densité fonctionnelle (DFT), un outil puissant pour la simulation de matériaux réels. Au cœur de la DMFT se trouve la solution d'un modèle d'impureté quantique. Seulement les solutionneurs d'impureté Monte Carlo en temps continu (CT-QMC) donnent des solutions exactes. En même temps, ces solutionneurs sont plutôt complexes et gourmands en temps de calcul. Le but de cette thèse est donc d'améliorer les solutionneurs d'impureté CT-QMC. Après une courte introduction à la DMFT et à ses extensions pour les amas, on commence par une revue de deux des solutionneurs CT-QMC, celui en développement d'interaction ou de "Rubtsov" (CT-INT) et celui en développement d'hybridation (CT-HYB). Mettant l'accent sur le dernier, qui est l'algorithme de choix dans le cadre des matériaux réels, on montre alors comment réduire un problème de signe, nous permettant ainsi de traiter la criticalité inhabituelle des organiques en couche. Avec la supraconductivité à haute température critique comme exemple, on discute ensuite comment assurer l'ergodicité du solutionneur CT-HYB dans le cadre des symétries brisées. Finalement, des optimisations algorithmiques sont présentées et combinées, amenant à des accélérations allant jusqu'à un facteur de 500 dans le contexte des matériaux réels.
Cette thèse étudie l'apport du calcul quantique bruité pour l’algorithme phare des fortes corrélations, la théorie du champ moyen dynamique (DMFT). Elle vise à mettre à profit les premiers dispositifs de calcul quantique, malgré leurs imperfections liées à un degré de contrôle expérimental encore limité. Dans un premier temps, une version améliorée de la méthode variationnelle de préparation de l'état fondamental du modèle d'impureté est proposée. Elle consiste en la réalisation de mises à jour de la base à une particule dans laquelle est décrit le Hamiltonien d'impureté. Ces mises à jour sont entrelacées avec des optimisations variationnelles de l'état, et guidées par la matrice densité à une particule de l'état variationnel optimisé courant. Cet algorithme nous a permis de réaliser la première implémentation hybride bruitée d'un schéma assimilé à la DMFT avec un système auxiliaire à deux impuretés. Aussi, nous montrons sur plusieurs exemples que cette méthode est capable d'augmenter la capacité d'un circuit variationnel donné à représenter l'état cible. Enfin, nous proposons de combiner les mises à jour de la base à une particule avec un algorithme variationnel dit adaptatif, qui construit le circuit itérativement. Nous montrons que cette approche permet de réduire, à précision donnée sur l'énergie de l'état optimisé, le nombre de portes du circuit. Dans un second temps, nous proposons de mettre à profit la dissipation qui affecte les qubits afin de diminuer les effets de la troncation du bain sur l'ajustement de l'hybridation du bain à celle de la DMFT. Nous montrons qu'une réduction en termes de sites de bain est bien à la portée d'une telle méthode. Cependant, nous faisons l'hypothèse d'un processus dissipatif qui n'est pas réaliste : la méthode doit donc encore être étudiée via un modèle plus proche des conditions expérimentales.
Cette thèse concerne le développement et l’utilisation de nouveaux algorithmes pour l’étude théorique de systèmes quantiques fortement corrélés. Il s’agit de systèmes dans lesquels les interactions, par exemple l’interaction de Coulomb répulsive entre les électrons d’un matériau, induisent des phénomènes physiques remarquables (comme la supraconductivité à haute température critique) qui ne peuvent être décrits dans le cadre d’une théorie standard à un corps. Pour comprendre l’origine de ces phénomènes, il est indispensable de traiter l’ensemble des particules et leurs interactions: on parle du problème à N -corps.Le projet de thèse consiste à continuer le développement, l’analyse et l’application des méthodes numériques dites diagrammatiques à ces systèmes. Si des progrès ont été faits, il reste beaucoup de questions fondamentales ouvertes quant à l’utilisation de méthodes perturbatives pour décrire un système qui est,par nature, dans un régime non-perturbatif. Quelles sont les limites de ces approches? Comment se manifestent les effets des corrélations dans la structure des séries perturbatives?Les développements algorithmiques seront utilisés pour l’étude de modèles de systèmes fortement corrélés, comme le modèle de Hubbard, ce qui permettra d’aborder des questions physiques d’intérêt actuel, en particulier en relation avec la physique des corrélations magnétiques et du pseudo-gap dans les cuprates supraconducteurs ou avec l’existence d’une transition de Mott sans phase ordonnée préemptive telle qu’elle est observée dans des expériences récentes sur des composés organiques.
Algorithmes hybrides quantiques-classiques pour le problème à N corps quantique Cette thèse étudie l'apport du calcul quantique bruité pour l’algorithme phare des fortes corrélations, la théorie du champ moyen dynamique (DMFT). Elle vise à mettre à profit les premiers dispositifs de calcul quantique, malgré leurs imperfections liées à un degré de contrôle expérimental encore limité. Dans un premier temps, une version améliorée de la méthode variationnelle de préparation de l'état fondamental du modèle d'impureté est proposée. Elle consiste en la réalisation de mises à jour de la base à une particule dans laquelle est décrit le Hamiltonien d'impureté. Ces mises à jour sont entrelacées avec des optimisations variationnelles de l'état, et guidées par la matrice densité à une particule de l'état variationnel optimisé courant. Cet algorithme nous a permis de réaliser la première implémentation hybride bruitée d'un schéma assimilé à la DMFT avec un système auxiliaire à deux impuretés. Aussi, nous montrons sur plusieurs exemples que cette méthode est capable d'augmenter la capacité d'un circuit variationnel donné à représenter l'état cible. Enfin, nous proposons de combiner les mises à jour de la base à une particule avec un algorithme variationnel dit adaptatif, qui construit le circuit itérativement. Nous montrons que cette approche permet de réduire, à précision donnée sur l'énergie de l'état optimisé, le nombre de portes du circuit. Dans un second temps, nous proposons de mettre à profit la dissipation qui affecte les qubits afin de diminuer les effets de la troncation du bain sur l'ajustement de l'hybridation du bain à celle de la DMFT. Nous montrons qu'une réduction en termes de sites de bain est bien à la portée d'une telle méthode. Cependant, nous faisons l'hypothèse d'un processus dissipatif qui n'est pas réaliste : la méthode doit donc encore être étudiée via un modèle plus proche des conditions expérimentales.
La compréhension de la supraconductivité à haute température critique exhibée par certains oxydes de métaux de transition demeure une question centrale en physique théorique de la matière condensée. Dans ce contexte et depuis la proposition historique de P. W. Anderson, le modèle de Hubbard répulsif en géométrie bidimensionnelle est devenu un paradigme pour tenter de capturer les propriétés essentielles des matériaux supraconducteurs non conventionnels. Cependant, la détermination de son état fondamental exact se heurte à la complexité exponentielle du problème quantique à N-corps. Dans cette thèse, nous poursuivons l’objectif de développer un schéma variationnel d’approximation s’affranchissantde toute hypothèse concernant les ordres magnétiques, de charge ou supraconducteurs susceptibles d’émerger de l’Hamiltonien à basse énergie. L’approche trouve son originalité dans l’introduction de corrélations par la restauration, avant variation, de symétries délibérément brisées dans un état d’essai en forme d’une superposition de fonctions d’onde versatiles de type Hartree-Fock et Bogoliubov-de Gennes. Pour des petits amas à deux etquatre sites, nous montrons analytiquement que cette méthode de champ moyen enchevêtré par les symétries de l’Hamiltonien permet de retrouver l’état fondamental exact quelle que soit l’intensité des interactions. Pour de plus grands réseaux dopés en trous et dans les régimes fortement corrélés, nous mettons en exergue un arrangement des moments magnétiques en spirale ou suivant une onde de densité de spin qui est alors accompagnée d’inhomogénéités sous forme de rayures périodiquement réparties. De tels ordres sont de plus entrelacés avec des corrélations d’appariement en onde d à longue distance qui, à la limite thermodynamique, signent la supraconductivité. Ces résultats sont obtenus grâce à des simulations systématiques dans une géométrie de tubes de plaquettes pouvant être réalisée expérimentalement à l’aide d’atomes froids piégés dans des réseaux optiques.
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).
Le domaine des problèmes quantiques à N-corps à l'équilibre et hors d'équilibre sont des sujets majeurs de la Physique et de la Physique de la matière condensée en particulier. Les propriétés d'équilibre de nombreux systèmes unidimensionnels en interaction sont bien comprises d'un point de vue théorique, des chaînes de spins aux théories quantiques des champs dans le continue. Ces progrès ont été rendus possibles par le développement de nombreuses techniques puissantes, comme, par exemple, l'ansatz de Bethe, le groupe de renormalisation, la bosonisation, les états produits de matrices ou la théorie des champs invariante conforme. Même si les propriétés à l'équilibre de nombreux modèles soient connues, ceci n'est en général pas suffisant pour décrire leurs comportements hors d'équilibre, et ces derniers restent moins explorés et beaucoup moins bien compris. Les modèles d'impuretés quantiques représentent certains des modèles à N-corps les plus simples. Mais malgré leur apparente simplicité ils peuvent capturer plusieurs phénomènes expérimentaux importants, de l'effet Kondo dans les métaux aux propriétés de transports dans les nanostructures, comme les points quantiques. Dans ce travail nous considérons un modèle d'impureté appelé "modèle de niveau résonnant en interaction" (IRLM). Ce modèle décrit des fermions sans spin se propageant dans deux fils semi-infinis qui sont couplés à un niveau résonant -- appelé point ou impureté quantique -- via un terme de saut et une répulsion Coulombienne. Nous nous intéressons aux situations hors d'équilibre où un courant de particules s'écoule à travers le point quantique, et étudions les propriétés de transport telles que le courant stationnaire (en fonction du voltage), la conductance différentielle, le courant réfléchi, le bruit du courant ou encore l'entropie d'intrication. Nous réalisons des simulations numériques de la dynamique du modèle avec la méthode du groupe de renormalisation de la matrice densité dépendent du temps (tDMRG), qui est basée sur une description des fonctions d'onde en terme d'états produits de matrices. Nous obtenons des résultats de grande précision concernant les courbes courant-voltage ou bruit-voltage de l'IRLM, dans un grand domaine de paramètres du modèle (voltage, force de l'interaction, amplitude de saut vers le dot, etc.). Ces résultats numériques sont analysés à la lumière de résultats exacts de théorie des champs hors d'équilibre qui ont été obtenus pour un modèle similaire à l'IRLM, le modèle de Sine-Gordon avec bord (BSG). Cette analyse est en particulier basée sur l'identification d'une échelle d'énergie Kondo et d'exposants décrivant les régimes de petit et grand voltage. Aux deux points particuliers où les modèles sont connus comme étant équivalents, nos résultats sont en accord parfait avec la solution exacte. En dehors de ces deux points particuliers nous trouvons que les courbes de transport de l'IRLM et du modèle BSG demeurent très proches, ce qui était inattendu et qui reste dans une certaine mesure inexpliqué.