Dopaje con Metales de Transición de Nanarcilla Caolinita: Perspectivas Estructurales y Electrónicas desde DFT

Resumen

Usando teoría de funcional de densidad corregida por dispersión (DFT-D2), investigamos nanarcillas de caolinita dopadas con metales de transición Cr, Mn, Fe y Co. El estudio examinó cómo estos dopantes modifican geometría cristalina y estructura electrónica, evaluando configuraciones ferromagnética (FM), antiferromagnética (AFM) y no magnética (NM). Cálculos de volumen cristalino, constantes de red, longitudes de enlace, distribución de carga y momentos de spin revelaron que Cr³⁺ y Fe³⁺ favorecen estado AFM, Mn³⁺ estabiliza ambos estados AFM y FM, mientras Co³⁺ prefiere estado NM. La sustitución metálica también expande la red e introduce estados localizados dentro del gap de banda, ofreciendo nuevas vías para funcionalizar caolinita en catálisis y almacenamiento de energía.

Antecedentes

La caolinita, un mineral de nanarcilla ampliamente distribuido, posee estructura en capas y abundantes grupos hidroxilo superficiales que confieren propiedades físicas únicas. Estas características han atraído investigadores en química de materiales, ciencia ambiental y física mineral. Su versatilidad ha llevado a aplicaciones en plásticos, catálisis y producción de cemento, y estudios recientes han explorado funcionalizar caolinita como soporte para nanopartículas en materiales de cambio de fase para energía solar o recubrimientos de óxidos conductores. Hibridación con nanopartículas funcionales ha mostrado mejorar actividad fotocatalítica de Pd-ZnO y luminiscencia de CdS vía efectos sinérgicos. Modificación superficial mediante anclaje de grupos funcionales o pretratamiento de activación ácida mejora aún más rendimiento de caolinita.

Investigaciones experimentales y teóricas han elucidado estructura y energética de minerales grupo caolin. Estudios de adsorción de metales pesados en superficies de caolinita, incluyendo Cd, Cu, Hg y Ni(II), así como Pb(II) y especies uranilo, han mapeado jerarquías de adsorción y perfiles energéticos. Dopaje con Mg, Ca y Fe ha mostrado influir en penetración de H₂O e interacciones intercapas. Análisis de estructura electrónica de caolinita prístina y defectuosa usando DFT estándar y funcionales híbridos revelaron impacto de defectos intrínsecos. Modelos DFT recientes también capturaron evolución de caolinita durante deshidroxilación, dealuminización y condensación de sílice, destacando cómo remoción de Al altera geometría y propiedades electrónicas.

Dopaje con metales de transición es estrategia establecida para ajustar propiedades de óxidos como Al₂O₃, TiO₂, MOFs y otros sólidos. Extender este enfoque a caolinita ofrece potencial para adaptar sus características estructurales y electrónicas para aplicaciones avanzadas. En este trabajo, exploramos sistemáticamente sustitución de Cr, Mn, Fe y Co en sitios Al³⁺ de caolinita, evaluando resultados geométricos, electrónicos y magnéticos usando cálculos DFT-D2.

Métodos

Todos los cálculos se realizaron con código CASTEP empleando funcional Perdew-Burke-Ernzerhof (PBE) GGA aumentado por correcciones de dispersión DFT-D2 de Grimme. Corte de onda plana de 500 eV y malla Monkhorst-Pack 2 × 2 × 3 k-puntos se usaron para optimizaciones geométricas y evaluaciones de estructura electrónica. Criterios de convergencia: 1,0 × 10⁻⁶ eV/átomo tolerancia SCF, 0,03 eV/Å para fuerzas atómicas, 0,05 GPa para esfuerzo y 0,001 Å para desplazamiento iónico. Estados de valencia incluyeron O(2s²2p⁴), Al(3s²3p¹), Cr(3s²3p⁶3d⁵4s¹), Mn(3d⁵4s²), Fe(3d⁶4s²), Co(3d⁷4s²). Se emplearon pseudopotenciales ultrasuaves (USPP) para Mn, Fe y Co, mientras pseudopotenciales norm-conservativos para elementos restantes. Optimizaciones geométricas usaron algoritmo Broyden-Fletcher-Goldfarb-Shanno (BFGS), y momentos magnéticos iniciales variaron para explorar estados FM, AFM y NM.

Resultados y Discusión

