Nouvelles contraintes expérimentales sur les constantes de couplage de l’interaction faible par la mesure en coïncidence de schémas de désintégration complexes

La caractérisation expérimentale précise des transitions bêta interdites non-uniques, représentant environ un tiers de toutes les transitions bêta connues, est un sujet à la fois important et très ardu. De fait, très peu d’études fiables existent dans la littérature. En effet, le spectre en énergie continu de ces transitions est difficile à mesurer précisément pour diverses raisons qui se cumulent les unes aux autres : grande diffusivité des électrons dans la matière et non-linéarité du système de détection, indisponibilité de certains radionucléides et présence d’impuretés, longues périodes de désintégration et complexité des schémas, etc. Des prédictions théoriques réalistes sont tout aussi difficiles car il est nécessaire de coupler des modélisations précises des structures atomiques et nucléaires des radionucléides à travers l’interaction faible, dans un même formalisme complètement relativiste. Pourtant, améliorer notre connaissance des transitions bêta interdites non-uniques est essentiel en métrologie de la radioactivité pour définir l’unité SI du becquerel dans le cas des émetteurs bêta purs. Cela peut avoir un impact fort en médecine nucléaire, pour l’industrie du nucléaire, et pour certaines thématiques de physique fondamentale, comme la recherche de matière noire et la physique des neutrinos de réacteurs.
Notre étude récente de la transition deuxième interdite non-unique du 99Tc, à la fois théorique et expérimentale, a mis en évidence que les transitions interdites non-uniques peuvent être particulièrement sensibles à la valeur effective des constantes de couplage de l’interaction faible. Ces dernières interviennent comme facteurs multiplicatifs des éléments de matrice nucléaires. L’utilisation de valeurs effectives permet de compenser les approximations employées dans les modèles de structure, telles que des corrélations simplifiées entre les nucléons dans l’espace de valence ou l’absence d’excitation du cœur. Cependant, leur ajustement ne peut se faire que par comparaison avec un spectre expérimental de grande précision. La prédictibilité des calculs théoriques, même les plus précis actuellement disponibles, est ainsi fortement remise en cause. S’il a déjà été démontré que des valeurs universelles ne peuvent être fixées, des valeurs pour chaque type de transition, ou pour un modèle nucléaire spécifique, sont possibles. Le but de ce sujet de thèse est donc d’établir de nouvelles contraintes expérimentales sur les constantes de couplage de l’interaction faible en mesurant précisément les spectres en énergie de transitions bêta. À terme, cela permettra d’établir des valeurs effectives moyennes robustes de ces constantes de couplage et d’obtenir un vrai pouvoir prédictif pour les calculs théoriques de désintégration bêta.
La plupart des transitions d’intérêt pour contraindre les constantes de couplage ont des énergies supérieures à 1 MeV et se situent au sein de schémas de désintégration complexes, avec émission de multiples photons gamma. Dans cette situation, la meilleure stratégie consiste en une détection bêta-gamma en coïncidence. Les techniques usuelles de détection en physique nucléaire sont appropriées mais nécessitent d’être extrêmement bien implémentées et contrôlées. Le doctorant pourra s’appuyer sur les résultats obtenus lors de deux thèses précédentes. Pour minimiser le phénomène d’auto-absorption des électrons dans la source, il devra améliorer une technique de préparation de sources radioactives ultra-minces développée au LNHB pour l’adapter aux activités importantes qui seront nécessaires. Il devra implémenter un nouveau dispositif de mesure, dans une chambre à vide dédiée, comprenant une détection en coïncidence de deux détecteurs silicium et deux détecteurs gamma. Plusieurs études seront nécessaires, mécaniques et par simulation Monte Carlo, pour optimiser la configuration géométrique en regard des différentes contraintes. L’optimisation de la chaîne d’électronique, l’acquisition, le traitement du signal, l’analyse des données, la déconvolution spectrale et l’élaboration d’un bilan d’incertitudes complet et robuste seront autant de sujets abordés. Ces développements instrumentaux permettront de mesurer avec une grande précision les spectres du 36Cl, du 59Fe, du 87Rb, du 141Ce, ou encore du 170Tm.
Ce sujet très complet permettra au doctorant d’acquérir des compétences instrumentales et d’analyse qui lui ouvriront de nombreuses opportunités de carrière. Le candidat devra posséder de bonnes connaissances en instrumentation nucléaire, en programmation et en simulations Monte Carlo, ainsi qu’une connaissance raisonnable des désintégrations nucléaires.

