Dissertations / Theses on the topic 'Hartree-Fock-Bogoliubov'

Cette thèse est consacrée à l'étude mathématique et numérique du modèle de Hartree-Fock-Bogoliubov (HFB) pour les systèmes quantiques attractifs, qui est abondamment utilisé en physique nucléaire. Après avoir présenté le modèle et ses principales caractéristiques, nous expliquons comment le discrétiser et nous montrons des résultats de convergence. Nous examinons tout particulièrement l'algorithme de point fixe (parfois appelé Roothaan) et montrons qu'il converge ou alors oscille entre deux états dont aucun n'est solution du problème. Ceci généralise au cadre HFB des résultats de Cancès et Le Bris pour le modèle plus simple de Hartree-Fock dans le cas répulsif. Suivant ces mêmes auteurs, nous proposons un algorithme basé sur la contrainte relâchée et pour lequel la convergence est garantie. Dans la dernière partie de la thèse, nous illustrons le comportement de ces algorithmes par des simulations numériques pour plusieurs modèles. Dans un premier temps nous considérons un système purement gravitationnel où les particules interagissent avec le potentiel de Newton. Nos simulations montrent que la matrice d'appariement est toujours non nulle, un fait qui n'a pas encore pu être démontré rigoureusement. Nous étudions ensuite un modèle très simplifié pour la description de protons et neutrons dans le noyau atomique.

La fission des noyaux est un phénomène extrêmement complexe et un modèle réaliste décrivant l'ensemble du processus n'existe pas encore à l'heure actuelle. Ce travail porte sur la description théorique des distributions de fragments de fission, en masse, charge, énergie et déformation. Il s'agit de réactualiser et d'enrichir le modèle de point de scission de B. D. Wilkins. Nous voulons tester si un modèle statistique, appliqué au point de scission et utilisant les résultats des calculs microscopiques modernes, permet une description quantitative des distributions des fragments de fission. Nous calculons les surfaces d'énergie disponible au point de scission, en fonction de la déformation des fragments. Ces surfaces sont obtenues à partir de calculs microscopiques Hartree-Fock-Bogoliubov contraints, qui assurent un traitement réaliste de la variation du potentiel de chaque fragment en fonction de sa déformation. La pondération statistique est décrite par les densités de niveaux des fragments. Nous nous attachons à éviter autant que possible d'inclure des paramètres empiriques dans notre modèle. Le choix de la distance entre les fragments au point de scission est discuté sur la base des configurations de scission obtenues dans les calculs dynamiques totalement microscopiques. La confrontation de nos résultats avec les données expérimentales, sur la base d'une large systématique, est très satisfaisante et nous amène à analyser les capacités et les limites de notre approche. Nous proposons diverses voies d'amélioration, notamment par l'ajout d'éléments de dynamique
Nucleus fission process is a very complex phenomenon and, even nowadays, no realistic models describing the overall process are available. The work presented here deals with a theoretical description of fission fragments distributions in mass, charge, energy and deformation. We have reconsidered and updated the B. D. Wilking Scission Point model. Our purpose was to test if this statistic model applied at the scission point and by introducing new results of modern microscopic calculations allows to describe quantitatively the fission fragments distributions. We calculate the surface energy available at the scission point as a function of the fragments deformations. This surface is obtained from a Hartree Fock Bogoliubov microscopic calculation which guarantee a realistic description of the potential dependence on the deformation for each fragment. The statistic balance is described by the level densities of the fragment. We have tried to avoid as much as possible the input of empirical parameters in the model. Our only parameter, the distance between each fragment at the scission point, is discussed by comparison with scission configuration obtained from full dynamical microscopic calculations. Also, the comparison between our results and experimental data is very satisfying and allow us to discuss the success and limitations of our approach. We finally proposed ideas to improve the model, in particular by applying dynamical corrections

