Les systèmes à électrons fortement corrélés sont d’intérêt particulier pour le calcul ab initio, cherchant à modéliser ces systèmes à partir des premiers principes. La théorie de la fonctionnelle de la densité associée à une prise en compte des corrélations locales en DFT+U ou en DFT+DMFT, permet de reproduire qualitativement la physique de ces systèmes. Cependant, ces méthodes font intervenir les paramètres d’interaction effective de Hubbard U et de Hund J. Ces derniers peuvent eux-mêmes être calculés de manière ab initio, notamment avec l’approximation de phase aléatoire contrainte (cRPA), ouvrant la voie au développement de schémas de calcul les plus prédictifs possible. Nous utilisons un schéma DFT+U/cRPA, dont le principe consiste à calculer les paramètres U et J en cRPA et la structure électronique en DFT+U de manière auto-cohérente. Nous appliquons ce schéma aux lanthanides allant du cérium au lutétium (en détaillant le cas du cérium dans ses phases gamma et alpha), et aux dioxydes des actinides allant de l'uranium au curium. Nous effectuons d'abord une étude de l'état fondamental en DFT+U en fonction de U, en détaillant l'influence des états métastables. Nous étudions le rôle de la localisation des orbitales corrélées sur l'interaction effective dans un cas particulier. Nous détaillons ensuite les valeurs de U obtenues en cRPA en fonction de celles utilisées pour le calcul DFT+U. Nous nous intéressons plus particulièrement aux effets des processus d’écrantage sur les valeurs de U obtenues. Nous montrons que les limitations du schéma observées (multiplicité des valeurs auto-cohérentes de U obtenues et/ou incompatibilité de ces valeurs avec les spectres de photoémission expérimentaux) sont causées par certains processus d’écrantage spécifiques. Cela suggère d'améliorer la description de ces processus d'écrantage en modifiant le modèle. / Strongly correlated electron systems are particularly interesting for ab initio calculus, which aims to model these systems from first principles. Density functional theory, improved by taking into account local correlations within DFT+U or DFT+DMFT, enables to reproduce qualitatively the physics of these systems. Nonetheless, these methods require the Hubbard and Hund effective interaction parameters U and J. These can be computed from first principles as well, within the constrained random phase approximation (cRPA), paving the way for numerical schemes as predictive as possible. Here, we use a DFT+U/cRPA scheme, which aims to compute self-consistently the U and J parameters with cRPA and the DFT+U electronic structure. We apply this scheme to lanthanides from cerium to lutetium (insisting on cerium and its gamma and alpha phases), and dioxides of actinides from uranium to curium. We study the DFT+U ground state in function of U, giving more details about the influence of metastable states. We study as well the influence of the localization of correlated orbitals on the effective interaction in a particular case. We then detail the values of U obtained with cRPA, in function of those used for the DFT+U calculation. We study more particularly the effects of screening processes on the obtained values of U. We show that the limitations of the scheme (multiplicity of self-consistent values of U and/or their incompatibility with experimental photoemission spectra) are caused by specific screening processes. This suggests to improve the description of these screening processes by modifying the model.

No description available.

