Magnons topologiques dans les matériaux quantiques
La topologie est devenue un paradigme essentiel en matière condensée, permettant de classer les phases de la matière selon des propriétés invariantes sous des déformations continues. Les premières recherches dans ce domaine se sont essentiellement concentrées sur les structures de bandes électroniques, conduisant par exemple à la découverte des « isolants topologiques ». Cependant, ces concepts topologiques ne sont pas restreints seulement aux électrons (fermions) et ainsi, l'application de tels concepts aux bosons, en particulier les magnons, suscite un intérêt croissant. Les magnons, qui sont des excitations collectives dans les matériaux magnétiques, illustrent comment la topologie influence la dynamique magnétique et affecte le transport de chaleur et de spin. Des analogues d'isolants topologiques et de semi-métaux apparaissent dans leurs structures de bandes de magnons. Les magnons offrent ainsi une plateforme pour étudier l'interaction entre symétries magnétiques et topologie, examiner l'effet des interactions sur les bandes topologiques, et générer des courants de spin protégés aux interfaces. La recherche de matériaux contenant des magnons topologiques est donc cruciale, surtout pour les applications en magnonique, qui exploitent les ondes de spin pour le stockage et le traitement rapide des données.
Ce projet de thèse se consacre à explorer ces aspects topologiques dans des matériaux quantiques candidats à l’aide de techniques de diffusion de neutrons et de rayons X dans les grandes infrastructures de recherche (ILL, ESRF, SOLEIL), pour analyser la structure de bande des magnons à la recherche de caractéristiques topologiques, comme les points de Dirac ou de Weyl. Les résultats expérimentaux seront soutenus par des calculs théoriques des bandes magnoniques intégrant des concepts topologiques.
Imagerie des champs de déformations dans les semi-conducteurs: du matériau au dispositif
Ce sujet traite de la visualisation et de la quantification des champs de déformation dans les matériaux semi-conducteurs, par des techniques utilisant le rayonnement synchrotron. Le contrôle de la déformation est fondamental pour optimiser les propriétés de transport électronique, mécaniques et thermiques. Dans une approche duale, nous combinerons la détermination du tenseur local de déformation déviatorique en balayant l'échantillon sous un nano faisceau polychromatique (µLaue) et une imagerie d'un champ de vu donné (microscopie aux rayons X en champ sombre, DFXM).
Des recherches originales s’intéresseront à améliorer l’analyse : (1) de la précision et de la vitesse de détermination quantitative des champs de déformation, (2) des distributions des gradients de déformation, et (3) du champ de déformation dynamique dans les matériaux piézoélectriques par des mesures stroboscopiques. Pour illustrer ces points, trois cas scientifiques correspondant à des matériaux microélectroniques pertinents et de complexité croissante seront étudiés :
1.Champs de déformation statiques entourant des contacts métalliques dans le Si, tels que les vias à travers le silicium (TSV) à haute densité dans la technologie CMOS.
2.Gradients de déformation dans des structures hétéroépitaxiales complexes Ge/GeSn avec des variations de composition le long de la direction de croissance.
3.Études de déformation dynamique de résonateurs acoustiques LiNbO3 en volume avec une fréquence de résonance dans la plage des MHz.
La validation de cette approche conceptuelle permettra une avancée significative dans le domaine de la microélectronique et l'ingénierie de déformation.
Etude de la diffusion de petits amas de xénon au sein du combustible nucléaire métallique UMo
Ce sujet de thèse est centré sur l’application de méthodes de calcul à l’échelle des atomes afin d’étudier la diffusion et la stabilité intra-granulaire d’amas de Xe au sein du combustible métallique UMo.
Les alliages d’uranium-molybdène UMo présentent d’excellentes propriétés thermiques et une bonne densité en uranium. C’est notamment pour ces propriétés que l’UMo monolithique est considéré comme l’un des potentiels combustibles candidats pour les réacteurs de recherche. Il est donc crucial pour le CEA de développer de nouveaux modèles de calcul permettant d’analyser l’évolution des propriétés thermomécaniques de l’UMo en conditions d’irradiation.
Au cours de cette thèse, votre travail consistera dans un premier temps à valider ou recalibrer si nécessaire les modèles de calcul à l’échelle atomique existants pour l’UMo dans la littérature. Vous devrez ensuite de les appliquer à la simulation de la stabilité et de la diffusion de petits amas de xénon au sein de cristaux d’UMo. Ces calculs seront effectués à l’aide de méthodes de dynamique moléculaire accélérée novatrices, et seront systématiquement comparés aux résultats obtenus pour le combustible nucléaire de référence UO2. Après avoir analysé vos résultats par comparaison aux mesures expérimentales de collaborateurs du département, vous serez en charge de transférer les données produites à d’autres chercheurs du département afin d’alimenter les codes de simulation des combustibles nucléaires à plus grande échelle. Vos résultats seront publiés au sein de publications scientifiques, et vous présenterez vos résultats dans le cadre de conférences scientifiques.
L’ensemble de ces travaux vous permettrons de compléter votre formation en acquérant des compétences applicables à de nombreux domaines de la science des matériaux: calculs ab initio, ajustement de potentiels interatomiques par techniques de « machine learning », dynamique moléculaire classique et accélérée, utilisation des super-calculateurs du CEA, ainsi que de nombreux éléments de physique statistique et de physique de la matière condensée, méthodes dont les membres de l’équipe encadrante sont des spécialistes.
Vous serez accueilli au sein du Laboratoire de Modélisation du Comportement des Combustibles (Institut IRESNE, CEA Cadarache). Il s’agit d’un groupe de recherche dynamique, au sein duquel vous serez amené à collaborer avec les autres doctorants présents au laboratoire. L’environnement de travail sera de plus riche en collaboration nationales et internationales (expérimentateurs du département, Institut ISAS (CEA Saclay), Laboratoire CINAM à Marseille, collaborations avec les laboratoires nationaux américains), qui vous permettront de vous insérer au sein de la communauté de la recherche en matériaux pour les sciences du nucléaire.
Prédire les propriétés thermodynamiques des défauts dans des alliages métalliques multi-composants à partir de l'échelle atomique par apprentissage statistique
Les propriétés et le comportement des matériaux dans des conditions extrêmes sont essentiels pour les systèmes énergétiques tels que les réacteurs de fission et de fusion. Cependant, prédire avec précision les propriétés des matériaux à haute température reste un défi. Les mesures directes de ces propriétés sont limitées par les instruments expérimentaux, et les simulations à l'échelle atomique basées sur des champs de force empiriques sont souvent peu fiables en raison d'un manque de précision. Ce problème peut être résolu à l'aide de techniques d'apprentissage statistique, qui ont récemment vu leur utilisation exploser en science des matériaux. Les champs de force construits par apprentissage statistique atteignent le degré de précision des calculs {it ab initio} ; cependant, leur mise en œuvre dans les méthodes d'échantillonnage est limitée par des coûts de calcul élevés, généralement supérieurs de plusieurs ordres de grandeur à ceux des champs de force traditionnels. Pour surmonter cette limitation, deux objectifs seront poursuivis dans cette thèse : (i) améliorer les champs de force par apprentissage statistique actif en trouvant un meilleur compromis précision-efficacité et (ii) créer des méthodes accélérées d'échantillonnage de l'énergie libre et des chemins cinétiques afin de faciliter l'utilisation de champs de force d'apprentissage statistique coûteux en termes de calcul. Pour le premier objectif, nous améliorons la construction des champs de force d'apprentissage statistique en nous concentrant sur trois facteurs clés : la base de données, le descripteur de l'environnement atomique local et le modèle de régression. Pour le deuxième objectif, nous mettrons en œuvre un schéma d'échantillonnage bayésien rapide et robuste pour estimer l'énergie libre anharmonique, qui est cruciale pour comprendre les effets de la température sur les solides cristallins, à l'aide d'une méthode de force de biais adaptative qui améliore considérablement la vitesse de convergence et la précision globale. Nous appliquerons les méthodes développées au calcul de l'énergie libre et de ses dérivées, des grandeurs physiques qui donnent accès aux propriétés thermo-élastiques des alliages et aux propriétés thermodynamiques des défauts ponctuels. Pour cela, nous utiliserons des extension algorithmiques qui permettent d'échantillonner un état métastable spécifique et aussi les chemins de transition vers d'autres bassins d'énergie et donc d'estimer les énergies libres de formation et de migration de défauts lacunaires. Les grandeurs thermodynamiques calculées seront ensuite utilisées comme données d'entrée de méthodes de Monte Carlo cinétique, qui permettra de mesurer les coefficients de diffusion dans les alliages complexes en fonction de la température. Un but sera d'essayer de relier les propriétés de transport atomique à la complexité de l'alliage. Notre approche ayant une rapidité considérablement supérieure à celle des méthodes standard, nous pourrons envisager de l'appliquer à des alliages complexes comprenant les éléments W, Ti, V, Mo et Ta à des températures et des compositions qui n'ont pas été étudiées expérimentalement.
Comportement de nanocavités sous chargement mécanique : de la compréhension des mécanismes physiques à l’homogénéisation de matériaux nanoporeux
Des nanocavités - typiquement de quelques nm à quelques dizaines de nm - sont souvent observées dans les métaux, par exemple dans les applications hautes températures suite à la condensation de lacunes ou dans les alliages métalliques utilisés dans les réacteurs nucléaires du fait de l’irradiation. La présence de ces nanocavités dégrade le comportement mécanique des matériaux et contribue à la rupture. Il est donc nécessaire de déterminer les mécanismes physiques associés au comportement de ces nanocavités sous chargement mécanique et d’obtenir des modèles homogénéisés décrivant le comportement macroscopique des matériaux nanoporeux. Les résultats disponibles dans la littérature restent à ce jour limités, notamment en ce qui concerne la représentativité des simulations réalisées et des modèles proposés vis-à-vis des applications d’intérêt. Cela inclut par exemple la prise en compte des défauts cristallins entourant les cavités, l’effet des chargements cycliques et la localisation des nanocavités au niveau des joints de grains. Les objectifs de cette thèse sont donc de déterminer le comportement de nanocavités sous chargement mécanique et les mécanismes physiques associés en considérant des situations réalistes vis-à-vis des applications, de développer des modèles analytiques à bases physiques permettant de décrire le comportement de nanocavités sous chargement mécanique, et enfin de proposer des modèles homogénéisés adaptés aux nanocavités et utilisables pour simuler la rupture par croissance et coalescence de cavités. Les cas d’applications visés sont ceux liés aux alliages métalliques sous irradiation mais les éléments de compréhension obtenus et les modèles développés pourront être utilisés dans un contexte plus large. Afin d’atteindre ces objectifs, des simulations de Dynamique Moléculaire (DM) seront réalisées, analysées à partir de la théorie élastique des dislocations et utilisées pour proposer des modèles homogénéisés pertinents pour les matériaux nanoporeux.
Modélisation multi-échelles de la diffusion de l’hydrogène dans un polycristal de nickel
Dans de nombreuses applications, des matériaux de structure métalliques sont en contact avec de l’hydrogène qui va pénétrer dans le métal, dégrader ses propriétés mécaniques, et parfois conduire à la rupture du matériau. Les mécanismes de fragilisation par l’hydrogène ont été très étudiés par le passé. Néanmoins, il n’existe toujours pas de modèle prédictif quantitatif de ces phénomènes. Cette thèse s’intéresse à la ségrégation de l’hydrogène aux joints de grains qui est un des mécanismes de fragilisation observé. L’objectif est de modéliser la cinétique de ségrégation en partant de l’échelle atomique, ce qui implique donc de trouver les structures d’équilibre des joints de grains, d’identifier les sites de ségrégation pour chaque joint de grain, et de quantifier l’effet d’un joint de grain sur le coefficient de diffusion de l’hydrogène. Ces données alimenteront ensuite un modèle en élément finis qui permettra de calculer la répartition de l’hydrogène au cours du temps, en prenant en compte la microstructure polycristalline de l’échantillon et les propriétés spécifiques à chaque joint de grain. Ces résultats seront comparés à des expériences de perméation d’hydrogène qui donnent accès à un coefficient de diffusion moyen, ainsi qu’à des mesures localisées sur un joint de grain particulier (méthodes PANI et SKPFM)
Etude à l’échelle atomique de la mobilité des dislocations dans le combustible MOX
La transition vers la neutralité carbone exige une augmentation rapide des énergies décarbonées, dont le nucléaire, qui nécessite une compréhension approfondie des matériaux irradiés. Le combustible à oxyde mixte (MOX) est particulièrement important, car il optimise l'utilisation des ressources nucléaires et réduit les déchets radioactifs. Le comportement mécanique du MOX sous irradiation est crucial pour garantir l’intégrité du combustible dans diverses conditions de fonctionnement.
L’objectif de la thèse est de réaliser des simulations atomistiques afin de comprendre la mobilité des dislocations, essentielle pour soutenir la modélisation multiéchelle du comportement mécanique du MOX. Des calculs de dynamique moléculaire permettront d'analyser la mobilité des dislocations en fonction de diverses conditions de température, de contraintes, de teneur en plutonium et de déviations à la stœchiométrie, avec pour but d’établir des lois de vitesse. Les résultats de ces simulations amélioreront la modélisation micromécanique dans la plateforme de simulation PLEIADES du CEA, dédiée à la simulation du cycle de vie complet du combustible nucléaire, de sa fabrication jusqu'à l’entreposage.
Le doctorant sera accueilli au Laboratoire de Modélisation du Comportement des Combustibles (Institut IRESNE, CEA-Cadarache), un environnement dynamique composé de 11 chercheurs et d'un nombre équivalent de doctorants. Situé en Provence, ce centre offre un cadre de travail agréable, entre les parcs naturels du Verdon et du Lubéron. La thèse se fera en collaboration avec l'IM2NP, un laboratoire à la pointe de la recherche en physique des matériaux.
Le candidat doit avoir de solides bases en physique des matériaux, idéalement en mécanique aux petites échelles. Ces compétences pourront être renforcées durant un stage de M2 au sein du laboratoire. Le doctorant valorisera son travail à travers des publications scientifiques et des présentations en conférences internationales, ouvrant ainsi des opportunités dans les domaines de la recherche et de l'industrie.
Structure et mobilité des agrégats et boucles d'interstitiels dans l'oxyde d'uranium
L’oxyde d’uranium (UO2) est le combustible usuel des centrales nucléaires à fission. A ce titre son comportement sous irradiation est très étudié. L’irradiation crée des défauts lacunaires ou interstitiels qui vont piloter l’évolution de la microstructure du matériau qui elle-même va impacter ses propriétés physiques (par exemple sa conductivité thermique) et mécaniques. Les agrégats d’interstitiels en particulier jouent un rôle prépondérant.
D’une part, aux plus petites tailles, la diffusion des interstitiels dans UO2 est encore assez mal comprise. En effet, expérimentalement, on observe l’apparition de boucles de dislocations constituées d’interstitiels de tailles pouvant atteindre la dizaine de nanomètres. A l’inverse on n’observe pas de cavités et les défauts lacunaires restent de tailles subnanométriques. Cela dénote une diffusion plus rapide des interstitiels que des lacunes, la diffusion permettant l’agglomération des interstitiels et la formation de boucles. Pourtant les calculs à l’échelle atomique ne montrent pas de différence majeure entre les coefficients de diffusion des lacunes et des interstitiels dans UO2. Une hypothèse pour expliquer cette contradiction apparente est que ce seraient les agrégats d’interstitiels qui diffuseraient rapidement (Garmon, Liu et al. 2023).
D’autre part, on s’attend à ce que les agrégats d’interstitiels tridimensionnels soient les germes des boucles de dislocations observées en microscopie électronique à transmission dans l’oxyde d’uranium irradié. Mais les mécanismes de transformations des agrégats en boucles et de changement de nature de boucles restent incompris dans l’oxyde d’uranium. Ces mécanismes ont très récemment été élucidée pour des métaux cubique à face centré (Jourdan, Goryaeva et al. 2024). Il est possible que des mécanismes comparables soient à l’œuvre dans UO2 avec la complication induite par l’existence deux sous-réseaux.
Nous proposons donc d’étudier par simulations à l’échelle atomique les agrégats d’interstitiels dans UO2.
On abordera d’abord la structure de ces agrégats subnanométriques tridimensionnels. Pour ce faire nous utiliserons les outils de classification des structures de défauts par intelligence artificielle mises au point au laboratoire (Goryaeva, Lapointe et al. 2020). On étudiera la diffusion de ces objets, par dynamique moléculaire et par recherche automatique de cols de migration à l’aide d’outils de type kinetic-ART (Béland, Brommer et al. 2011). Dans un deuxième temps, on étudiera la stabilité relative des agrégats 3D et des boucles de dislocations fautées et parfaites et les transformations entre ces différents objets.
Cette étude devra reposer sur des potentiels d’interaction interatomiques. On commencera par utiliser les potentiels empiriques disponibles dans la littérature avant nous tourner vers les potentiels de type Machine Learning (Dubois, Tranchida et al. 2024) en développement au Département d’Etudes du Combustibles du CEA Cadarache.
Béland, L. K., et al. (2011). "Kinetic activation-relaxation technique." Physical Review E 84(4): 046704.
Chartier, A., et al. (2016). "Early stages of irradiation induced dislocations in urania." Applied Physics Letters 109(18).
Dubois, E. T., et al. (2024). "Atomistic simulations of nuclear fuel UO2 with machine learning interatomic potentials." Physical Review Materials 8(2).
Garmon, A., et al. (2023). "Diffusion of small anti-Schottky clusters in UO2." Journal of Nuclear Materials 585: 154630.
Goryaeva, A. M., et al. (2020). "Reinforcing materials modelling by encoding the structures of defects in crystalline solids into distortion scores." Nature Communications 11(1).
Jourdan, T., et al. (2024). "Preferential Nucleation of Dislocation Loops under Stress Explained by A15 Frank-Kasper Nanophases in Aluminum." Physical Review Letters 132(22).
Compréhension des mécanismes de piégeage de l’hélium dans des nouvelles nuances d’alliages base nickel développées pour les réacteurs à sel fondu
Les alliages base nickel sont les matériaux de structure naturels envisagés pour les Réacteurs à Sel Fondu (RSF). Ils présentent d’excellentes caractéristiques mécaniques et une bonne résistance à la corrosion. Dans ces matériaux, la production d’hélium, principalement causée par la transmutation du nickel par les neutrons rapides peut atteindre des teneurs suffisantes pour fragiliser fortement le matériau ou provoquer son gonflement sous irradiation. L’hélium est très peu soluble dans le matériau et condense sous forme de bulles ou ségrége aux joints de grains. Pour limiter ces phénomènes et réussir à piéger l’hélium, une solution consiste à introduire dans le matériau qui sera irradié, une densité importante de nano-précipités dont les interfaces serviront de site de germination pour des bulles nanométriques à même de piéger l’hélium « sur place » pour empêcher ce dernier de migrer vers les joints de grains et de dégrader les performances du matériau. Il s'agira d'étudier par microscopie électronique en transmission corrigée la cinétique de précipitation des phases thermodynamiquement attendue ainsi que la structure atomique des interfaces formées entre les précipités et la matrice. Une simulation de la précipitation par champ de phase sera également envisagée. Enfin, les mécanismes de piégeage de l'He aux interfaces sera étudié à l'aide de la spectroscopie de perte d'énergie des électrons (EELS).
Simulation atomistique de la rupture de verres borosilicatés hétérogènes
Les verres borosilicatés hétérogènes contiennent des précipités cristallins ou amorphes qui forment des phases secondaires incrustées dans la matrice vitreuse. Ces matériaux sont appréciés pour leur résistance élevée au choc thermique et leur excellente durabilité chimique, les rendant idéaux pour diverses applications telles que les ustensiles de cuisine et le matériel de laboratoire. En particulier, dans l'industrie nucléaire, de nombreuses matrices vitreuses de conditionnement de déchets radioactifs contiennent des précipités en raison de la présence d'éléments peu solubles.
Il a été démontré que des phases secondaires peuvent affecter considérablement les propriétés mécaniques, en particulier la résistance à la fracture. Cependant, les mécanismes spécifiques liés à ce phénomène à l'échelle atomique restent mal expliqués. En particulier, il est crucial de comprendre l'effet de la nature de ces phases (cristallines ou amorphes) et de leur interface avec la matrice vitreuse.
L'objectif principal de ce projet est d'étudier les mécanismes spécifiques par lesquels les précipités influencent les propriétés mécaniques à l'échelle atomique. Il vise également à comprendre comment ces précipités affectent la propagation de fissures. Pour cela, des outils de modélisation numérique basés sur la dynamique moléculaire seront utilisés. Cette technique simule le comportement individuel des atomes au fil du temps sous différentes conditions de test. Ainsi, elle permet d'examiner la structure locale des pointes de fissure et leur interaction avec les précipités à l'échelle atomique, fournissant des informations précieuses sur les mécanismes sous-jacents de résistance aux fissures dans les verres hétérogènes.