La théorie covariance de la fonctionnelle de la densité (CDF), basée sur un petit nombre de paramètres ajustables, a été utilisée avec succès pour décrire l’état fondamental et les états excités des noyaux de la carte nucléaire, pour A>12. Cette approche permet de décrire les systèmes nucléaires finis avec un Lagrangien hadronique universel résolu dans le cadre de l'approche Relativiste-Hartree-Fock-Bologuibov (RHFB). Ce modèle est également utilisé pour l'étude des étoiles compactes, car il peut être étendu à des densités élevées où la relativité restreinte ne peut pas être ignoré. Ce modèle peut également être étendu pour inclure la contribution des hypérons et ainsi que d'autres particules exotiques. Dans ce travail, la description et des prédictions basées sur l'approche RHFB pour les noyaux dans des conditions extrêmes de la masse, d'isospin et de température sont présentés.Dans la première partie de cette thèse, nous explorons l'apparition de nouvelles fermetures de couches sphériques pour des noyaux super-lourds, où les fermetures de couches sont caractérisées en termes de gap à deux nucléons. Les résultats dépendent légèrement des Lagrangians effectifs utilisés, mais les nombres magiques au-delà de ^{208} Pb sont prédit pour un nombre de protons Z=120 et 138, et pour un nombre de neutrons N=172, 184, 228, et 258. Les effets de couche sont sensibles à différents termes de champ de moyen, tels que le couplage spin-orbite, la masse scalaire et la masses effective, ainsi que l'interaction de tensorielle de Lorentz. Ces termes ont des poids différents dans les Lagrangians effectifs employées, expliquant les variations, somme toute petites, dans leurs prédictions. Employant le modèle RHFB le plus avancé, nous avons trouvé que le nucléide ^{304} 120 est favorisée comme étant le prochain noyau sphérique doublement magique au-delà de ^{208} Pb.Dans la deuxième partie de cette thèse, nous étudions l'apparition de nouveaux nombres magiques pour les noyaux de masse intermédiaire riches en neutrons, et nous analysons le rôle des interactions pseudo-vecteur et de tensorielle de Lorentz. Basé sur la transformation de Foldy-Wouthuysen, nous discutons en détail le rôle joué par les différents termes des interactions pseudo-vecteur et de tensorielle de Lorentz. Dans l'apparition des nouveaux nombres magiques N=16 , 32 et 34. Les noyaux ^{24} O ^{48} Si et ^{52,54} Ca sont prédits avec un grand gap au niveau de Fermi et un gap d'appariement zéro (^{24} O,^{54}Ca ) ou quasi-nul (^{48} Si,^{54} Ca), les rendant candidats pour de nouveaux nombres magiques des noyaux riches en neutrons. Nous constatons que les interactions de Lorentz pseudo-vecteur et tensorielle induisent des évolutions très spécifiques des énergies à une particule, ce qui pourrait signer la présence et la nécessité d'approches relativistes avec des interactions d'échanges de mésons.Dans la dernière partie de cette thèse, nous étudions les transitions de phase et excitations thermiques des deux noyaux stables et faiblement liés. Les prédictions de divers Lagrangiens relativistes et des différentes interactions d'appariement sont discutées. La température critique de la transition d'appariement dépend linéairement du gap d'appariement à température nulle, et cette dépendance est similaire pour une interaction de portée nulle ou bien finie. Les calculs présentés montrent des caractéristiques intéressantes des corrélations d'appariement à température finie, comme la persistance d'appariement et les phénomènes de re-entrance superfluide. En outre, nous analysons la réponse thermique de certains noyaux.En conclusion, le travail présenté dans cette thèse montre des résultats très intéressants et nouveaux pour trois des questions les plus importantes en physique nucléaire: la quête d'un nouvel îlot de stabilité dans la région des super-lourds, l'apparition de nouveaux nombres magiques dans les noyaux exotiques, et la réponse d'un système finis aux excitations thermiques
The covariant density functional (CDF) theory with a few number of parameters has been successfully employed to describe ground-state and excited-states of nuclei over the entire nuclear landscape for A > 12. It describes finite nuclear systems with a universal hadronic Lagrangian, which is solved considering the relativistic-Hartree-Fock-Bologuibov approach (RHFB). This model is also employed for the study of compact stars, since it can be extended to high densities where special relativity cannot be ignore. This model can also be extended to include the contribution of hyperons and as well as other exotic particles. In this work, the description and some predictions based on RHFB approach for nuclei under extreme conditions of mass, isospin and temperature are presented.In the first part, we explore the occurrence of spherical shell closures for superheavy nuclei, where shell closures are characterized in terms of two-nucleon gaps. The results depend slightly on the effective Lagrangians used, but the magic numbers beyond ^{208}Pb are generally predicted to be Z = 120 and 138 for protons, and N = 172, 184, 228, and 258 for neutrons. Shell effects are sensitive to various terms of the mean-field, such as the spin-orbit coupling, the scalar and the effective masses, as well as the Lorentz-tensor interaction. These terms have different weights in the effective Lagrangians employed, explaining the (relatively small) variations in the predictions. Employing the most advanced RHFB model, we founded that the nuclide ^{304}120 is favored as being the next spherical doubly-magic nucleus beyond ^{208}Pb.In the second part, we investigate the formation of new shell gaps in intermediate mass neutron-rich nuclei, and analyze the role of the Lorentz pseudo-vector and tensor interactions. Based on the Foldy-Wouthuysen transformation, we discuss in detail the role played by the different terms of the Lorentz pseudo-vector and tensor interactions in the appearance of the N=16, 32 and 34 shell gaps. The nuclei ^{24}O, ^{48}Si and ^{52,54}Ca are predicted with a large shell gap and zero (^{24}O, ^{52}Ca) or almost zero (^{48}Si, ^{54}Ca) pairing gap, making them candidates for new magic numbers in neutron rich nuclei. We find that the Lorentz pseudo-vector and tensor interactions induce very specific evolutions of single-particle energies, which could clearly sign their presence and reveal the need for relativistic approaches with exchange interactions.In the last part, we study the phase transitions and thermal excitations of both stable and weakly-bound nuclei. The predictions of various relativistic Lagrangians and different pairing interactions are discussed. The critical temperature of the pairing transition is found to depend linearly on the zero-temperature pairing gap, and this dependence is similar for a zero-range or a finite-range pairing interaction. The present calculations show interesting features of the pairing correlations at finite temperature, such as the pairing persistence and pairing re-entrance phenomena. Also, we analyze the thermal response of some nuclei.In conclusion, the work presented in this thesis shown interesting and new results for three of the most important questions in nuclear physics: the quest for a new island of stability in the superheavy region, the appearance of new magic numbers in exotic nuclei, and the response of finite-systems to thermal excitations

