Dopage aux métaux de transition de la nanoclay kaolinite : insights structuraux et électroniques par DFT

Résumé

En utilisant la théorie de la fonctionnelle de la densité avec correction de dispersion (DFT-D2), nous avons étudié des nanoclays kaolinite dopés aux métaux de transition Cr, Mn, Fe et Co. L'étude examine comment ces dopants modifient la géométrie cristalline et la structure électronique, évaluant les configurations ferromagnétique (FM), antiferromagnétique (AFM) et non magnétique (NM). Les calculs de volume cristallin, constantes de maille, longueurs de liaison, distribution de charge et moments de spin révèlent que Cr³⁺ et Fe³⁺ favorisent l'état AFM, Mn³⁺ stabilise les états AFM et FM, tandis que Co³⁺ préfère l'état NM. La substitution métallique élargit aussi le réseau et introduit des états localisés dans le gap de bande, ouvrant de nouvelles voies pour fonctionnaliser la kaolinite en catalyse et stockage d'énergie.

Contexte

La kaolinite, un minéral nanoclay largement distribué, possède une structure en feuillets et d'abondants groupes hydroxyles de surface qui confèrent des propriétés physiques uniques. Ces caractéristiques ont attiré des chercheurs en chimie des matériaux, sciences environnementales et physique minérale [1–11]. Sa polyvalence a mené à des applications dans les plastiques, la catalyse et la production de ciment, et des études récentes ont exploré la fonctionnalisation de la kaolinite comme support pour nanoparticules dans les matériaux à changement de phase pour l'énergie solaire [4,5] ou les revêtements d'oxydes conducteurs [9,12]. L'hybridation avec des nanoparticules fonctionnelles a amélioré l'activité photocatalytique de Pd-ZnO et la luminescence de CdS via effets synergétiques [6,7]. La modification de surface par ancrage de groupes fonctionnels [13,14] ou prétraitement par activation acide [2] améliore encore les performances de la kaolinite.

Des investigations expérimentales et théoriques ont élucidé la structure et l'énergétique des minéraux du groupe kaolin [15–22]. Les études d'adsorption de métaux lourds sur surfaces de kaolinite, incluant Cd, Cu, Hg et Ni(II) [23], ainsi que Pb(II) et espèces uranyles [24–28], ont cartographié les hiérarchies d'adsorption et profils énergétiques. Le dopage au Mg, Ca et Fe influence la pénétration d'H₂O et les interactions inter-feuillets [29]. Les analyses de structure électronique de kaolinite pure et défectueuse par DFT standard et fonctionnelles hybrides [30] ont révélé l'impact des défauts intrinsèques. Des modèles DFT récents ont aussi capturé l'évolution de la kaolinite lors de la déhydroxylation, déalumination et condensation de silice, soulignant comment le retrait d'Al modifie la géométrie et les propriétés électroniques [1,31,32].

Le dopage aux métaux de transition est une stratégie bien établie pour moduler les propriétés d'oxydes comme Al₂O₃ [33], TiO₂ [34], MOFs [35] et autres solides [36]. Étendre cette approche à la kaolinite offre le potentiel d'adapter ses caractéristiques structurales et électroniques pour des applications avancées. Dans ce travail, nous avons exploré systématiquement la substitution de Cr, Mn, Fe et Co aux sites Al³⁺ de la kaolinite, évaluant les résultats géométriques, électroniques et magnétiques par calculs DFT-D2.

Méthodes

Tous les calculs ont été effectués avec le code CASTEP [37] utilisant la fonctionnelle PBE (Perdew-Burke-Ernzerhof) GGA [38] augmentée des corrections de dispersion DFT-D2 de Grimme [39]. Un cutoff d'onde plane de 500 eV et un maillage Monkhorst-Pack 2 × 2 × 3 ont été utilisés pour les optimisations géométriques et évaluations de structure électronique. Les critères de convergence étaient fixés à 1,0 × 10⁻⁶ eV/atome pour la tolérance SCF, 0,03 eV/Å pour les forces atomiques, 0,05 GPa pour la contrainte et 0,001 Å pour le déplacement ionique. Les états de valence incluaient O(2s²2p⁴), Al(3s²3p¹), Cr(3s²3p⁶3d⁵4s¹), Mn(3d⁵4s²), Fe(3d⁶4s²) et Co(3d⁷4s²). Des pseudopotentiels ultrasoft (USPP) étaient employés pour Mn, Fe et Co, tandis que des pseudopotentiels norm-conservant étaient utilisés pour les autres éléments. Les optimisations géométriques utilisaient l'algorithme BFGS (Broyden-Fletcher-Goldfarb-Shanno), et les moments magnétiques initiaux étaient variés pour explorer les états FM, AFM et NM.

Résultats et Discussion

La structure de kaolinite de départ correspond à notre travail précédent [1]. La Figure 1 montre la supercellule 2 × 2 × 1 relaxée (quatre unités de kaolinite). Le réseau de kaolinite, Al₂Si₂O₅(OH)₄, consiste en une feuille octaédrique Al-O et une feuille tétraédrique Si-O liées via l'oxygène apical (Oₐ). Le tétraèdre Si-O comporte un Oₐ et trois O_b basaux, tandis que l'octaèdre Al-O incorpore deux atomes Oₐ et quatre oxygènes hydroxyles partagés avec les octaèdres voisins. Les groupes hydroxyles se divisent en OH_inter (surface) et OH_inner (internes). Par conséquent, le cristal contient des liaisons distinctes Si-Oₐ, Si-O_b, Al-O_inter, Al-O_inner et Al-O, comme illustré.

