Le modèle en couches constitue aujourd'hui un cadre théorique de référence pour appréhender les propriétés des noyaux à basse énergie. Son applicabilité demeure toutefois limitée par une croissance rédhibitoire de la dimension de l'espace des états aussi bien avec le nombre de couches de valence qu'avec le nombre de nucléons. Les méthodes Monte-Carlo quantiques (QMC) permettent a priori de contourner une telle difficulté en offrant une alternative à la diagonalisation du hamiltonien. Elles reposent sur une reformulation stochastique de l'équation de Schrödinger qui ramène le problème à N-corps à un ensemble de problèmes à un corps, numériquement solubles, et décrivant des particules indépendantes évoluant chacune dans un champ extérieur fluctuant. L'originalité de l'échantillonnage QMC proposé dans ce mémoire réside dans l'utilisation d'une approche variationnelle, avec restauration des symétries avant variation, pour guider le mouvement brownien et pour le contraindre afin de contrôler le problème de phase inhérent aux schémas QMC pour des fermions en interaction. Dans les couches sd et fp et avec les interactions résiduelles réalistes, nous avons ainsi pu obtenir une spectroscopie " yrast " en excellent accord avec les résultats issus de la diagonalisation du hamiltonien. En outre, une ouverture vers les systèmes d'électrons fortement corrélés est présentée au travers de nouveaux schémas QMC récemment suggérés pour le modèle de Hubbard en géométrie bidimensionnelle. Contrairement aux échantillonnages traditionnels, ils garantissent des trajectoires à poids positifs quel que soit le régime considéré. Nous avons prouvé que ces méthodes sont en réalité reliées à l'approche stochastique mise en oeuvre pour le modèle en couches. L'origine des erreurs systématiques qu'exhibent ces schémas, pourtant exempts de problème de signe avec le hamiltonien de Hubbard, est par ailleurs élucidée.
Le modèle en couches constitue aujourd'hui un cadre théorique de référence pour appréhender les propriétés des noyaux à basse énergie. Son applicabilité demeure toutefois limitée par une croissance rédhibitoire de la dimension de l'espace des états aussi bien avec le nombre de couches de valence qu'avec le nombre de nucléons. Les méthodes Monte-Carlo quantiques (QMC) permettent a priori de contourner une telle difficulté en offrant une alternative à la diagonalisation du hamiltonien. Elles reposent sur une reformulation stochastique de l'équation de Schrödinger qui ramène le problème à N-corps à un ensemble de problèmes à un corps, numériquement solubles, et décrivant des particules indépendantes évoluant chacune dans un champ extérieur fluctuant. L'originalité de l'échantillonnage QMC proposé dans ce mémoire réside dans l'utilisation d'une approche variationnelle, avec restauration des symétries avant variation, pour guider le mouvement brownien et pour le contraindre afin de contrôler le problème de phase inhérent aux schémas QMC pour des fermions en interaction. Dans les couches sd et fp et avec les interactions résiduelles réalistes, nous avons ainsi pu obtenir une spectroscopie " yrast " en excellent accord avec les résultats issus de la diagonalisation du hamiltonien. En outre, une ouverture vers les systèmes d'électrons fortement corrélés est présentée au travers de nouveaux schémas QMC récemment suggérés pour le modèle de Hubbard en géométrie bidimensionnelle. Contrairement aux échantillonnages traditionnels, ils garantissent des trajectoires à poids positifs quel que soit le régime considéré. Nous avons prouvé que ces méthodes sont en réalité reliées à l'approche stochastique mise en oeuvre pour le modèle en couches. L'origine des erreurs systématiques qu'exhibent ces schémas, pourtant exempts de problème de signe avec le hamiltonien de Hubbard, est par ailleurs élucidée.
Simulation des noyaux et de leurs constituants avec des ordinateurs quantiques La méthode de Brisure de Symétrie/Restauration de Symétrie est un outil largement utilisé pour traiter un ensemble de particules en interaction. Cette méthode permet souvent d'améliorer l'approximation de l'énergie de l'état fondamental d'un Hamiltonien dans un cadre variationnel. Dans cette technique, certaines symétries du système sont brisées intentionnellement dans la fonction d'onde utilisée lors de la recherche de l'état fondamental. Ce processus permet de capturer des corrélations complexes entre les particules à un coût numérique plus faible par rapport à un cadre qui préserve ces symétries. Étant donné que le véritable état fondamental doit respecter les symétries du problème, une étape de restauration de symétrie est a priori nécessaire pour revenir dans le sous-espace approprié de l'espace de Hilbert total.Cette thèse examine comment mettre en œuvre le schéma de Brisure de Symétrie/Restauration de Symétrie à l'aide d'ordinateurs quantiques. L'algorithme quantique que nous utilisons pour la composante variationnelle est la méthode de ``Variational Quantum Eigensolver (VQE)". Son but est de préparer une approximation de l'état fondamental en utilisant une fonction d'onde paramétrisée afin de minimiser l'énergie associée au Hamiltonien. L'ansatz quantique de Bardeen-Cooper-Schrieffer (BCS) a été utilisé comme ansatz brisant les symétries, et les techniques développées ont été testées en utilisant les Hamiltoniens d'appariement et le modèle de Hubbard. Nous avons également utilisé la transformation de Jordan-Wigner pour encoder les opérateurs présents dans le Hamiltonien. Avec ces outils, à l'instar des calculs sur ordinateurs classiques, deux manières d'utiliser le processus de Brisure de Symétrie/Restauration de Symétrie en conjonction avec la méthode VQE sont développés : nous pouvons varier les paramètres d'un ansatz de Brisure de Symétrie avant ou après la restauration de symétrie. La première approche est identifiée comme Projection Quantique Après Variation, tandis que la dernière est nommée Variation Quantique Après Projection.L'un des aboutissements de cette thèse est le développement de diverses techniques pour la restauration de symétrie. Certaines d'entre elles sont basées sur l'algorithme d'Estimation de Phase Quantique (QPE), ce qui comprend les méthodes de QPE itérative et de Rodéo. Nous avons également étudié d'autres techniques basées sur le concept ``d'Oracles quantiques", telles que les méthodes Oracle+Hadamard et Grover-Hoyer. L'algorithme de Combinaison Linéaire d'opérateurs Unitaires a été utilisé pour implémenter directement soit la projection, soit l'oracle. Nous avons également adapté le formalisme de ``Classical Shadows" pour la restauration de symétries dans le but d'optimiser les ressources quantiques nécessaires. Les différentes méthodes de restauration de symétrie sont illustrées à travers la thèse pour les symétries associées au nombre de particules et au spin total.Dans la dernière partie de la thèse, des techniques dites hybrides quantiques-classiques sont présentées, qui peuvent aider à extraire des informations précieuses sur le spectre de basse énergie du Hamiltonien. En supposant qu'il est possible d'extraire avec précision les moments du Hamiltonien à partir de leur fonction génératrice à l'aide d'un ordinateur quantique, nous introduisons deux méthodes d'analyse des spectres : la méthode de t-expansion, qui combine l'évolution en temps imaginaire et l'approximation de Padé, et la méthode de Krylov, une méthode d'expansion du sous-espace quantique qui peut fournir des informations sur les valeurs propres du Hamiltonien et sur l'évolution de la probabilité de survie. De plus, nous présentons également des résultats de la méthode Krylov Quantique. Cette technique donne des informations similaires à la méthode de Krylov, mais sans avoir besoin d'estimer les moments du Hamiltonien, une tâche qui reste difficile sur un ordinateur quantique actuellement.
La méthode de Brisure de Symétrie/Restauration de Symétrie est un outil largement utilisé pour traiter un ensemble de particules en interaction. Cette méthode permet souvent d'améliorer l'approximation de l'énergie de l'état fondamental d'un Hamiltonien dans un cadre variationnel. Dans cette technique, certaines symétries du système sont brisées intentionnellement dans la fonction d'onde utilisée lors de la recherche de l'état fondamental. Ce processus permet de capturer des corrélations complexes entre les particules à un coût numérique plus faible par rapport à un cadre qui préserve ces symétries. Étant donné que le véritable état fondamental doit respecter les symétries du problème, une étape de restauration de symétrie est a priori nécessaire pour revenir dans le sous-espace approprié de l'espace de Hilbert total.Cette thèse examine comment mettre en œuvre le schéma de Brisure de Symétrie/Restauration de Symétrie à l'aide d'ordinateurs quantiques. L'algorithme quantique que nous utilisons pour la composante variationnelle est la méthode de ``Variational Quantum Eigensolver (VQE)". Son but est de préparer une approximation de l'état fondamental en utilisant une fonction d'onde paramétrisée afin de minimiser l'énergie associée au Hamiltonien. L'ansatz quantique de Bardeen-Cooper-Schrieffer (BCS) a été utilisé comme ansatz brisant les symétries, et les techniques développées ont été testées en utilisant les Hamiltoniens d'appariement et le modèle de Hubbard. Nous avons également utilisé la transformation de Jordan-Wigner pour encoder les opérateurs présents dans le Hamiltonien. Avec ces outils, à l'instar des calculs sur ordinateurs classiques, deux manières d'utiliser le processus de Brisure de Symétrie/Restauration de Symétrie en conjonction avec la méthode VQE sont développés : nous pouvons varier les paramètres d'un ansatz de Brisure de Symétrie avant ou après la restauration de symétrie. La première approche est identifiée comme Projection Quantique Après Variation, tandis que la dernière est nommée Variation Quantique Après Projection.L'un des aboutissements de cette thèse est le développement de diverses techniques pour la restauration de symétrie. Certaines d'entre elles sont basées sur l'algorithme d'Estimation de Phase Quantique (QPE), ce qui comprend les méthodes de QPE itérative et de Rodéo. Nous avons également étudié d'autres techniques basées sur le concept ``d'Oracles quantiques", telles que les méthodes Oracle+Hadamard et Grover-Hoyer. L'algorithme de Combinaison Linéaire d'opérateurs Unitaires a été utilisé pour implémenter directement soit la projection, soit l'oracle. Nous avons également adapté le formalisme de ``Classical Shadows" pour la restauration de symétries dans le but d'optimiser les ressources quantiques nécessaires. Les différentes méthodes de restauration de symétrie sont illustrées à travers la thèse pour les symétries associées au nombre de particules et au spin total.Dans la dernière partie de la thèse, des techniques dites hybrides quantiques-classiques sont présentées, qui peuvent aider à extraire des informations précieuses sur le spectre de basse énergie du Hamiltonien. En supposant qu'il est possible d'extraire avec précision les moments du Hamiltonien à partir de leur fonction génératrice à l'aide d'un ordinateur quantique, nous introduisons deux méthodes d'analyse des spectres : la méthode de t-expansion, qui combine l'évolution en temps imaginaire et l'approximation de Padé, et la méthode de Krylov, une méthode d'expansion du sous-espace quantique qui peut fournir des informations sur les valeurs propres du Hamiltonien et sur l'évolution de la probabilité de survie. De plus, nous présentons également des résultats de la méthode Krylov Quantique. Cette technique donne des informations similaires à la méthode de Krylov, mais sans avoir besoin d'estimer les moments du Hamiltonien, une tâche qui reste difficile sur un ordinateur quantique actuellement.
Les calculateurs quantiques promettent de résoudre des problèmes insolvables pour les ordinateurs classiques. En particulier, manipuler des états superposés et intriqués de bits quantiques est une voie naturelle pour étudier la physique des systèmes quantiques à N particules. Avec les avancées rapides des technologies quantiques, nous avons pu voir émerger les premiers processeurs quantiques, nommés processeurs "NISQ" (de l'anglais "Noisy-Intermediate-Scale-Quantum"). Les progrès de ces machines ouvrent de nouveaux chemins aux physiciens pour simuler des systèmes quantiques. Cependant, leur taille limitée et les erreurs liées à la décohérence limitent fortement leur usage aujourd'hui. La conception d'algorithmes quantiques efficaces, résistant aux erreurs, est un défi considérable. Ce travail de thèse s'attelle à comprendre les capacités des processeurs quantiques à simuler des systèmes quantiques à N corps, avec une attention particulière aux applications pour la physique de la matière condensée. Dans un premier temps, nous abordons le sujet du calcul de l'état fondamental d'Hamiltoniens, à l'aide d'un algorithme variationnel quantique-classique (nommée l'algorithme VQE, de l'anglais « Variational Quantum Eigensolver »). Nous évaluons les performances de cette méthode sur le modèle de Hubbard, grâce à différentes métriques, comme la fidélité, l'énergie et les propriétés physiques des solutions obtenues. Nous démontrons, malgré des performances prometteuses en simulant un calculateur quantique sans bruit, que le niveau de bruit des machines contemporaines rend impossible la réalisation d'un calcul significatif. En revanche, certains états quantiques aux propriétés intéressantes peuvent avoir une forme simple. Leur préparation sur un calculateur quantique est un problème pertinent dans le but d'initialiser un calcul quantique avec ces états pour les amener vers des régimes où ils perdent leur forme simple. Pour illustrer ceci, nous prenons l'exemple du modèle Affleck-Lieb-Kennedy-Tasaki (AKLT), le premier modèle de chaîne de spin entier antiferromagnétique possédant une solution simple vérifiant la conjecture d'Haldane. Nous proposons deux schémas pour préparer les états AKLT, l'un suivant leur construction naturelle en projetant des paires de qubits sur leur sous-espace triplet, l'autre utilisant leur forme de réseaux de tenseurs. Nous analysons l'implémentation de ces méthodes en fonction de la taille du système, avec une réalisation expérimentale sur un processeur quantique. Enfin, nous explorons la synergie entre les méthodes de réseaux de tenseurs et les calculateurs quantiques bruités pour simuler la dynamique hors-équilibre de systèmes de spins. La simulation de dynamique quantique est une tâche naturelle pour les calculateurs quantiques, mais atteindre des temps longs est fortement limité par le bruit. D'autre part, les réseaux de tenseurs sont également limités aux temps courts, contraints par la quantité d'intrication qu'ils peuvent encoder. Nous proposons d'utiliser le meilleur des deux mondes : la dynamique à temps court est simulée efficacement par un réseau de tenseur, compilée vers des circuits quantiques de faible profondeur pour poursuivre l'évolution temporelle sur le calculateur quantique. Nous étudions les performances de ce schéma hybride quantique-classique en regardant la fidélité des états obtenus et la production d'intrication, avec une validation expérimentale. Finalement, nous proposons d'optimiser les circuits quantiques approximant l'opérateur évolution temporelle grâce à des réseaux de tenseurs, résultant en une simulation quantique classiquement améliorée utilisant moins de ressources quantiques.
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.
Nouveaux algorithmes pour l’étude des propriétés d’équilibre et hors d’équilibre des systèmes quantiques fortement corrélés Quel est le point commun entre les étoiles formant une galaxie, les gouttes d'eau s'écoulant dans une rivière, et les électrons d'une céramique superconductrice lévitant au-dessus d'un aimant ? Tous ces systèmes ne peuvent être décrits par le mouvement isolé d'une seule de leurs composantes. C'est l'ensemble des particules et de leurs interactions qui fait émerger leurs singulières propriétés : on parle du problème à N corps.Dans cette Thèse, nous nous intéressons aux propriétés des systèmes d'électrons fortement corrélés, dont la physique est gouvernée par les principes de la mécanique quantique. Les méthodes analytiques étant rapidement limitées, nous développons de nouvelles approches numériques afin de quantifier précisément les propriétés de matériaux dans lesquels les interactions entre particules deviennent importantes.Nous nous intéressons tout d'abord aux propriétés d'équilibre de la pérovskite Sr2IrO4, un matériau structurellementéquivalent au cuprate supraconducteur La2CuO4. Nous mettons en évidence l'existence d'un pseudogapet décrivons la structure électronique de ce matériau en fonction du dopage.Nous développons ensuite des extensions aux algorithmes de Monte Carlo déterminantaux pour l'étude de quantités dynamiquescomme l'énergie propre, et nous montrons qu'il est possible de regrouper un nombre factoriel de diagrammes en une somme de déterminants, réduisantainsi fortement le problème de signe fermionique.Dans un deuxième temps, nous décrivons les systèmes fortement corrélés hors d'équilibre.Nous commençons par revisiter le Monte Carlo diagrammatique en temps réel dans une nouvelle base qui permet aux diagrammes du vide de s'annulerdirectement. Au cours d'un échantillonnage statistique, ceci permetd'atteindre la limite de long temps nécessaire à l'étude des états stationnaires des systèmes hors d'équilibre.Pour terminer, nous étudions la transition métal-isolant induite par un champ électrique de Ca2RuO4, qui coexiste avec une transition structurelle.Un algorithme basé sur l'approximation sans croisement nous permettent de calculer le courant en fonction du champ crystallin dans ce matériauet de comparer nos résultats aux données expérimentales.
Étude numérique des systèmes à N corps contraints Cette thèse présente une étude numérique de trois modèles de physique à N corps soumis à des contraintes directes ou indirectes. La présence de contraintes dans un système complexe à plusieurs corps requiert des algorithmes numériques non triviaux pour les étudier. Les deux premiers modèles comportent une forme directe de contrainte locale et sont le modèle de Dimères Quantiques (QDM) de Rokhsar-Kivelson et un modèle de mécanique statistique classique de boucles sans croisement avec des interactions attractives, tous les deux étudiés sur un réseau carré. L'étude de tels modèles contraints a connu un regain d'intérêt récemment avec leur réalisations directes sur des simulateurs quantiques de réseaux d'atomes de Rydberg. Les simulations du modèle classique utilisent l'algorithme Monte Carlo de boucles dirigées (Directed loop algorithm) tandis que le QDM nécessite un nouveau schéma de Monte Carlo quantique basé sur le cadre des expansions en séries stochastiques appelé Sweeping Cluster Algorithm (SCA). Nous présentons une modification du SCA pour rendre les simulations ergodiques à température finie. Notre étude numérique montre l'existence d'une phase critique dans les deux modèles séparée par une transition de phase à température finie vers une phase ordonnée de dimères ou de boucles qui brise spontanément certaines symétries du réseau. Nous montrons que dans le cas où les interactions sont attractives, cette transition de phase est de type Kosterlitz-Thouless et peut être comprise en construisant une théorie des champs de hauteurs. Le diagramme de phase à température finie du QDM présente quant à lui une ré-entrance inhabituelle dans la phase critique. La partie finale de cette thèse traite du rôle des symétries non-abéliennes dans les processus de thermalisation des systèmes quantiques. Nous étudions les états propres de haute énergies d'une chaîne de spins désordonnée avec une symétrie SU(3). Le modèle n'a pas directement de contraintes mais nous effectuons une diagonalisation exacte dans une base contrainte de tableaux de Young en utilisant la symétrie SU(3) complète de ce modèle. En examinant les sondes couramment utilisées pour la thermalisation (statistiques spectrales, distribution des observables locales et entropie d'intrication), nous montrons que le modèle présente un régime non ergodique sur une grande plage de tailles de système pour un désordre fort, contrastant avec la thermalisation rapide observée à faible désordre.