Le travail présenté dans ce mémoire de thèse consiste au développement d'un modèle Hartree-Fock-Bogoliubov relativiste en symétrie axiale. Pour la première fois, un modèle de type champ moyen autocohérent relativiste incorporant explicitement les contributions d'échange est à même de décrire les noyaux atomiques déformés et superfluides. Le choix d'une formulation covariante n'est en aucun cas lié au besoin d'une cinématique nucléaire relativiste mais à l'exploitation des symétries de Lorentz. En particulier, le maintien d'une distinction entre scalaire de Lorentz et composante temporelle des 4-vecteurs est source d'efficacité et de simplification en ce qui concerne la description des systèmes nucléaires. Dans ce cadre, les nucléons sont considérés comme des particules ponctuelles représentées par un spineur de Dirac. Leur interaction dans le milieu nucléaire est decrit en terme d'échange de mésons effectifs. La dynamique de ces degrés de liberté est spécifiée par une densité lagrangienne phénoménologique, caractérisée par huit paramètres libres contraints par l'ajustement sur la masse de noyaux sphériques ainsi que par la reproduction de propriétés de la matière nucléaire. Une fonctionnelle de la densité pour l'énergie s'obtient à partir de la densité lagrangienne traitée dans l'approximation Hartree-Fock-Bogoliubov. En particulier, la prise en compte explicite des termes d'échange permet l'introduction du pion, qui ne contribue que par son terme de Fock à l'échelle du champ moyen. La minimisation de la fonctionnelle de la densité pour l'énergie conduit aux équations Hartree-Fock-Bogoliubov relativistes considérées en symétrie axiale. Ces dernières sont résolues dans une base d'oscillateur harmonique déformé. Le modèle ainsi développé est utilisé pour décrire l'état fondamental des isotopes de carbone, néon et magnésium.

L'interaction effective entre les nucléons est un ingrédient essentiel des
calculs microscopiques de structure nucléaire. L'une des formes utilisée
depuis les années 1970 est la force effective phénoménologique proposée par
D. Gogny. Cette force donne d'excellents résultats dans les noyaux à
l'approximation du champ moyen. La présence de composantes de contact ne
permet pas cependant de l'employer en toute généralité pour décrire les
corrélations au-delà du champ moyen qui se manifestent fréquemment dans les
noyaux.

