Density functional study of native point defects in CaO

Resumen Técnico: Estudio de Funcional de la Densidad sobre Defectos de Punto Nativos en CaO

Planteamiento del Problema
El óxido de calcio (CaO) es un material candidato para óxidos de puerta en estructuras de metal–óxido–semiconductor debido a su amplio ancho de banda (~7.1 eV) y su alto constante dieléctrica. Sin embargo, su fiabilidad eléctrica está fuertemente influenciada por defectos de punto nativos, particularmente las vacantes de oxígeno ($V_O$). Por el contrario, se ha propuesto que los defectos en el CaO son candidatos prometedores para defectos cuánticos debido a sus largos tiempos de coherencia de espín, sugiriéndose complejos específicos como $X_{Ca}–V_O$ como análogos al centro NV del diamante. A pesar de estudios previos de la Teoría del Funcional de la Densidad (DFT), aspectos clave de los defectos nativos en CaO siguen sin estar claros, incluyendo la estabilidad relativa de defectos aislados frente a complejos de defectos, las barreras de migración que gobiernan la difusión y el recocido, y los orígenes específicos de los picos de absorción y emisión óptica observados experimentalmente. Además, existen discrepancias en la literatura con respecto a los estados de carga estables de las vacantes de calcio ($V_{Ca}$), donde algunos estudios de funcionales híbridos predicen estados neutros o positivamente cargados estables, lo cual es contrario al comportamiento típico de los óxidos.

Metodología
Los autores realizaron cálculos sistemáticos de primer principio mediante DFT utilizando el paquete de simulación Vienna Ab initio (VASP).

• Funcionales: Los cálculos principales utilizaron el funcional meta-GGA fuertemente restringido y apropiadamente normalizado (SCAN) para determinar las relajaciones estructurales, las energías de formación y las barreras de migración. Para abordar las energías de transición óptica y corregir la subestimación del ancho de banda inherente a los funcionales semilocales, se realizaron cálculos de un solo disparo (one-shot) con el funcional híbrido (HSE) sobre las geometrías relajadas con SCAN. El parámetro de mezcla ($\alpha$) de HSE se ajustó a 0.504 para reproducir el ancho de banda experimental de 7.09 eV.
• Modelos: Las propiedades de la red (bulk) se calcularon utilizando celdas unitarias primitivas. Los cálculos de defectos emplearon superceldas de 216 átomos (expansión 3×3×3) con muestreo de punto $\Gamma$ para HSE y mallas de Monkhorst-Pack de 2×2×2 para SCAN.
• Termodinámica: Las energías de formación se calcularon a través del ancho de banda para diversos niveles de Fermi bajo condiciones de potencial químico pobre en O, rico en O e intermedio. Se determinaron los niveles de transición de estado de carga.
• Dinámica: Las barreras de migración se calcularon utilizando el método de nudged elastic band de imagen ascendente (cNEB) para estimar las temperaturas de recocido donde los defectos se vuelven móviles.
• Óptica: Las energías de transición óptica se evaluaron dentro del marco de Franck–Condon, combinando las geometrías de SCAN con las estructuras electrónicas de HSE.

Contribuciones Clave y Resultados

1. Estabilidad y Formación de Defectos:

   • Vacantes de Oxígeno ($V_O$): Las $V_O$ son favorecidas bajo condiciones pobres en O. Actúan como donadores profundos con niveles de transición de estado de carga en 0.95 eV (+/0) y 2.65 eV (2+/+) por debajo del mínimo de la banda de conducción (CBM). El estado neutro ($V_O^0$) exhibe un estado fundamental de espín-1/2.
   • Vacantes de Calcio ($V_{Ca}$): Contrario a algunos estudios recientes de HSE que predicen una $V_{Ca}$ neutra o positiva estable, los cálculos de SCAN (y HSE con ajuste estándar) indican que la $V_{Ca}$ es termodinámicamente estable solo en el estado de carga 2- a través de todo el ancho de banda. No se encuentran estados de defecto dentro de la banda para $V_{Ca}^{2-}$. Los autores atribuyen la discrepancia en otros estudios a la alta sensibilidad de la energética de los defectos al parámetro exacto de mezcla de intercambio y al alineamiento de los bordes de banda.
   • Complejos: Bajo condiciones de crecimiento realistas (potenciales químicos intermedios), tanto las $V_O$ como las $V_{Ca}$ tienen bajas energías de formación. El estudio confirma la estabilidad termodinámica de los complejos de vacantes, específicamente $(V_{Ca}–V_O)^0$ y $(V_{Ca}–V_O)^-$. Las energías de unión son positivas (1.95 eV y 1.07 eV, respectivamente), lo que indica una formación de complejos favorable. El complejo $(V_{Ca}–V_O)^-$ es paramagnético.

2. Migración y Recocido:

   • Las barreras de migración para vacantes aisladas son altas: 1.88 eV para $V_O^{2+}$ y 2.09 eV para $V_{Ca}^{2-}$.
   • Las temperaturas de recocido ($T_a$) calculadas para estas vacantes aisladas oscilan entre ~677 K y 1138 K. Esto sugiere que las vacantes aisladas son inmóviles a temperatura ambiente pero pueden difundirse durante el recocido a alta temperatura.
   • Los intersticiales ($O_i$, $Ca_i$) y los defectos tipo antisitio fueron encontrados con barreras de migración bajas (por ejemplo, 0.13 eV para $O_i^0$) o altas energías de formación, lo que los hace improbables de persistir como defectos aislados o formar antisitios estables; en su lugar, tienden a formar complejos de intersticial-vacante.

3. Propiedades Ópticas:

   • Al combinar las geometrías de SCAN con las energías de transición de HSE, los autores asignan los picos ópticos experimentales:
     • Absorción a 340 nm (3.65 eV): Atribuida a la transición $(V_{Ca}–V_O)^- \to (V_{Ca}–V_O)^0 + e^-$.
     • Absorción a 260 nm (4.77 eV): Atribuida al complejo $(V_{Ca}–V_O)^0$.
     • Absorción a 400 nm (3.11 eV): Aunque está cerca de la energía de ionización de la $V_O^0$ aislada (2.69 eV), los autores argumentan que este pico se interpreta mejor como proveniente de un estado tipo donador hidrogenoide ($V_O^+$ vinculando débilmente un electrón) en lugar de una $V_O^0$ aislada termodinámicamente estable, lo cual es consistente con la dependencia de la temperatura de la señal del centro F.
     • Emisión a 605 nm (2.05 eV): Atribuida a la transición $(V_{Ca}–V_O)^0 + e^- \to (V_{Ca}–V_O)^-$.
     • Emisión a 370 nm (3.35 eV): Consistente con transiciones que involucran a $V_O^+$ y $(V_{Ca}–V_O)^-$.

Significancia
El artículo proporciona una comprensión microscópica consistente de los defectos nativos en CaO, resolviendo las ambigüedades respecto a la estabilidad de los estados de carga de la vacante de calcio y la naturaleza de los complejos de defectos. Establece que los complejos de vacantes son termodinámicamente estables y probablemente se formen bajo condiciones experimentales, pudiendo actuar como fuentes de ruido magnético en aplicaciones cuánticas debido a su naturaleza paramagnética. Además, el estudio vincula con éxito especies de defectos y complejos específicos con los espectros ópticos observados experimentalmente, ofreciendo una base para la futura ingeniería de defectos y el desarrollo de defectos cuánticos en CaO. El trabajo también destaca la importancia de la consistencia metodológica (por ejemplo, usar SCAN para la energética y HSE para las transiciones) sobre el simple ajuste del ancho de banda al predecir las propiedades de los defectos.