Significance of the dispersion force for ferroelectric switching in ZnO and related materials

Resumen Técnico: Significado de la Fuerza de Dispersión para la Conmutación Ferroeléctrica en ZnO y Materiales Relacionados

Planteamiento del Problema
El óxido de zinc (ZnO) con estructura de wurtzita es un material polar de banda prohibida ancha, pero su barrera intrínseca de conmutación ferroeléctrica es prohibitivamente alta para aplicaciones prácticas. Aunque los materiales sustituidos como el Zn$_{0.5}$Mg$_{0.5}$O exhiben una polarización conmutable, los mecanismos a escala atómica que gobiernan esta conmutación de polarización concertada siguen siendo objeto de debate. Específicamente, la estabilidad de la fase intermedia hexagonal plana (h-ZnO o estructura C) es controvertida. Las observaciones experimentales de nanocristales de h-ZnO metaestables entran en conflicto con muchas predicciones computacionales que identifican a C como un estado de transición inestable en cristales infinitos. Una hipótesis central de este trabajo es que los enfoques computacionales estándar han tratado inadecuadamente las fuerzas de dispersión de van der Waals, que a menudo se omiten en sistemas de "átomos duros" como el ZnO y el MgO, lo que potencialmente conduce a discrepancias entre la teoría y el experimento respecto a la estabilidad de estas fases intermedias.