Dans ce travail, nous étudions des extensions de la force de Gogny, notamment
une généralisation où le terme dépendant de la densité de portée nulle est
remplacé par un terme de portée finie. Les paramètres intervenant dans la forme
analytique de la force sont ajustés sur les propriétés de la matière nucléaire
infinie symétrique et de la matière neutronique, et sur les observables globales
de quelques noyaux stables, en particulier celles associées aux corrélations
d'appariement. Nous présentons la méthode permettant d'inclure ce type de
force dans les codes de calcul Hartree-Fock-Bogoliubov et nous analysons les
résultats obtenus dans de nombreux noyaux. Les nouvelles versions de la force
de Gogny apparaissent reproduire la structure nucléaire avec une qualité égale
ou supérieure à la version traditionnelle.

We study the application of Bose-Einstein condensates (BECs) to simulations of phenomena across a number of disciplines in physics, using theoretical and computational methods. ¶ Collapsing condensates as created by E. Donley et al. [Nature 415, 39 (2002)] exhibit potentially useful parallels to an inflationary universe. To enable the exploitation of this analogy, we check if current quantum field theories describe collapsing condensates quantitatively, by targeting the discrepancy between experimental and theoretical values for the time to collapse. To this end, we couple the lowest order quantum field correlation functions to the condensate wavefunction, and solve the resulting Hartree-Fock-Bogoliubov equations numerically. Complementarily, we perform stochastic truncated Wigner simulations of the collapse. Both methods also allow us to study finite temperature effects. ¶ We find with neither method that quantum corrections lead to a faster collapse than is predicted by Gross-Pitaevskii theory. We conclude that the discrepancy between the experimental and theoretical values of the collapse time cannot be explained by Gaussian quantum fluctuations or finite temperature effects. Further studies are thus required before the full analogue cosmology potential of collapsing condensates can be utilised. ¶ As the next project, we find experimental parameter regimes in which stable three-dimensional Skyrmions can exist in a condensate. We show that their stability in a harmonic trap depends critically on scattering lengths, atom numbers, trap rotation and trap anisotropy. In particular, for the Rb87 |F=1,m_f=-1>, |F=2,m_f=1> hyperfine states, stability is sensitive to the scattering lengths at the 2% level. We find stable Skyrmions with slightly more than 2*10^6 atoms, which can be stabilised against drifting out of the trap by laser pinning. ¶ As a stepping stone towards Skyrmions, we propose a method for the stabilisation of a stack of parallel vortex rings in a Bose-Einstein condensate. The method makes use of a ``hollow'' laser beam containing an optical vortex, which realises an optical tunnel for the condensate. Using realistic experimental parameters, we demonstrate numerically that our method can stabilise up to 9 vortex rings. ¶ Finally, we focus on analogue gravity, further exploiting the analogy between flowing condensates and general relativistic curved space time. We compare several realistic setups, investigating their suitability for the observation of analogue Hawking radiation. We link our proposal of stable ring flows to analogue gravity, by studying supersonic flows in the optical tunnel. We show that long-living immobile condensate solitons generated in the tunnel exhibit sonic horizons, and discuss whether these could be employed to study extreme cases in analogue gravity. ¶ Beyond these, our survey indicates that for conventional analogue Hawking radiation, simple outflow from a condensate reservoir, in effectively one dimension, has the best properties. We show with three dimensional simulations that stable sonic horizons exist under realistic conditions. However, we highlight that three-body losses impose limitations on the achievable analogue Hawking temperatures. These limitations vary between the atomic species and favour light atoms. ¶ Our results indicate that Bose-Einstein condensates will soon be useful for interdisciplinary studies by analogy, but also show that the experiments will be challenging.