Simulation quantique des noyaux atomiques

Les noyaux atomiques forment des systèmes quantiques à N corps fortement corrélés régis par l'interaction forte de la QCD. Le modèle en couches nucléaire, qui diagonalise le hamiltonien dans une base dont la dimension croît exponentiellement avec le nombre de nucléons, constitue une approche bien établie pour décrire leur structure. Cependant, cette explosion combinatoire limite les calculs classiques haute performance à une fraction restreinte de la carte des noyaux.
Les calculateurs quantiques offrent une alternative prometteuse grâce à leur capacité naturelle à manipuler des espaces de Hilbert exponentiellement grands. Bien que nous demeurions dans l'ère NISQ avec ses qubits bruités, ils pourraient révolutionner les applications du modèle en couches.
Cette thèse vise à développer une approche complète pour la simulation quantique des systèmes nucléaires complexes. Un premier jalon crucial consiste à créer une interface logicielle intégrant les données de structure nucléaire (orbitales nucléoniques, interactions nucléaires) aux plateformes de calcul quantique, facilitant ainsi les applications futures en physique nucléaire.
Le projet explore deux classes d'algorithmes : les approches variationnelles et non-variationnelles. Pour les premières, l'expressivité des ansätze quantiques sera analysée systématiquement, notamment dans le contexte de la brisure et de la restauration de symétries. Les solveurs quantiques variationnels (VQE), particulièrement prometteurs pour les systèmes hamiltoniens, seront implémentés avec un accent sur la technique ADAPT-VQE adaptée au problème à N corps nucléaire.
Un défi majeur réside dans l'accès aux états excités, aussi cruciaux que l'état fondamental en structure nucléaire, alors que les VQE se concentrent principalement sur ce dernier. La thèse développera donc des algorithmes quantiques dédiés aux états excités, en testant diverses méthodes : expansion de l'espace de Hilbert (Quantum Krylov), techniques de fonctions de réponse (équations du mouvement quantiques), et méthodes basées sur l'estimation de phase. L'objectif final est d'identifier les approches les plus adaptées en termes de passage à l'échelle et de résistance au bruit, pour des applications avec des hamiltoniens nucléaires réalistes.

Description de phénomènes collectifs au sein des noyaux d’atome au-delà de l’approche de la fonctionnelle de la densité dépendant du temps

Contexte :
Prédire l'organisation et la dynamique des neutrons et protons au sein des noyaux d'atome est un véritable
défi scientifique, crucial pour concevoir les technologiques nucléaires du futur mais aussi pour répondre à
des questions fondamentales comme l'origine des atomes lourds dans notre univers. Dans ce cadre, le CEA,
DAM, DIF développe des approches théoriques visant à simuler la dynamique des constituants élémentaires
du noyau atomique. Les équations du mouvement obtenues, dans le cadre de la mécanique quantique, sont
résolues sur nos supercalculateurs. Les années 2010 ont vu une montée en puissance de l’approche dite de la
fonctionnelle de la densité dépendant du temps (TDDFT) pour traiter ce problème. Malgré sa description
révolutionnaire de certains phénomènes tels que les résonances géantes observées dans les noyaux d’atome
ou encore la fission nucléaire, cette approximation possède des limites intrinsèques.

Objectifs :
Cette thèse vise à développer et explorer une nouvelle approche théorique permettant de décrire les
phénomènes de mouvements collectifs des protons et neutrons au sein du noyau atomique. L’idée est de
généraliser l’approche TDDFT afin de mieux prédire certaines propriétés des réactions nucléaires comme la
répartition de l’énergie entre les fragments issus d’une fission nucléaire. Partant de premier travaux allant
dans cette direction, le/la doctorante devra dériver les équations du mouvement de cette nouvelle approche
puis les implémenter sous forme d’une librairie C++ optimisée et tirant profit des supercalculateurs du CEA.
Le but sera enfin d’étudier comment cette nouvelle approche améliore les prédictions de certains
phénomènes comme l’amortissement des résonances géantes dans les noyaux d’atomes ou encore la
naissance des fragments générés lors d’une fission nucléaire.

Description microscopique des propriétés des fragments de fission évaluées à la scission

