Phase stability and ionic transport in post-spinel CaV$_2$O$_4$ cathode

Resumen Técnico: Estabilidad de Fase y Transporte Iónico en el Cátodo de Post-Espinela CaV₂O₄

Planteamiento del Problema
Las baterías de calcio (CB, por sus siglas en inglés) ofrecen una alternativa sostenible y de alta densidad energética a la tecnología de ion-litio, aprovechando la abundancia terrestre y la alta capacidad teórica del calcio. Sin embargo, el despliegue práctico de las CB se ve actualmente obstaculizado por la falta de cátodos de intercalación de alto rendimiento. Aunque el cribado computacional identificó al post-espinela CaV₂O₄ como un candidato prometedor, la validación experimental reveló limitaciones significativas: el material exhibe solo una extracción parcial de calcio (alcanzando ~65 mAh g⁻¹ a 298 K y ~155 mAh g⁻¹ a 323 K) y no logra la desintercalación completa incluso a temperaturas elevadas. Las causas subyacentes de esta capacidad electroquímica restringida siguen sin estar claras. Específicamente, es una cuestión abierta si la incapacidad de alcanzar la capacidad teórica total se debe a una inestabilidad termodinámica (imposibilidad de ciertas fases) o a un atrapamiento cinético (altas barreras de migración y separación de fases). Además, la reestructuración dinámica de la red huésped durante los ciclos, indicada por distintos escalones de voltaje y cambios de pendiente, requiere una explicación sistemática tanto termodinámica como cinética.

Metodología
Este estudio emplea un marco computacional integral que integra cálculos de primeros principios con mecánica estadística para investigar la termodinámica y el transporte iónico de CaxV₂O₄ (0 ≤ x ≤ 1).

• Cálculos de Primeros Principios: Se realizaron relajaciones de geometría con espín polarizado utilizando la teoría del funcional de la densidad (DFT) con el funcional de intercambio-correlación SCAN+U (U = 1.0 eV para los electrones V 3d) para calcular las energías de formación y las barreras de migración.
• Expansión de Clústeres (CE): Se construyó un Hamiltoniano de CE ajustándolo a un conjunto de entrenamiento de 262 energías de formación calculadas mediante DFT a través de diversas configuraciones de vacantes de Ca. El modelo utilizó superceldas (1×1×1, 1×1×2, 1×1×3) para capturar interacciones de largo alcance.
• Simulaciones de Monte Carlo: Se realizaron simulaciones de Monte Carlo en el Gran Canónico (GCMC) utilizando el Hamiltoniano de CE para generar el diagrama de fase temperatura-composición y los perfiles de voltaje en equilibrio. Las simulaciones de Monte Carlo Canónico (CMC) con recocido simulado se utilizaron para verificar las configuraciones de estado fundamental.
• Barreras de Migración: Se realizaron cálculos de Nudged Elastic Band (NEB) para determinar las barreras de migración de Ca²⁺ en regiones de una sola fase termodinámicamente estables. Para las composiciones donde el DFT-NEB era computacionalmente prohibitivo, potenciales interatómicos aprendidos por máquina (MLIP) y modelos de aprendizaje transferido (TL) proporcionaron estimaciones cualitativas.

Resultos Clave

1. Estabilidad de Fase y Diagrama: El estudio construye un diagrama de fase detallado temperatura-composición que revela seis regiones de una sola fase distintas (α a ζ) correspondientes a contenidos específicos de Ca (x = 0, 0.25, 0.5, 0.67, 0.83 y 1.0).

   • El casco convexo (convex hull) a 0 K identifica estados fundamentales estables en x = 0.25, 0.5 y 0.67.
   • La fase ε (x ∼ 0.83) es única; se estabiliza solo a temperaturas elevadas (370–590 K) mediante reacciones invariantes eutécticas y peritécticas, careciendo de un estado fundamental a 0 K.
   • Existen regiones de dos fases significativas, particularmente la región γ + δ (0.5 ≤ x ≤ 0.67), que persiste hasta temperaturas superiores a 1000 K, lo que indica que es inevitable bajo condiciones operativas electroquímicas.
   • La solubilidad de Ca en V₂O₄ (fase α) y de las vacantes en CaV₂O₄ (fase ζ) es insignificante a 298 K, pero aumenta a temperaturas más altas.