La fission induite par neutron, découverte il y a plus de 70 ans, a de nombreuses applications, par exemple industrielles pour la production d'énergie, et intervient dans la nucléosynthèse. Cependant, sa description microscopique reste un problème ouvert. En effet, les degrés de liberté qui interviennent dans ce processus dynamique sont complexes. De plus, les noyaux fissiles ont un nombre élevé de nucléons en interaction (>200). Il s'agit donc d'un problème à N-corps quantique. Or, une résolution directe de ce dernier n'est pas possible à l'heure actuelle. Dans ce contexte, la description microscopique de la fission considérée ici est la suivante : la première étape consiste à déterminer un ensemble de configurations de champ moyen qui représentent différentes déformations du noyau, incluant ainsi explicitement les degrés de liberté collectifs qui leur sont associés. Dans la seconde étape, la dynamique est décrite dans cet espace de configurations en utilisant la méthode de la coordonnée génératrice dépendante du temps (TDGCM). L'approximation des recouvrements gaussiens (GOA) est alors utilisée. Cependant, elle introduit une erreur de modèle et limite les extensions comme par exemple la prise en compte explicite de degrés de liberté intrinsèques. Ce travail de thèse a pour objectif de décrire le processus de fission avec la TDGCM sans recourir à la GOA. Cela implique de résoudre l'équation de la dynamique en TDGCM appelée équation de Hill-Wheeler dépendante du temps (TD-HW). Les méthodes d'évaluations des matrices des recouvrements et du hamiltonien collectif sont présentées dans le cas d'une interaction de Gogny. La matrice des recouvrements représente la métrique de l'espace des configurations, et la matrice du hamiltonien collectif contient les couplages énergétiques entre les configurations. Les configurations sont exprimées dans des bases de particules deux à deux distinctes, introduisant des instabilités numériques dans les méthodes d'évaluation standard. Un formalisme adapté à ces bases est proposé permettant d'éliminer ces instabilités. Deux méthodes de résolution de TD-HW sont présentées. La première consiste à calculer l'opérateur d'évolution associé à l'équation de Hill-Wheeler dépendante du temps. Elle est adaptée à un faible nombre de configurations. La seconde utilise un schéma de discrétisation en temps permettant l'inclusion d'un plus grand nombre de configurations dans le modèle. Ce formalisme est ensuite appliqué à la description de la réaction de fission induite par neutron sur le plutonium 239, et une comparaison avec la TDGCM+GOA est effectuée
Nuclear fission, where an atomic nucleus separates into two fragments while emitting a large amount of energy, is at the core of many applications in society (energy production) and national security (deterrence, non-proliferation). It is also a key ingredient of the mechanisms of formation of elements in the universe. Yet, nearly 80 years after its experimental discovery its theoretical description in terms of the basic constituents of the nucleus (protons and neutrons) and their interaction remains a challenge. In this thesis, we describe the fission process as follows. In a first step, we use large supercomputers to compute the deformation properties of the nucleus based on our knowledge of nuclear forces. In a second step, we simulate the time evolution of the system from its ground state up to the fragments separation with a fully quantum-mechanical approach called the time-dependent generator coordinate method (TDGCM). While results are in good qualitative agreement with experimental data, the implementation of the TDGCM so far had been greatly simplified using what is known as the Gaussian overlap approximation (GOA). We also developed the formalism and a numerical implementation of the exact TDGCM - without the GOA. This will allow the first systematic validation of that approximation and an assessment of the resulting theoretical uncertainties. The second chapter presents the description of the neutron induced fission process using the TDGCM+GOA. The third one introduces the developments carried out in this thesis allowing the description of the fission process with the TDGCM without the GOA. The last chapter shows the first results obtained with this approach