La fission est l’une des réactions nucléaires les plus difficiles à décrire, reflétant la diversité des aspects dynamiques du problème à N-corps. Au cours de ce processus, le noyau explore des états de déformations extrêmes aboutissant à la formation de deux fragments. Alors que la nombre de degrés de liberté (DLC) mis en jeu est extrêmement grand, l’approximation de champ moyen est un bon point de départ qui opère une réduction drastique des DLC, l’élongation et l’asymétrie étant incontournables. Cette réduction introduit des discontinuités dans la génération successive des états par lesquels le noyau transite, la continuité en énergie n’assurant pas la continuité des états issus d’un principe variationnel. Récemment, une nouvelle méthode basée sur des contraintes associées aux recouvrements de fonctions d’onde, a été mise en oeuvre afin d’assurer cette continuité jusqu’à la scission et au-delà (vallée de Coulomb). Cette continuité est capitale pour décrire la dynamique du processus.
L’objectif de la thèse que nous proposons est de réaliser pour la première fois la mise en oeuvre à deux dimensions de cette nouvelle approche afin de prendre en compte l’ensemble de la collectivité générée par les DLC d’élongation et d’asymétrie. Les développements théoriques et numériques s’inscrivent dans le cadre de la méthode de la coordonnée génératice dépendant du temps. Ce type d’approche contient une première étape statique qui consiste en la génération de surfaces d’énergie potentielle (PES) obtenues par des calculs Hartree-Fock-Bogoliubov sous contraintes et une seconde étape dynamique qui décrit la propagation dynamique d’un paquet d’onde sur ces surfaces via la résolution de l’équation de Schrödinger dépendant du temps. C’est à partir de cette deuxième étape que les observables sont en général extraites.
Dans le cadre de cette thèse, le doctorant devra:
- en première étape, construire des PES à deux dimensions continues pour l’état adiabatique et les premiers états excités. Cela mettra en oeuvre les trois algorithmes Link, drop et Deflation
- en deuxième étape, extraire des observables accessibles par ce type d’approches: les rendements, le bilan énergétique à la scission, la déformation des fragments, le nombre moyen de neutrons émis. Nous souhaitons en particulier étudier l’impact des excitations intrinsèques sur les observables de fission, qui se manifestent essentiellement dans la descente du point scelle vers la scission.
Enfin, ces résulats seront confrontés aux données expérimentales, dans des actinides et des pré-actinides d’intérêt. Notamment, les mesures récentes très précises obtenues par les expériences SOFIA pour des noyaux peu à très exotiques devraient contribuer à tester la précision et la prédictivité de nos approches, et à guider les futurs développements des approches à N-corps et de l’interaction nucléaire dans le cadre de la fission.

Vers une plateforme d’irradiation photonique multimodale : fondements et conceptualisation

Les techniques d’irradiation photonique exploitent les interactions entre un faisceau de photons de haute énergie et la matière pour effectuer des mesures non destructives. En induisant des réactions photonucléaires, telles que l’activation photonique, les résonances de fluorescence nucléaire (NRF) et la photofission, ces techniques d’irradiation permettent de sonder la matière en profondeur. L’association de ces différentes techniques de mesure nucléaire au sein d’une plateforme d’irradiation unique permettrait une identification précise et quantitative d’une grande variété d’éléments, en sondant le volume des matériaux ou objets étudiés. Le faisceau de photons de haute énergie est généralement produit par rayonnement de freinage (phénomène de Bremsstrahlung) au sein d’une cible de conversion d’un accélérateur linéaire d’électrons. Une alternative innovante consiste à exploiter les électrons de haute énergie délivrés par une source laser-plasma, convertis par rayonnement de freinage ou par diffusion Compton inverse. Une plateforme basée sur une telle source offrirait de nouvelles possibilités, car les sources laser-plasma peuvent atteindre des énergies significativement supérieures, permettant ainsi l'accès à de nouvelles techniques et applications d'imagerie avancées. L’objectif de cette thèse est d’établir les fondements et de conceptualiser une plateforme d’irradiation photonique multimodale. Un tel dispositif viserait à se baser sur une source laser-plasma et permettrait la combinaison des techniques d’activation photonique, des résonances de fluorescence nucléaire (NRF) et de la photofission. En repoussant les limites des mesures nucléaires non destructives, cette plateforme offrirait des solutions innovantes à des défis majeurs dans des secteurs stratégiques tels que la sécurité et le contrôle aux frontières, la gestion des colis de déchets radioactifs, ainsi que l'industrie du recyclage.

Réactions nucléaires induites par des anti-ions légers – apport du modèle INCL

