First principles study of chalcogen vacancy effect on the optoelectronic and photocatalytic properties of transition metal dichalcogenides monolayers

Resumen Técnico: Estudio de Primeros Principios sobre los Efectos de las Vacancias de Calcógeno en Monocapas de TMDC

Planteamiento del Problema
Los dicalcogenuros de metales de transición bidimensionales con la fórmula MX₂ (donde M = Mo, W y X = S, Se) son candidatos prometedores para dispositivos optoelectrónicos y aplicaciones fotocatalíticas, tales como la división del agua y la reducción de CO₂. Sin embargo, su eficiencia práctica se ve a menudo obstaculizada por altas tasas de recombinación de portadores, baja capacidad de captación de luz y una cinética de reacción lenta. Además, durante la síntesis mediante métodos como la deposición química de vapor (CVD) o la exfoliación mecánica, las vacancias de calcógeno son inevitables. Si bien los planos basales prístinos en TMDCs estables suelen ser inactivos, el impacto específico de estos defectos inevitables en las propiedades estructurales, optoelectrónicas y fotocatalíticas de las monocapas requiere una investigación sistemática para determinar si degradan el rendimiento o si ofrecen ventajas funcionales.

Metodología
Los autores emplearon cálculos de mecánica cuántica de primeros principios para investigar prototípicas monocapas de MX₂ (MoS₂, MoSe₂, WS₂ y WSe₂).

• Modelado de Defectos: Utilizando el Alloy Theoretical Automated Toolkit (ATAT) junto con el Vienna ab initio Simulation Package (VASP), los autores utilizaron aproximaciones de expansión de clúster para identificar las concentraciones de vacancia energéticamente más favorables mediante la construcción de envolventes convexas de energía (energy convex hulls). Las concentraciones identificadas fueron 6.25% para MoS₂, MoSe₂ y WSe₂, y 1.56% para WS₂.
• Detalles Computacionales: Los cálculos se realizaron mediante la Teoría del Funcional de la Densidad (DFT) con la aproximación de gradiente generalizado Perdew-Burke-Ernzerhof (PBE) y el funcional híbrido Heyd-Scuseria-Ernzerhof (HSE06) para predicciones más precisas de la brecha de banda (band gap). Las interacciones de Van der Waals se corrigieron utilizando el método DFT-D2 de Grimme.
• Análisis de Estabilidad: La estabilidad estructural se evaluó mediante cálculos de dispersión de fonones (usando PHONOPY) y simulaciones de dinámica molecular ab initio (AIMD) a 300 K utilizando un termostato de Nosé-Hoover sobre una trayectoria de 8 ps.
• Evaluación de Propiedades: El estudio analizó estructuras de bandas electrónicas, la asimetría del potencial electrostático (promedio planar a lo largo del eje c), funciones de trabajo, coeficientes de absorción óptica y la alineación de los bordes de banda con respecto a los potenciales redox del agua y los potenciales de reducción de CO₂.

Resultos Clave

• Estabilidad Estructural: Las optimizaciones estructurales revelaron que las longitudes de enlace M–X alrededor de las vacancias difieren solo ligeramente de las estructuras prístinas, lo que indica una deformación local mínima. La dispersión de fonones y las simulaciones AIMD confirmaron que las estructuras defectuosas son tanto dinámica como térmicamente estables a temperaturas de operación, sin que se observara ruptura de enlaces.
• Propiedades Electrónicas: La introducción de vacancias de calcógeno crea estados electrónicos dentro de la brecha de banda, lo que conduce a una reducción en el ancho de la brecha de banda en comparación con las monocapas prístinas. Los cálculos HSE06 indicaron que el Máximo de la Banda de Valencia (VBM) y el Mínimo de la Banda de Conducción (CBM) derivan principalmente de los orbitales d del metal de transición alrededor de la vacancia, con contribuciones menores de los átomos de calcógeno.
• Asimetría Electrostática y Separación de Carga: Un hallazgo crítico es la inducción de asimetría en el potencial electrostático. La presencia de vacancias en una capa de calcógeno crea una diferencia en la función de trabajo ($\delta\phi$) entre el lado de la vacancia y el lado prístino. Esta asimetría genera un campo eléctrico interno que facilita la separación de carga y reduce la recombinación, una característica ausente en las monocapas simétricas y prístinas. El efecto fue más pronunciado en las estructuras basadas en Se debido a la mayor polarizabilidad del Selenio.
• Propiedades Ópticas: La introducción de vacancias causó un desplazamiento al rojo (red shift) en el primer pico de absorción para las monocapas de MoX₂, permitiendo la absorción de fotones de menor energía, incluyendo la radiación infrarroja, en comparación con sus contrapartes prístinas.
• Desempeño Fotocatalítico:
  • División del Agua: Todas las monocapas defectuosas poseen brechas de banda superiores a 1.23 eV y posiciones de CBM adecuadas para la Reacción de Evolución de Hidrógeno (HER). Sin embargo, el Máximo de la Banda de Valencia (VBM) se sitúa por encima del potencial requerido para la Reacción de Evolución de Oxígeno (OER), lo que sugiere que, si bien la HER es factible, se requiere un cocatalizador para asistir la OER.
  • Reducción de CO₂: El estudio evaluó la capacidad de reducir CO₂ en diversos productos (HCOOH, CO, HCHO, CH₃OH, CH₄). Se identificó que las monocapas defectuosas de WS₂ y WSe₂ son capaces de reducir CO₂ en todos los productos químicos útiles enumerados. MoSe₂ fue capaz de producir HCHO, CH₃OH y CH₄, mientras que MoS₂ se limitó a CH₃OH y CH₄.

Significancia y Afirmaciones
El artículo argumenta que las monocapas de TMDC con vacancias no deben considerarse productos defectuosos para ser descartados. En su lugar, los resultados demuestran que estos defectos pueden mejorar las capacidades fotocatalíticas en comparación con las contrapartes prístinas al reducir las brechas de banda, inducir una asimetría electrostática beneficiosa para la separación de carga y permitir la absorción de luz visible. Los autores afirman que sus hallazgos proporcionan directrices para la síntesis experimental de monocapas de MX₂ con ingeniería de vacancias específicamente diseñadas para dispositivos optoelectrónicos y aplicaciones fotocatalíticas, destacando particularmente a WS₂ y WSe₂ como materiales eficientes para la reducción de CO₂ y la división del agua (con asistencia de cocatalizador para la OER).