Si le noyau atomique peut sembler aujourd'hui un système bien connu, sa description microscopique théorique comporte toujours des lacunes : décrire la structure nucléaire des noyaux tout le long de la charte des noyaux de manière satisfaisante par une seule interaction s'avère encore compliqué aujourd'hui.Les interactions les plus performantes et les plus utilisées pour atteindre ce but sont les interactions effectives. On peut généralement classer celles-ci en deux catégories: celles de portée finie et celles de portée nulle. L'interaction de Skyrme, dont il est principalement question dans ce manuscrit, appartient à la deuxième catégorie. Apparue dans les années 50 comme un développement en gradients à l'ordre 2, sa simplicité en a longtemps fait une interaction de choix. Une collaboration internationale, UNEDF, a toutefois montré que l'interaction en l'état était trop limitée pour reproduire correctement les données expérimentales actuelles et qu'il était nécessaire d'y trouver des extensions. Ce manuscrit propose une telle extension, décrivant le processus aboutissant à l'utilisation d'un développement en gradients jusqu'à l'ordre 6.Après avoir présenté les principales forces effectives utilisées dans le cadre de la structure nucléaire, nous montrerons le comportement de cette nouvelle interaction dans le modèle idéal de la matière infinie et son application dans la description des étoiles à neutrons. Nous verrons ainsi que les performances de notre interaction se sont avérées prometteuses. La comparaison des mêmes propriétés pour des potentiels de portée finie a validé ce constat. Notre extension à l'interaction de Skyrme offre à minima des performances comparables dans ce modèle.En conséquence, nous détaillerons également dans la suite du manuscrit la manière dont nous avons poursuivi le développement de notre nouvelle interaction en l'adaptant pour la description de noyaux sphériques. La dernière partie présente le processus d'ajustement de notre interaction, visant à une description des noyaux le long de la charte nucléaire. Une paramétrisation issue d'un tel ajustement sera enfin proposée / The nucleus may seem well known, but its microscopic description still isn't perfect. Describing accurately the nuclei along the nuclear chart remain complicated.Nowadays, the effective interactions are the best choice to achieve such a goal, and can be sorted in two categories: finite-range and zero-range. Skyrme interaction, the main subject of this report, belongs to the latter. Used since the 50s because of its simplicity, it is an expansion of gradients up to the 2th order. However, an international collaboration, UNEDF, showed recently that this interaction was too limited to describe all the nuclear constraints given by experimental datas. The conclusion was the necessity of Skyrme interaction extensions.This thesis is describing the process from the design of the interaction to its use in nuclei and in fitting procedures.First, we review the requirements and the construction of a general nuclear interactions, and the main interactions used in nuclear structure, both finite and zero range.We then show the behaviour of these interactions in a perfect model: the infinite nuclear matter. This is considered a preliminary test for every interaction aiming to describe nuclei, and the performances of our interactions are discussed, along the performances of the other models.The last part of the thesis is the adaptation of our extension for calculations of nuclei. Furthermore, the fitting procedure of our extended Skyrme interaction on both infinite nuclear matter and nuclei quantities is detailed. A resulting parametrisation is presented at the end of the thesis, aiming to describe more accurately the nuclear chart and the additional physic provided by the new terms

Cette thèse présente trois aspects centrés autour de la QRPA (Quasiparticle Random Phase Approximation).Le premier consiste en l'utilisation d'un code à symétrie axiale pour confronter des données calculéesà des résultats expérimentaux, et pour alimenter un code microscopique de réactions. Cette étapeest l'occasion d'analyser la spectroscopie du noyau à basse énergie (quelques dizaines de MeV), etplus spécifiquement (mais pas uniquement) la chaîne isotopique de l'étain (Z=50). Le second facetterepose sur l'amélioration d'un formalisme de calcul des opérateurs de transitions éléctromagnétiquesmultipolaires, et d'une méthode de généralisation du calcul de ces opérateurs permettant de faciliterla programmation en uniformisant le code pour les différentes multipolarités. Finalement, afin dedépasser la contrainte de la symétrie axiale, un nouveau code de calcul en symétrie triaxiale a étédéveloppé. Ses caractéristiques et son développement sont présentés, suivis des premiers résultatsdus à son exploitation. / This thesis presents three aspects centered around the QRPA (Quasiparticle Random Phase Approximation).The first consists in the use of an axial code to confront computed data with experimental results andto feed a microscopic reaction code. This step is a chance to analyse low-energy spectroscopy (fewtens of MeV) of some nuclei, and more precisely (but not exclusively) the tin isotopic chain (Z=50).The second one relies on the improvement of the formalism to calculate multipolar electromagnetictransition operators, and a method to consolidate the computation of these operators, allowing toease the programming by unifying the code for different multipolarities. Finally, in order to overcomethe axial symmetry constraint, a new triaxial code has been developed. Its assets and developmentare presented, followed by the first batch of results.

Materials intended for disposal of spent nuclear fuel require a particular combination of physical and chemical properties. The driving forces and mechanisms underlying the material’s behavior must be scientifically understood in order to enable modeling at the relevant time- and length-scales. The processes that determine the mechanical behavior of copper canisters and iron inserts, as well as the evolution of their mechanical properties, are strongly dependent on the properties of various defects in the bulk copper and iron alloys. The first part of the present thesis deals with precipitation in the cast iron insert. A nodular cast iron insert will be used as the inner container of the spent nuclear fuel. Precipitation is investigated by computing effective interaction energies for point defect pairs (solute–solute and vacancy–solute) in bcc iron using first-principles calculations. The main considered impurities in the iron matrix include 3sp (Si, P, S) and 3d (Cr, Mn, Ni, Cu) solute elements. By computing interaction energies possibility of formation of different second phase particles such as late blooming phases (LBPs) in the cast iron insert is evaluated. The second part is devoted to the fundamentals of dislocations and their role in plastic deformation of metals. Deformation of single-crystal copper under high strain rates is simulated by employing dislocation dynamics (DD) method to examine the effect of strain rate on mechanical properties as well as dislocation microstructure development. Creep deformation of copper canister at low temperatures is studied. The copper canister will be used in the long-term storage of spent nuclear fuel as the outer shell of the waste package to provide corrosion protection. A glide rate is derived based on the assumption that at low temperatures it is controlled by the climb rate of jogs on the dislocations. Using DD simulation creep deformation of copper at low temperatures is modeled by taking glide but not climb into account. Moreover, effective stresses acting on dislocations are computed using the data extracted from DD simulations. / <p>QC 20170428</p>