La estructura inicial de caolinita corresponde a nuestro trabajo previo. Figura 1 muestra supercelda 2 × 2 × 1 relajada (cuatro unidades caolinita). Marco de caolinita, Al₂Si₂O₅(OH)₄, consiste en lámina octaédrica Al-O y lámina tetraédrica Si-O unidas vía oxígeno apical (Oₐ). Tetraedro Si-O presenta un Oₐ y tres O_b basales, mientras octaedro Al-O incorpora dos átomos Oₐ y cuatro oxígenos hidroxilo compartidos con octaedros vecinos. Grupos hidroxilo dividen en OH_inter (superficial) y OH_inner (interno). Consecuentemente, cristal contiene enlaces distintos Si-Oₐ, Si-O_b, Al-O_inter, Al-O_inner y Al-O, como ilustra figura.

Fuerzas de dispersión dominan estabilización intercapas en arcillas. Entre funcionales híbridos, PBE-D2 reproduce con precisión parámetros de red experimentales manteniéndose computacionalmente eficiente. Comparación con datos experimentales (Tabla 1) muestra PBE-D2 reduce error relativo en volumen celda unidad a ~0,4 % (vs ~3,4 % para PBE simple). Vectores de red a y b exhiben ~0,4 % error, y distancia intercapas c se contrae 0,17 Å (~2 %) bajo correcciones de dispersión. Longitudes de enlace Si-Oₐ, Al-OH_inner y Al-O permanecen largely unchanged, mientras Al-OH_inter (crítico para química superficial) muestra mejora sustancial. Dos enlaces Al-O distintos (~1,95 y 2,00 Å) surgen de desajuste de red entre láminas Si-O y Al-O; tanto PBE como PBE-D2 ligeramente sobreestiman estas distancias (promedio ~1,96 Å).

Sustitución de metal de transición involucró reemplazar Al³⁺ con TM³⁺ (Cr³⁺, Mn³⁺, Fe³⁺, Co³⁺). Respuestas estructurales son similares entre PBE y PBE-D2, pero nos enfocamos en resultados PBE-D2. Tabla 1 resume parámetros de red, longitudes de enlace, cargas, spins y estados magnéticos. Diferencias de energía por átomo TM entre estados AFM y FM son 0,022 eV (Cr), –0,006 eV (Mn) y 0,094 eV (Fe). Caolinita dopada con Co estabiliza solo en estado NM.

Volúmenes de celda unidad se expanden en orden Mn- > Fe- > Cr- > prístina > Co-, principalmente debido a enlaces M-O más largos respecto a Al-O, agrandando vectores a y b. Enlaces Si-O_b también se alargan, distorsionando sutilmente ángulos α y β. Caolinita-Mn muestra aumento 1,4 % volumen en FM versus AFM, mientras volúmenes Caolinita-Cr y -Fe son largely unaffected por ordenamiento magnético. Momentos magnéticos de dopantes reflejan los en Al₂O₃ dopado con TM, mientras cargas Mulliken son ligeramente mayores, sugiriendo reactividad增强.

Figura 2 destaca evolución de longitudes de enlace tras dopaje. Distancias M-O y Si-O_b aumentan, y enlaces Al-O divididos desaparecen tras sustitución Cr y Co. Ordenamiento magnético afecta significativamente distribuciones de enlace Mn-, pero enlaces Cr y Fe son relativamente insensibles a estado de spin.

Densidad de estados proyectada (PDOS) para Cr³⁺(d³), Mn³⁺(d⁴), Fe³⁺(d⁵) y Co³⁺(d⁶) (Fig. 3) y densidades de carga correspondientes (Fig. 4) ilustran división d-estados en niveles t₂g y e_g bajo simetría octaédrica. Distorsiones Jahn-Teller evidentes para Mn³⁺, donde ocupación orbital d_z² repele ligandos a lo largo eje z, alargando enlaces M-O. Cr³⁺, Fe³⁺ y Co³⁺ muestran ocupaciones balanceadas t₂g/e_g, produciendo efectos Jahn-Teller mínimos. Estas modificaciones electrónicas introducen estados 3d localizados dentro del gap de banda, potencialmente mejorando capacidades catalíticas, captura CO, carga fármacos y almacenamiento energético de caolinita. Estrategias similares podrían extenderse a montmorillonita, perlita y talco.

Conclusiones

Cálculos DFT-D2 revelan que dopantes Cr³⁺ y Fe³⁺ estabilizan preferentemente estado AFM, Mn³⁺ favorece estado FM y Co³⁺ permanece no magnético. Sustitución de metales de transición expande la red y reorganiza distribuciones de enlaces M-O, mientras introduce estados 3d que se dividen más fuertemente dentro del gap de banda de caolinita. Estos insights proporcionan base sólida para adaptar estructura y propiedades electrónicas de caolinita para aplicaciones catalíticas, adsortivas y de almacenamiento de energía.