Les noyaux riches en neutrons de masse A=100-110 constituent une région de grand intérêt pour l'étude de la structure nucléaire loin de la stabilité. De précédentes études de cette région de masse ont déjà révélé la complexité de l'évolution de la collectivité et de la déformation dans les chaînes isotopiques de Zr, Mo, Ru et Pd. Afin d'étendre les données expérimentales sur la collectivité à des états de plus haut spin et à des noyaux plus riches en neutrons, des temps de vie d'états excités ont été mesurés dans des noyaux produits par une réaction de fusion-fission en cinématique inverse au GANIL. Les fragments de fission étaient séparés et identifiés en A et Z grâce au spectromètre magnétique de grande acceptance VAMOS tandis que le rayonnement gamma était détecté dans l'ensemble de détecteurs germanium EXOGAM. Environ vingt temps de vie d'états 2+, 4+ et 6+ ont été extraits à l'aide du plunger de Cologne. Cette expérience représente la première mesure RDDS dans des fragments de fission identifiés évènement par évènement à la fois en A et Z.Cette étude des noyaux de masse A=100-110 est complétée par des calculs auto-cohérents de champ moyen et au-delà avec la force de Gogny (D1S). La structure des états fondamentaux et excités est décrite dans le cadre du modèle de Hartree-Fock-Bogoliubov avec des contraintes sur les déformations axiale et triaxiale. Les excitations individuelles sont étudiées par des calculs bloqués et les états de haut spin sont décrits dans l'approximation du champ tournant. Enfin, la méthode de la coordonnée génératrice approchée par un hamiltonien collectif en cinq dimensions (5DCH) est appliquée aux états collectifs de basse énergie. Les résultats sont comparés aux mesures de la collectivité
Neutron-rich nuclei of mass A=100-110 are of great interest for the study of nuclear structure far from stability. Previous experimental and theoretical studies suggest a complex evolution of deformation and collectivity in the isotopic chains of Zr, Mo, Ru and Pd.In order to extend information on the evolution of the collectivity towards higher spin states and more neutron-rich nuclei, lifetimes of excited states were measured in nuclei produced through a fusion-fission reaction in inverse kinematic at GANIL. Fission fragments were separated and identified in both A and Z with the high acceptance magnetic spectrometer VAMOS while the EXOGAM germanium detectors array was used for the coincident gamma-ray detection. Lifetimes of about twenty excited states were extracted using the plunger device of Cologne. This is the first RDDS measurement on fission fragments which are identified in A and Z on an event-by-event basis. The study of this mass region is completed by theoretical calculations using self consistent mean field and beyond mean field methods implemented with the Gogny force (D1S). The structure of the ground states and the excited states is described with Hartree-Fock-Bogoliubov calculations with constraints placed on the axial and triaxial deformations. Individual excitations are investigated through blocking calculations and the high spin states are studied through cranking calculations. Finally, an approximated generator coordinate method (GCM+GOA) using the 5DCH hamiltonian is used to describe the low energy collective states and to interpret the experimental evolution of the collectivity

Ce mémoire traite de différents aspects des phénomènes collectifs rencontrés en dynamique nucléaire. Le formalisme de la densité fonctionnelle d'énergie dépendante du temps (TDEDF) constitue un des modèles les plus aboutis pour réaliser de telles études. Une des premières implémentations du formalisme TDEDF permettant de traiter l'appariement nucléaire, basée sur la résolution des équations Hartree-Fock-Bogoliubov dépendantes du temps, est présentée. Dans une seconde partie, le formalisme TDEDF est appliqué à l'étude de certaines excitations collectives du noyau. Tout d'abord, l'accent est mis sur les propriétés de structure et de décroissance de la résonance géante monopolaire pour des noyaux à (sous-)couche fermée. Un autre mode d'excitation du noyau atomique est ensuite étudié : les vibrations d'appariement. L'impact des fonctionnelles utilisées a été discuté. La troisième partie est dédiée à l'étude de la synthèse des éléments les plus lourds par fusion d'ions lourds presque symétriques. Ces réactions sont sujettes au phénomène de suppression de la fusion. Ce dernier est illustré au moyen du formalisme TDEDF. Une expérience test sur la fusion par voie symétrique, menée au GANIL sur le système Xe{136}+Sn{124}->Rf{260}*, est ensuite analysée. Aucun résidu d'évaporation n'a été observé. Une section efficace maximale de fusion-évaporation de 81 pb a été obtenue pour la voie 2n. Les difficultés expérimentales rencontrées sont discutées et mises en perspectives avec l'avènement d'une nouvelle génération de faisceaux d'ions radioactifs tels que ceux qui seront disponibles avec Spiral2.

La détermination de la masse nucléaire permet d'accéder à l'énergie de liaison donc au bilan de toutes les forces agissant dans le noyau. L'exploration de la structure nucléaire requiert la mesure des masses loin de la stabilité. Néanmoins la production des ions radioactifs est peu sélective et l'ion d'intérêt est souvent noyé dans un fond de contaminants isobariques. Ces dernières années, l'utilisation des pièges à ions linéaires et du refroidissement par gaz tampon a été fortement développée permettant la mesure de masse des noyaux exotiques avec une grande précision. Cette thèse porte sur le développement innovant d'un piège de Paul circulaire, "le cirque d'ions". Sa géométrie lui permet de confiner les ions sur un grand nombre de tours. Ainsi, ils peuvent être refroidis avec une pression de gaz tampon de l'ordre de 10-4 mbar. Par ailleurs, il est capable de trier en masse les ions piégés et refroidis avec un pouvoir de résolution suffisant pour séparer des isobares. Nous avons réalisé et testé le premier prototype de ce spectromètre à Orsay. Cette thèse étudie également les prédictions des masses par les différents modèles théoriques, en les comparant avec l'approche microscopique HFB-17. Actuellement, la masse des noyaux très exotiques peut seulement être prédite. L'étude de l'énergie de séparation de deux neutrons avec HFB-17 met en évidence des zones d'instabilité numérique. Afin de l'améliorer, nous avons apporté une correction aux énergies de déformation calculées par ce modèle, en utilisant comme référence les masses des chaînes isotopiques 143-146Xe et 223-229Rn que nous avons mesurées avec le spectromètre ISOLTRAP au CERN.