L’interaction d’une antiparticule avec un noyau atomique est un type de réaction qu’il faut savoir simuler pour pouvoir répondre à des questions fondamentales. On peut citer comme exemples, la collaboration PANDA (FAIR) avec des faisceaux d’antiproton de l’ordre du GeV qui envisage l’étude des interactions nucléon-hypéron, ainsi que celle de la peau de neutron, par la production d’hypérons et d’antihypérons. Cette même peau de neutron est aussi étudiée avec des antiprotons au repos avec l’expérience PUMA (AD - Cern). Au même endroit nous collaborons avec l’expérience ASACUSA pour l’étude de la production des particules chargées. Pour répondre à ces études, notre code de réactions nucléaires INCL a été étendu aux antiprotons (thèse D. Zharenov soutenue fin 2023). Au-delà de l’antiproton il y a les antideutérons et antiHe-3. Ces antiparticules sont d’un intérêt plus récent, avec notamment l'expérience GAPS (General AntiParticle Spectrometer) qui vise à mesurer les flux de ces particules dans le rayonnement cosmique. L’idée est de mettre en évidence la matière noire, dont ces particules seraient des produits de décroissance, et dont la quantité mesurée doit ressortir plus facilement du bruit de fond astrophysique que dans le cas des antiprotons. Le sujet proposé est donc l’implantation des anti-noyaux légers dans INCL avec comparaisons à des données expérimentales.

Mesure du flux elliptique des quarks charmés dans les collisions Pb-Pb semi-centrales à 5 TeV au CERN avec LHCb.

Les collisions d'ions lourds offrent une opportunité unique d'étudier le plasma de quarks et de gluons (QGP), un état exotique de la matière dans lequel les quarks et les gluons ne sont plus confinés dans les hadrons, et qui aurait existé quelques microsecondes après le Big Bang. Parmi les sondes clés pour l'étude du QGP figurent les quarks charmés. En effet, ces derniers conservent l'histoire de leurs interactions avec le QGP, les rendant essentiels pour comprendre les propriétés du QGP. La production de quarks charmés et leurs interactions avec le QGP sont étudiées à travers les mesures des hadrons, mésons et baryons contenant au moins un quark ou antiquark charm, tels que les mésons D0 ou les baryons Lambda_c. Cependant, le processus d'hadronisation — la manière dont les quarks charm se confinent dans des baryons ou mésons incolores — reste encore mal compris.

Une approche prometteuse pour approfondir la compréhension de l'hadronisation des quarks charmés réside dans la mesure de leur écoulement elliptique, une mesure de corrélations angulaires à longue distance, une signature des effets collectifs dus à la thermalisation du QGP. En comparant l'écoulement elliptique des mésons D0 et des baryons Lambda_c, sensible aux propriétés du milieu créé, les chercheurs peuvent approfondir leurs connaissance sur le mécanisme d'hadronisation des quarks charmés.

Pour mesurer cet écoulement elliptique, l'étudiant.e sélectionné.e développera une méthode innovante exploitant pleinement les capacités du détecteur LHCb. Cette méthode, jamais appliquée auparavant, permet une interprétation plus intuitive et théoriquement robustes des mesures d'écoulement elliptique par rapport aux méthodes traditionnelles. Le/la candidat.e adaptera cette technique pour le détecteur LHCb afin de mesurer, comparer et interpréter l'écoulement elliptique des baryons charmés Lambda_c et des mésons D0 dans les nouvelles données PbPb collectés par LHCb en 2024.

Etude des mécanismes de réaction pour la synthèse d’éléments super-lourds

Cette thèse a pour but d’étudier les mécanismes de réaction menant à la synthèse de noyaux super-lourds en appui au programme de recherche expérimental qui sera développé auprès de l’installation S3 de Spiral2 au GANIL à Caen. Elle vise à renforcer la précision de la modélisation de la réaction en évaluant, grâce à l’analyse d’incertitudes, les expériences les plus prometteuses pour contraindre les paramètres et la modélisation.

L'une des principales activités de la physique nucléaire est l'étude des propriétés des noyaux exotiques jusqu'aux limites d’existence des noyaux, dans les régions où les rapports proton-neutron sont extrêmes (driplines proton/neutron) et aux nombres de masse A et atomiques Z les plus élevés. Les noyaux dits super-lourds ne peuvent exister au-delà de la limite établie par le modèle de la goutte liquide – définie par une barrière de fission qui disparaît –, que grâce aux effets en couches de la mécanique quantique. Ces noyaux sont particulièrement intéressants parce qu'ils se situent à la limite entre la physique des petits systèmes à quelques corps et celle des systèmes à grand nombre de corps : les nombres magiques de protons et de neutrons, Z et N, sont remplacés par une région ou un îlot de magicité étendu en Z et N.