Vues de dessus (haut) et de côté (bas) de la kaolinite. Les liaisons Si-Oₐ (noir), Si-O_b (noir), M-OH_inter (rouge), M-OH_inner (vert) et M-O (bleu) sont indiquées par des pointillés.

Les forces de dispersion dominent la stabilisation inter-feuillets dans les argiles. Parmi les fonctionnelles hybrides, PBE-D2 reproduit précisément les paramètres de maille expérimentaux tout en restant efficiente [21,44,45]. La comparaison avec les données expérimentales (Tableau 1) montre que PBE-D2 réduit l'erreur relative du volume de maille à ~0,4 % (vs ~3,4 % pour PBE seul). Les vecteurs de maille a et b présentent une erreur ~0,4 %, et la distance inter-feuillets c se contracte de 0,17 Å (~2 %) sous corrections de dispersion. Les longueurs de liaison Si-Oₐ, Al-OH_inner et Al-O sont largement inchangées, tandis que Al-OH_inter (critique pour la chimie de surface) montre une amélioration substantielle. Les deux liaisons Al-O distinctes (~1,95 et 2,00 Å) proviennent du désaccord de maille entre feuilles Si-O et Al-O ; PBE et PBE-D2 surestiment légèrement ces distances (moyenne ~1,96 Å).

Distributions de liaisons dans la kaolinite dopée Cr, Mn, Fe, Co. Différents états magnétiques sont montrés pour chaque kaolinite dopée TM. Les longueurs moyennes de liaisons Si-O et M-O sont indiquées par pointillés ; les traits pleins soulignent les liaisons M-OH_inter (rouge), M-OH_inner (vert) et M-O (bleu) pour les configurations les plus stables.

La substitution aux métaux de transition a impliqué le remplacement d'Al³⁺ par TM³⁺ (Cr³⁺, Mn³⁺, Fe³⁺, Co³⁺). Les réponses structurales sont similaires entre PBE et PBE-D2, mais nous nous concentrons sur les résultats PBE-D2. Le Tableau 1 résume les paramètres de maille, longueurs de liaison, charges, spins et états magnétiques. Les différences d'énergie par atome TM entre états AFM et FM sont de 0,022 eV (Cr), –0,006 eV (Mn) et 0,094 eV (Fe). La kaolinite dopée au Co se stabilise uniquement en état NM.

Les volumes de maille s'élargissent dans l'ordre Mn- > Fe- > Cr- > pur > Co-, principalement en raison de liaisons M-O plus longues qu'Al-O, élargissant les vecteurs a et b. Les liaisons Si-O_b s'allongent aussi, distordant subtilement les angles α et β. La kaolinite-Mn montre une augmentation de volume de 1,4 % en FM versus AFM, alors que les volumes de kaolinite-Cr et -Fe sont largement insensibles à l'ordre magnétique. Les moments magnétiques des dopants reflètent ceux dans Al₂O₃ dopée TM [33], tandis que les charges de Mulliken sont légèrement plus élevées, suggérant une réactivité accrue.

La Figure 2 met en évidence l'évolution des longueurs de liaison lors du dopage. Les distances M-O et Si-O_b augmentent, et les liaisons Al-O scindées disparaissent après substitution par Cr et Co. L'ordre magnétique affecte significativement les distributions de liaisons Mn, mais les liaisons Cr et Fe sont relativement insensibles à l'état de spin.

Les densités d'états projetées (PDOS) pour Cr³⁺(d³), Mn³⁺(d⁴), Fe³⁺(d⁵) et Co³⁺(d⁶) (Fig. 3) et les densités de charge correspondantes (Fig. 4) illustrent la scission des états d en niveaux t₂g et e_g sous symétrie octaédrique. Les distorsions de Jahn-Teller sont évidentes pour Mn³⁺, où l'occupation de l'orbital d_z² repousse les ligands le long de l'axe z, allongeant les liaisons M-O. Cr³⁺, Fe³⁺ et Co³⁺ affichent des occupations t₂g/e_g équilibrées, produisant des effets Jahn-Teller minimes. Ces modifications électroniques introduisent des états 3d localisés dans le gap de bande, potentiellement améliorant les capacités catalytiques, de capture de CO, de chargement de médicaments et de stockage d'énergie de la kaolinite [50–57]. Des stratégies similaires pourraient s'étendre à la montmorillonite, perlite et talc.

PDOS des états les plus stables de kaolinite dopée TM. Le jaune surligne les orbitaux 3d occupés les plus hauts ; les flèches brunes et violettes pointent vers les premiers et seconds orbitaux 3d inoccupés les plus bas près du niveau de Fermi.

Densité de charge partielle des orbitaux 3d TM dans la kaolinite dopée TM, correspondant aux états PDOS. Isosurface à 0,02 e/ų.

Conclusions

Les calculs DFT-D2 révèlent que les dopants Cr³⁺ et Fe³⁺ stabilisent préférentiellement l'état AFM, Mn³⁺ favorise l'état FM, et Co³⁺ reste non magnétique. La substitution aux métaux de transition élargit le réseau et réorganise les distributions de liaisons M-O, tout en introduisant des états 3d qui se scindent plus fortement dans le gap de bande de la kaolinite. Ces insights fournissent une base solide pour adapter la structure et les propriétés électroniques de la kaolinite pour des applications catalytiques, d'adsorption et de stockage d'énergie.