Ce mémoire est consacré à l'étude mathématique de divers modèles variationnels permettant la description de systèmes quantiques, en particulier infinis. Les outils mathématiques utilisés sont ceux de l'analyse non linéaire, du calcul des variations, des équations aux dérivées partielles, de la théorie spectrale et du calcul scientifique.

Une première partie contient quelques résultats pour des systèmes finis. Nous étudions des approximations de l'équation de Schrödinger pour N électrons dans une molécule ou un atome, puis le modèle de Hartree-Fock-Bogoliubov pour un système de fermions interagissant avec une force de type gravitationnelle.

Dans une seconde partie nous proposons une nouvelle méthode pour démontrer l'existence de la limite thermodynamique pour des systèmes quantiques interagissant avec la force de Coulomb.

Ensuite, nous construisons deux modèles de type Hartree-Fock pour des systèmes infinis. Le premier est un modèle relativiste, déduit de l'électrodynamique quantique, et qui permet de décrire le comportement d'électrons, couplés avec celui du vide de Dirac qui peut se polariser. Le second modèle décrit l'état d'un cristal non relativiste en présence d'un défaut chargé ; il est complété par une nouvelle approche numérique.

La dernière partie du mémoire est consacrée au problème de pollution spectrale, un phénomène observé lorsque l'on cherche à calculer des valeurs propres au milieu du spectre essentiel, par exemple pour des opérateurs de Dirac ou de Schrödinger périodique.

Strong electron correlations in electronic structure theory are purely quantum effects arising as a result of degeneracies in molecules and materials, and exhibit significantly different yet interesting characters than do weak correlations. Although weak correlations have recently been able to be described very efficiently and accurately within single particle pictures, less known are good prescriptions for treating strong correlations efficiently. Brute-force calculations of strong correlations in wave function theories tend to be very computationally-intensive, and are usually limited to small molecules for applications. Breaking symmetry in a mean-field approximation is an efficient alternative to acquire strong correlations with, in many cases, qualitatively accurate results. The symmetry broken in quantum chemistry has been traditionally of spin, in so-called unrestricted methods, which typically break spatial symmetry as a consequence, and vice versa, in most situations. In this work, we present a novel approach to accurately describing strong correlations with a mean-field cost by means of Hartree- Fock-Bogoliubov (HFB) theory. We are inspired by the number-symmetry-breaking in HFB, which, with an attractive particle interaction, accounts for strong correlations, while maintaining spin and spatial symmetry. We show that this attractive interaction must be restricted to the chemically-relevant orbitals in an active space to obtain physically meaningful results. With such constraints, our constrained pairing mean-field theory (CPMFT) can accurately describe potential energy curves of various strongly-correlated molecular systems, by cleanly separating strong and weak correlations. To achieve the correct dissociation limits in hetero-atomic molecules, we have modified our CPMFT functional by adding asymptotic constraints. We also include weak correlations by combining CPMFT with density functional theory for chemically accurate results, and reveal the connection between CPMFT and traditional unrestricted methods. The similarity between CPMFT and unrestricted methods leads us to the idea of constrained active space unrestricted mean-field approaches. Motivated by CPMFT, we partially retrieve spin-symmetry that has been fully broken in unrestricted methods. We allow symmetry breaking only in an active space. This constrained unrestricted Hartree-Fock (CUHF) is an interpolation between two extrema: the fully broken-symmetry solution and the symmetry preserved solution. This thesis defines the theory behind and reports the results of CUHF. We first show that, if an active space is chosen to include only open-shell electrons, CUHF reduces to restricted open-shell Hartree-Fock (ROHF), and such CUHF proves in many ways significantly