2. Perfiles de Voltaje: Los perfiles de voltaje calculados exhiben escalones distintos correspondientes a las transiciones de fase.

   • Se observan cuatro escalones bien definidos a 0 K y 300 K, correspondientes a las regiones de coexistencia de dos fases: α+β (4.60 V), β+γ (4.59 V), γ+δ (1.40 V) y δ+ζ (0.60 V).
   • Un sutil escalón de voltaje aparece cerca de x ∼ 0.83 a temperaturas superiores a 370 K, correspondiente a la emergencia de la fase ε. Esto se alinea con las observaciones experimentales de cambios sutiles de pendiente en el rango de x ∼ 0.80–0.85.
   • El voltaje promedio calculado (2.73 V) es inferior a los promedios experimentales, una discrepancia atribuida a reacciones secundarias y sobrepotenciales en los montajes experimentales.

3. Transporte Iónico y Barreras Cinéticas: Los cálculos de las barreras de migración revelan una dependencia fuerte y no monotónica de la composición.

   • Barreras Altas: La fase γ (x = 0.5) exhibe una barrera de migración máxima de ~1160 meV, y la fase α (x = 0) muestra una barrera de ~1050 meV. Estos valores altos se atribuyen al fuerte ordenamiento de Ca-vacante (arreglo tipo cúbico de diamante) impulsado por la atracción electrostática, lo que crea un severo cuello de botella cinético.
   • Barreras Bajas: Las fases δ (x = 0.67) y ζ (x = 1.0) exhiben barreras más bajas (~661 meV y ~600 meV, respectivamente), lo que sugiere una movilidad favorable en estas regiones.
   • La fase ε (x = 0.83) se predice con barreras intermedias similares a la fase ζ, lo que implica que no representa un cuello de botella cinético.

Significancia y Reivindicaciones
El artículo establece una base teórica para la baja capacidad electroquímica observada experimentalmente en CaV₂O₄. Los autores concluyen que la capacidad limitada no es una imposibilidad termodinámica, sino el resultado de limitaciones cinéticas derivadas de la interacción entre la separación de fases persistente y el transporte iónico dependiente de la composición.

Específicamente, el estudio identifica que durante la extracción de Ca, el sistema entra en la región de dos fases γ + δ. Mientras que la fase δ permite el transporte de iones, la fase γ (x = 0.5) coexistente actúa como un trampa cinética debido a su extremadamente alta barrera de migración (~1.16 eV) causada por el fuerte ordenamiento de Ca-vacante. En consecuencia, la capacidad electroquímica accesible a 298 K está limitada cinéticamente a aproximadamente la mitad de la capacidad teórica (correspondiente a una extracción hasta x ≈ 0.5), lo cual concuerda con los datos experimentales que muestran ~0.3 mol de extracción de Ca.

Los autores sugieren que estrategias como el dopaje de cationes (para romper el ordenamiento de Ca-vacante y suavizar el paisaje de energía potencial) y la reducción del tamaño de partícula (para acceder a soluciones sólidas fuera del equilibrio y sortear las barreras de nucleación) podrían mitigar estos cuellos de botella cinéticos. Además, operar a temperaturas elevadas puede mejorar la capacidad al aumentar la solubilidad y reducir la polarización, de acuerdo con las tendencias experimentales. Este trabajo proporciona un marco cuantitativo para comprender las limitaciones de capacidad en huéspedes de intercalación basados en Ca y ofrece estrategias de diseño específicas para permitir el uso práctico de CaV₂O₄ como cátodo.