Alors qu’il existe de nombreux calculs microscopiques de barrières de fission pour des noyaux composés pair-pairs, il n’y a cependant que relativement peu de tels calculs pour des noyaux de masse impaire. Ceci est dû aux complications induites par la brisure de la symétric de reversement du sens du temps au niveau du champ moyen qui est engendrée par la présence d’un nucleon non apparié. Pour éviter cette difficulté, des calculs existants pour des noyaux de masse impaire ont tout simplement négligé ces effets de brisure de la symétrie de reversement du sens du temps.Dans ce travail, on se donne pour but d’améliorer la description des barrières de fission, aussi bien que des propriétés spectroscopiques du niveau fondamental et de l’état isomérique de fission,pour quelques isotopes de masse impaire dans la région des actinides en prenant en compte de tels effets. Ceci a été réalisé dans le cadre du formalisme de Skyrme–Hartree–Fock plus BCS avec blocking en adaptant ce formalisme à la brisure de la symétrie considérée. L’interaction résiduelle d’appariement a été approchée par une force de séniorité dont les paramètres ont été ajustés pour reproduire les différences de masse pair-impair de quelques noyaux de la région des actinides.Les énergies des têtes de bande rotationnelle de basse énergie ont été calculées dans le cadre du modèle unifié de Bohr-Mottelson pour quatre noyaux bien déformés (235U, 239Pu, 237Np, 241Am)produisant un bon accord qualitatif avec les données pour les noyaux impairs en neutrons. L’accord significativement moins bon obtenu pour les noyaux impairs en protons pourrait résulter de l’usage de l’approximation de Slater pour l’interaction d’échange de Coulomb. Les énergies de déformation de deux noyaux impairs en neutrons (235U, 239Pu) ont été calculées pour quelques configurations de particule individuelle, jusqu’après la barrières de fission externe. La symétrie axiale a été imposée tandis que la brisure de la symétrie droite-gauche (ou de parité intrinsèque) a été permise dans la région de la seconde barrière. Les hauteurs des barrières de fission pour ces noyaux impairs dépendent significativement des configurations de particule individuelle. Un accord qualitatif avec les données disponibles pour les hauteurs de barrières des noyaux impairs considérés et leurs voisins pair-pairs a été généralement obtenu. / While there have been numerous microscopic calculations on fission barriers of even-even compoundnuclei, there are however, relatively few such work dedicated to odd-mass nuclei. This is dueto the complications posed by the breaking of the time-reversal symmetry at the mean-field leveldue to the presence of an unpaired nucleon. In order to circumvent this difficulty, previous fission barriercalculations of odd-mass nuclei have been performed by neglecting the effect of time-reversalsymmetry breaking. This work aims to improve on the description of fission barriers as well asthe spectroscopic properties of ground and fission-isomeric state, of some odd-mass actinide nucleiby taking the effect of time-reversal symmetry breaking into account. This has been perfomedwithin a Skyrme-Hartree-Fock-plus-BCS framework with blocking, where the BCS formalism hasbeen adapted to accomodate this symmetry breaking. The Skyrme nucleon-nucleon effective forcehas been used with various sets of parameters (SIII, SkM*, SLy5*). The residual pairing interactionhas been approximated by seniority forces whose neutron and proton parameters have beenfitted to reproduce the odd-even mass differences of some actinide nuclei. The low-lying rotationalband-head energies evaluated within the Bohr-Mottelson unified model have been determined forfour well-deformed odd-nuclei (235U, 239Pu, 237Np, 241Am) yielding a good qualitative agreementto the data for odd-neutron nuclei. The agreement was significantly less good for the odd-protonnuclei, possibly due to the use of the Slater approximation for the exchange Coulomb interaction.The deformation energies of two odd-neutron nuclei (235U and 239Pu) have been calculated forsome single-particle configurations up to a point beyond the outer fission-barrier. Axial symmetrynuclear shape has been assumed while a breaking of the left-right (or intrinsic parity) symmetryhas been allowed around the outer fission-barrier. The fission-barrier heights of such odd-neutronnuclei depend significantly on the particle configurations. A special attention has been paid tothe very important rotational correction to deformation energies. In particular, the correction ofthe moment of inertia calculated from the usual Belyaev expression was considered. Overall, aqualitative agreement with available data on fission-barrier heights for the considered odd-neutronnuclei and their even neighbours has been obtained.

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.