Modelling of metal-hydrogen systems by DFT, CVM and Calphad method
Modélisation de systèmes métal-hydrogène par couplage des méthodes DFT, CVM et Calphad
Résumé
Hydrogen absorption in the interstitial sites of metals is crucial for major issues such as alloy embrittlement or hydrogen storage for energy applications. This phenomenon modifies the physicochemical properties of the host metal and may lead to the formation of ordered MHy compounds called hydrides. Within this framework, the Calphad modeling method (CALculation of PHAse Diagrams) is a relevant tool for understanding and predicting the behavior of metals and alloys in the presence of hydrogen. However, there is no Calphad database centered on hydrogen for calculating phase equilibria in multi-constituent systems (ternary, quaternary…).The present thesis proposes to use a multi-scale modeling approach to study metal-hydrogen (M-H) binary systems, which are the first step in designing such a Calphad database. First, systematic DFT (Density Functional Theory) calculations were carried out for 31 binary M-H systems considering 30 potential crystal structures, resulting in 30 × 31 = 930 hydrides, stable or metastable. This high throughput approach allowed in particular to determine the enthalpies of formation at 0 K, which represent important input data for the Calphad method. New hydrides that have never been experimentally observed could be predicted at high pressure (TaH2, ZrH3 ...).Then, phonon calculations in the harmonic approximation were performed on the most stable hydrides. They allow, on the one hand, to correct the DFT calculated enthalpies of formation by considering the energy and entropy due to the atom vibrations, which are not negligible for the light hydrogen atom. On the other hand, a large-scale study focused on the modification of the free energy of formation due to hydrogen substitution by its isotopes, known as "isotopic effect". Predictions were made on the nature of this effect as function of temperature. Moreover, the random insertion of hydrogen atoms in solid solution was studied using statistical thermodynamic methods: CVM (Cluster Variation Method) and Monte-Carlo simulation. These methods have been implemented in calculation codes, applied to face centered cubic (fcc) and body-centered cubic (bcc) metals. The input data are the interactions energies between nearest neighbor atoms. They are provided by the CEM (Cluster Expansion Method) coupled with DFT calculations. A comparative study of the Ni-H and Pd-H systems revealed the specificities of the thermodynamic behaviors of both solid interstitial solutions. Furthermore, dihydrogen pressure is an important parameter because many hydrides form only at very high pressure. To improve the Calphad model accuracy at high pressure, the model of Lu et al. was applied to the condensed phases of the Ni-H, Rh-H and Mg-H systems. This model allows to determine the contribution to the free enthalpy due to the pressure force work. The input data may be both quasi-harmonic phonon calculation results and experimental data. Finally, a comprehensive Calphad model of the Ni-H system was carried out by integrating the model of Lu et al. The DFT enthalpy of formation and the mixing enthalpy determined by CEM were used as input data, to complement the available experimental data
L'absorption d'hydrogène dans les sites interstitiels des métaux se situe au cœur de problématiques majeures comme la fragilisation des alliages ou le stockage de l'hydrogène à des fins énergétiques. En effet, ce phénomène modifie les propriétés physico-chimiques du métal hôte et peut conduire à la formation de composés ordonnés MHy appelés hydrures. Dans ce cadre, la méthode de modélisation Calphad (CALculation of PHAse Diagrams) constitue un outil pertinent pour comprendre et prédire le comportement des métaux et alliages en présence d'hydrogène. Toutefois, il n'existe pas de base de données Calphad centrée sur l'hydrogène, permettant de calculer des équilibres de phases dans des systèmes multi-constituants (ternaires, quaternaires…).Ainsi, la présente thèse est consacrée à l'étude de systèmes binaires métal-hydrogène (M-H), par une approche modélisatrice multi-échelle. Ces systèmes binaires représentent la première étape à la conception d'une telle base de données Calphad. Tout d’abord, des calculs DFT (Density Functional Theory) systématiques ont été réalisés pour 31 systèmes binaires M-H considérant 30 structures cristallographiques potentielles, soit 30 × 31 = 930 hydrures, stables ou métastables. Cette approche de criblage a permis notamment de déterminer les enthalpies de formation à 0 K, données d’entrée essentielles de la méthode Calphad. De nouveaux hydrures n'ayant jamais été observés expérimentalement ont été prédits à haute pression (TaH2, ZrH3…).Puis, des calculs de phonon dans l'approximation harmonique ont été réalisés sur les hydrures les plus stables. D’une part, ils permettent de corriger les enthalpies de formation issues de la DFT, en tenant compte de l'énergie et de l'entropie dues aux vibrations des atomes, non négligeables pour l'atome léger d'hydrogène. D’autre part, une étude à grande échelle a porté sur les modifications de l'enthalpie libre de formation résultant de la substitution de l'hydrogène par ses isotopes, dit « effet isotopique ». Des prédictions ont pu être réalisées sur la nature de cet effet en fonction de la température. De plus, pour étudier l'insertion aléatoire des atomes d'hydrogène en solution solide, nous avons utilisé des méthodes de thermodynamique statistique : CVM (Cluster Variation Method) et simulation de Monte-Carlo. Ces méthodes ont été implémentées dans des codes de calcul, appliqués aux métaux cubiques faces centrées (cfc) et cubiques centrés (cc). Les données d'entrée de ces codes sont les énergies associées aux interactions locales entre atomes plus proches voisins. Elles sont fournies par la CEM (Cluster Expansion Method) couplée à la DFT. Une étude comparée des systèmes Ni-H et Pd-H a mis en évidence les spécificités des comportements thermodynamiques des solutions solides interstitielles de chacun de ces systèmes. Par ailleurs, la pression en dihydrogène constitue un paramètre important car de nombreux hydrures ne se forment qu'à très haute pression. Afin d'améliorer la précision des modèles Calphad à pression élevée, le modèle de Lu et al. a été appliqué aux phases condensées des systèmes Ni-H, Rh-H et Mg-H. Ce modèle permet de déterminer la contribution à l'enthalpie libre du travail des forces de pression, en admettant en entrée aussi bien des calculs de phonon quasi-harmoniques que des données expérimentales. Enfin, la modélisation Calphad complète du système Ni-H a été réalisée en intégrant le modèle de Lu et al. L'enthalpie de formation calculée par DFT et l'enthalpie de mélange déterminée par CEM ont également été utilisées en données d'entrée, en complément des données expérimentales disponibles
Origine : Version validée par le jury (STAR)