Metodología
Los autores emplearon un conjunto integral de métodos ab initio y de la teoría del funcional de la densidad (DFT) para mapear las superficies de energía potencial (PES) para la conmutación ferroeléctrica concertada en ZnO y Zn$_{0.5}$Mg$_{0.5}$O. La vía de conmutación se modeló como una transición entre dos estructuras de wurtzita (A y A') a través de un intermedio hexagonal plano (C) y posibles estados de transición (B).

El enfoque computacional incluyó:

1. Métodos de Función de Onda de Alto Nivel: Hartree-Fock (HF), teoría de la perturbación de segundo orden de Møller-Plesset (MP2) y la aproximación de fase aleatoria (RPA). Estos métodos fueron seleccionados para proporcionar puntos de referencia con tratamientos variables de la correlación electrónica y la dispersión, donde MP2 y RPA ofrecen mejoras sistemáticas sobre HF.
2. Teoría del Funcional de la Densidad (DFT): Una gama de funcionales que incluye PBESol, PBE0, HSE06 y r2SCAN.
3. Correcciones de Dispersión: Para abordar la omisión de la dispersión de largo alcance en los funcionales estándar, el estudio aplicó correcciones empíricas: D3(BJ) para HF, PBESol, PBE0 y HSE06; rVV10 para r2SCAN; y se probaron otras correcciones (D2, D3, ULG, optB88).
4. Análisis Estructural: Los cálculos se realizaron en celdas unitarias de 4 átomos (Zn$_2$O$_2$ y ZnMgO$_2$). La coordenada de reacción se definió mediante el desplazamiento fraccional de los cationes a lo largo del eje c. Se realizaron tanto optimizaciones estructurales completas como optimizaciones restringidas a vectores de red fijos (específicamente los vectores de red de la wurtzita aplicados a la estructura plana, denotada como C$_v$).
5. Convergencia: Se realizaron pruebas extensas de mallas de puntos k (hasta 12 $\times$ 12 $\times$ 8) y cortes de la onda plana (hasta 700 eV), particularmente para MP2 y RPA, para asegurar la precisión cuantitativa.

Resultos Clave

• Estabilidad de h-ZnO (ZnO): El consenso de los cálculos MP2 y RPA es que la estructura hexagonal plana C no es una fase metaestable en cristales de ZnO infinitos; es un estado de transición. Aunque las correcciones de dispersión en algunos funcionales de DFT (por ejemplo, PBE0/D3, HSE06/D3) estabilizan artificialmente a C lo suficiente como para que parezca metaestable, los métodos de alto nivel indican que sigue siendo inestable. Esto sugiere que las observaciones experimentales de nanocristales de h-ZnO estables se deben probablemente a efectos de tamaño finito o de solvente, más que a una estabilidad intrínseca del bulto (bulk).
• Estabilidad de h-Zn$_{0.5}$Mg$_{0.5}$O: En contraste, para el sistema sustituido con Mg, MP2 y RPA predicen consistentemente que la fase hexagonal plana C es, al menos, una fase metaestable y, en algunos casos, más estable que la fase de wurtzita A. Sin embargo, los métodos de DFT con correcciones de dispersión estándar (D3(BJ), rVV10) tienden a sobreestimar significativamente la estabilidad de C en comparación con MP2 y RPA. Solo el enfoque r2SCAN/rVV10 predice que A es ligeramente más estable que C, alineándose más cerca de los puntos de referencia de alto nivel.
• Papel de la Dispersión: Se encontró que las fuerzas de dispersión son más significativas para las estructuras hexagonales planas hipotéticas restringidas a los vectores de red de la fase de wurtzita (C$_v$). Para las estructuras optimizadas, las correcciones de dispersión generalmente estabilizan C con respecto a A, pero la magnitud de este efecto varía drásticamente entre diferentes funcionales de DFT.
• Desempeño de los Métodos:
  • MP2 y RPA proporcionan una descripción auto-consistente y equilibrada de las interacciones iónicas, covalentes y de dispersión, sirviendo como los puntos de referencia más fiables.
  • La DFT sin dispersión a menudo concuerda cualitativamente con MP2/RPA para ZnO, pero falla para Zn$_{0.5}$Mg$_{0.5}$O.
  • Las Correcciones de Dispersión Empírica en DFT (específicamente D3(BJ) y rVV10) a menudo sobre-estabilizan la fase plana, lo que genera inconsistencias cualitativas con la teoría de alto nivel. La corrección ULG y el funcional optB88 mostraron un mejor acuerdo con los resultados de MP2/RPA que D3(BJ).
• Parámetros de Red: Si bien las correcciones de dispersión son críticas para la energética, juegan un papel menor en la predicción de los parámetros de red de equilibrio, manteniéndose las desviaciones respecto al experimento pequeñas en la mayoría de los métodos.

Significancia y Reivindicaciones
El artículo afirma que un tratamiento preciso de las fuerzas de dispersión es esencial para describir la conmutación de la polarización y el comportamiento ferroeléctrico en materiales con estructura de wurtzita, incluso aquellos compuestos de "átomos duros" donde la dispersión se tradicionalmente omite. El estudio demuestra que:

1. La estabilidad percibida de h-ZnO está fuertemente influenciada por la elección de la corrección de dispersión, existiendo un consenso de la teoría de alto nivel (MP2/RPA) de que el h-ZnO es inestable en el bulto, resolviendo el conflicto con los datos de nanocristales experimentales al atribuir la estabilidad a factores ambientales.
2. Para Zn$_{0.5}$Mg$_{0.5}$O, la reducción de las barreras de conmutación y la estabilización de la fase intermedia son fenómenos reales, pero su descripción cuantitativa es altamente sensible al método computacional. La DFT estándar con correcciones de dispersión empíricas a menudo produce predicciones cualitativas poco fiables para este sistema.
3. MP2 y RPA ofrecen la precisión cuantitativa necesaria para distinguir entre la metaestabilidad real y los estados de transición en estos sistemas, resaltando que la suposición de "átomo duro" para omitir la dispersión es insuficiente para comprender los mecanismos de conmutación ferroeléctrica.

Los autores concluyen que, si bien la dispersión es un contribuyente significativo a la energética de la conmutación, no explica por sí sola la reducción de la barrera observada tras la sustitución con Mg; la interacción entre los cambios en el enlace iónico y covalente inducidos por la sustitución también es crítica y requiere métodos que traten la correlación electrónica de manera consistente.