Bogoliubov Many-Body Perturbation Theory for Nuclei : Systematic Generation and Evaluation of Diagrams and First ab initio Calculations Les dernières décennies ont donné lieu à un développement rapide des théories ab initio visant à décrire les propriétés des noyaux à partir de l'interaction nucléonique. Un tel développement a été rendu possible à la fois par la très importante croissance de la puissance de calcul et de nouveaux développements formels. Le présent travail se consacre au développement de la théorie de perturbation à N corps de Bogolioubov récemment proposée, qui repose sur l'usage d'un état de référence brisant la symétrie associée au nombre de particules pour permettre une description des noyaux à simple couche ouverte. Le formalisme est tout d'abord décrit en détails, son lien avec la théorie de perturbation à N corps standard est établi, tout comme sa connexion avec la théorie de cluster couplés de Bogolioubov. L'extension du formalisme à des ordres plus élevés à partir de méthodes de théorie des graphes est ensuite présentée ainsi que le programme ADG qui génère et évalue les diagrammes BMBPT à un ordre quelconque. Les implications de ce développement formel dépassent le cadre du présent travail, les méthodes développées pouvant être appliqués à d’autres méthodes à N corps. Pour terminer, de premiers résultats numériques pour les isotopes de l'oxygène, du calcium et du nickel sont présentés. Ces résultats établissent la théorie de perturbation à N corps de Bogolioubov comme une méthode de premier intérêt pour des calculs à grande échelle sur les chaînes isotopiques et isotoniques de masse moyenne.
Bogoliubov Many-Body Perturbation Theory for Nuclei : Systematic Generation and Evaluation of Diagrams and First ab initio Calculations Les dernières décennies ont donné lieu à un développement rapide des théories ab initio visant à décrire les propriétés des noyaux à partir de l'interaction nucléonique. Un tel développement a été rendu possible à la fois par la très importante croissance de la puissance de calcul et de nouveaux développements formels. Le présent travail se consacre au développement de la théorie de perturbation à N corps de Bogolioubov récemment proposée, qui repose sur l'usage d'un état de référence brisant la symétrie associée au nombre de particules pour permettre une description des noyaux à simple couche ouverte. Le formalisme est tout d'abord décrit en détails, son lien avec la théorie de perturbation à N corps standard est établi, tout comme sa connexion avec la théorie de cluster couplés de Bogolioubov. L'extension du formalisme à des ordres plus élevés à partir de méthodes de théorie des graphes est ensuite présentée ainsi que le programme ADG qui génère et évalue les diagrammes BMBPT à un ordre quelconque. Les implications de ce développement formel dépassent le cadre du présent travail, les méthodes développées pouvant être appliqués à d’autres méthodes à N corps. Pour terminer, de premiers résultats numériques pour les isotopes de l'oxygène, du calcium et du nickel sont présentés. Ces résultats établissent la théorie de perturbation à N corps de Bogolioubov comme une méthode de premier intérêt pour des calculs à grande échelle sur les chaînes isotopiques et isotoniques de masse moyenne.
Théorie de perturbation à N corps de Bogoliubov projetée : Solutions aux difficultés formelles et techniques Étant donnée l’interaction entre les nucléons,la résolution de l’équation de Schrödinger à A corps permet d’accéder aux propriétés des états quantiques des noyaux. La théorie des perturbations à N corps (MBPT) basée sur un état de champ moyen d’Hartree-Fock permet de traiter les corrélations dynamiques mises en jeu dans les noyaux à doubles couches fermées. Au-delà des fermetures de couches, la brisure spontanée de la symétrie U(1)associée à la conservation du nombre de protons et de neutrons permet en sus l’inclusion des corrélations statiques au niveau du champ moyen et la formulation de la théorie des perturbations à N corps de Bogoliubov (BMBPT). Néanmoins, les résultats des calculs BMBPT présentent des contaminations du fait que la brisure de symétrie n’est qu’émergente dans les systèmes quantiques finis tels que le noyau atomique. Ainsi, la restauration de la symétrie U(1) au-delà du champ moyen est nécessaire pour une description correcte et donne lieuà la formulation de la méthode BMBPT projetée(PBMBPT). Le but est d’implémenter PBMBPTafin d’effectuer des calculs ab initio pour les noyaux à couche ouverte et de masse intermédiaire. Le présent travail apporte des solutions systématiques aux problèmes formels et techniques qui se posent lorsde l’implémentation numérique de PBMBPT : une méthode automatique et sûre permettant de générer les diagrammes PBMBPT ainsi que les expressions correspondantes, un formalisme permettant de calculer d’autres observables que l’énergie, une extension de l’approximation en produit normal aux méthodes avec brisures de symétries ainsi qu’un outil automatique et sûr de réduction sphérique de réseaux de tenseurs.
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.
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.
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.
Ce mémoire est une synthèse, sous forme de sélection d'articles, des travaux développés depuis ma soutenance de thèse en 2001. De nombreuses lignes de recherche se sont dégagées dans les années qui ont suivi ma thèse et mes activités scientifiques se sont développées en suivant différentes directions. Le noyau atomique est un système quantique composé par N fermions qui nécessite un traitement théorique basé sur des modèles microscopiques. Les travaux développés et l'expérience acquise dans ce domaine se prêtent donc de manière très naturelle à des extensions interdisciplinaires pour l'étude d'autres systèmes quantiques formés par N fermions comme, par exemple, les agrégats métalliques ou les gaz d'atomes fermioniques piégés. J'ai exploré quelques-uns de ces liens interdisciplinaires et l'ensemble de ces activités a contribué à former ma vision générale du problème à N corps. Cela a par ailleurs enrichi et élargi mes perspectives de recherche en physique nucléaire. La plus longue partie du mémoire est consacrée à la réflexion autour des noyaux exotiques effectuée avec des modèles basés sur le champ moyen pour étudier : l'évolution de la structure loin de la stabilité et certains phénomènes exotiques liés à cette évolution, les corrélations d'appariement, les halos, l'effet sur les propriétés de l'état fondamental du terme tenseur et d'autres termes non standards dans l'interaction de Skyrme. Des applications dans les domaines des réactions et de l'astrophysique nucléaire sont aussi discutées. Quelques extensions de la méthode RPA sont détaillées. Ces travaux visent à une formulation complètement auto-consistante de la méthode pour un meilleur traitement des corrélations dans l'état fondamental des systèmes à N corps. Dans les conclusions et perspectives du mémoire, les limites de validité de l'approximation de champ moyen sont mises en évidence et un certain nombre de possibilités d'amélioration du pouvoir prédictif et de raffinement des modèles sont présentées et discutées.
Les methodes variationnelles se sont revelees etre de puissants outils pour evaluer des quantites physiques dans des systemes de mecanique quantique ayant un nombre fini des degres de liberte, dans les cas ou la theorie des perturbations est inapplicable. Les applications de ces methodes, notamment en physique atomique, nucleaire et du solide, ont donne d'excellents resultats. Par consequent, il est tentant d'essayer de les adapter dans le cadre de la theorie quantique des champs. La premiere partie de ce travail consiste a introduire une nouvelle approche variationnelle-perturbative applique a un systeme de mecanique quantique bien connu: l'oscillateur anharmonique. De maniere tout a fait remarquable, cette methode a fourni une amelioration de sept ordres de grandeur a la precision caracteristique des resultats de hartree-fock. D'un autre cote, cette nouvelle methode semble etre bien adaptee pour des calculs de theorie quantique de champs. On espere qu'en combinant les estimations variationnelles avec le developpement perturbatif et la theorie de renormalisation, on pourra obtenir des reponses variationnelles finies. La deuxieme partie de ce travail est consacree a cet objectif. On a donc applique cette strategie a un modele non trivial et renormalisable a deux dimensions a savoir le modele de gross-neveu. On a trouve donc que les reponses variationnelles de diverses quantites physiques ne presentent pas des divergences. De plus, ces resultats sont assez raisonnables, si on les compare avec les resultats exacts de la litterature, et peuvent s'ameliorer de facon systematique