La synthèse de ces noyaux très- et super-lourds par des réactions de fusion-évaporation est un défi expérimental en raison des sections efficaces extrêmement faibles. La modélisation de la réaction complète afin de guider les expériences est également un défi difficile, car les modèles développés pour les noyaux plus légers ne peuvent pas être simplement extrapolés. Les réactions de fusion sont entravées par rapport à ce qui est observé avec les noyaux légers en raison de la très forte interaction Coulombienne, qui est renforcée par la forte répulsion causée par le grand nombre de charges positives (protons) dans le système, en concurrence avec la force d'attraction forte (nucléaire) dans un régime hautement dynamique. Le pouvoir prédictif des modèles doit être amélioré, bien que l'origine du phénomène d'entrave soit qualitativement bien comprise. Les ambiguïtés quantitatives sont suffisamment importantes pour observer des différences de quelques ordres de grandeur dans les probabilités de fusion calculées par différents modèles. Une petite modification de la section efficace pourrait nécessiter de nombreux mois pour réaliser des expériences réussies.

Au GANIL, en collaboration avec d'autres instituts, nous avons développé un modèle qui décrit les trois étapes de la réaction de synthèse des noyaux super-lourds. Les développements futurs se concentreront sur la recherche de moyens d'évaluer les modèles afin d'améliorer leur pouvoir prédictif, notamment en concevant des expériences spécifiques afin de contraindre l’amplitude de l'entrave à la fusion. Bien entendu, une analyse minutieuse de l'incertitude permettra d'améliorer le pouvoir prédictif des modèles. Des méthodes standard ainsi que des méthodes d'analyse de données de pointe telles que l'analyse bayésienne peuvent être utilisées.

Ce travail de doctorat sera effectué en collaboration avec le groupe expérimental du GANIL et une équipe de recherche à Varsovie (Pologne). En fonction des compétences de l'étudiant, la thèse sera plus orientée vers des développements formels ou vers les expériences menées dans la nouvelle installation S3 sur Spiral2. La participation aux expériences est possible.

Etudes expérimentales et théoriques de la génération du moment angulaire nucléaire et de l’énergie d’excitation des fragments de fission

La découverte de la fission en 1939 a profondément modifié notre connaissance de la physique nucléaire. Cette réaction permet de diviser des noyaux lourds comme l'uranium 235, en deux noyaux (fragments) plus légers, tout en libérant une grande quantité d'énergie. Les travaux de recherche sur la fission prennent la forme de modèles nucléaires servant à produire des bases de données nucléaires, qui sont essentiels pour simuler les réacteurs nucléaires. La qualité de ces données est encore insuffisante aujourd’hui, car notre compréhension fine de la fission reste très fragmentaire.
Ce travail de thèse vise à mieux décrire la génération du moment angulaire et l'énergie d'excitation des fragments de fission d’un point de vue expérimental et théorique. Ces recherches permettront à la fois de mieux comprendre le processus sous-jacent et d’améliorer le pouvoir de prédiction des outils de simulations, notamment les modèles utilisés pour calculer les échauffements gamma au sein d’un réacteur. Une partie du travail du doctorant consistera en l’exploitation des données acquises durant une thèse récente. Une autre partie sera la participation à des campagnes expérimentales complémentaires auprès du réacteur nucléaire de l’Institut Laue-Langevin (ILL), à l’aide du spectromètre LOHENGRIN afin de mesurer les rapports isomériques et les distributions en énergie cinétique des fragments de fission.
Le doctorant sera positionné au sein d’un laboratoire de physique nucléaire et de physique des réacteurs. Il développera des compétences en analyse de données, en physique nucléaire ainsi qu’en programmation informatique. Les langages utilisés seront C++ et python. Les débouchés sont la recherche en milieu académique ou industriel, également des postes de Data Scientist.

Test d’invariance de renversement du temps dans la désintégration beta nucléaire : analyse des données de MORA à JYFL

L’expérience MORA recherche des signes de violation de CP dans la désintégration beta d’ions piégés polarisés. Elle emploie des techniques de pointe afin d’atteindre une sensibilité jamais atteinte pour la mesure de la corrélation D (<10-4). Cette corrélation est sensible à une Nouvelle Physique qui pourrait expliquer l’asymétrie matière antimatière observée dans l’univers. La thèse consiste en l’analyse des données de la campagne qui se poursuit à Jyväskylä pour 23Mg+ et 39Ca+, en Finlande.

Top