Rate Equation for the Transfer of Interstitials across Interfaces between Equilibrated Crystals

Resumen Técnico: Ecuación de Tasa para la Transferencia de Intersticiales entre Interfaces entre Cristales en Equilibrio

Planteamiento del Problema
El artículo aborda una discrepancia entre las observaciones experimentales y los modelos cinéticos existentes respecto a la carga y descarga de soluciones sólidas intersticiales, específicamente hidruros metálicos (p. ej., Pd, PdPt). Los experimentos indican que el tiempo característico ($T$) para los cambios de composición aumenta drásticamente (en órdenes de magnitud) cuando el rango de composición incluye una transformación de fase o la vecindad de un punto crítico de brecha de miscibilidad. La ecuación estándar de Butler-Volmer (BV), ampliamente utilizada para analizar la cinética de hidruros controlada electroquímicamente, no logra explicar esta observación. El marco de BV asume la transferencia entre sitios con ocupación y energía de estado fundamental fijas a priori, tratando la fuerza impulsora únicamente como la diferencia en los potenciales químicos ($\Delta \mu_S$). Sin embargo, en los hidruros metálicos a nanoescala, la difusión interna es órdenes de magnitud más rápida que la transferencia interfacial, lo que significa que las fases voluminosas alcanzan el equilibrio interno de forma casi instantánea. En consecuencia, la transferencia interfacial ocurre entre conjuntos ergonicos extendidos de partículas de soluto donde las ocupaciones y las energías de estado fundamental fluctúan, un escenario que no es capturado por los modelos estándar de probabilidad de salto.

Metodología
El autor deriva una nueva ley de tasa cinética basada en la mecánica estadística y el principio de balance detallado, en lugar de probabilidades de salto fenomenológicas. El sistema se modela como un compuesto de dos fases de solución sólida intersticial (A y B) separadas por una barrera de energía con un estado de transición ($T$).

1. Marco Estadístico: Las fases se tratan como sistemas canónicos a temperatura constante ($T$) con números conservados de partículas de soluto ($N_S$) y sitios de la matriz ($N_M$). Se asume que la difusión volumétrica es lo suficientemente rápida como para mantener la ergodicidad interna dentro de cada fase.
2. Análisis del Estado de Transición: La probabilidad de una "configuración de transferencia" ($C_p$) —donde una partícula se encuentra en un estado activado en el sitio de origen y el sitio de destino está vacío— se calcula sumando todos los microestados del sistema restringido.
3. Descomposición de la Energía Libre: El cambio de energía libre ($\delta F$) requerido para alcanzar la configuración de transferencia se computa mediante un proceso virtual de cuatro pasos:
   • Comprimir el conjunto de soluto en la fase de origen para liberar un sitio (costo entrópico).
   • Comprimir el conjunto de soluto en la fase de destino para liberar un sitio (costo entrópico).
   • Eliminar una partícula de la fase de origen hacia un estado de referencia (trabajo de potencial químico).
   • Insertar la partícula en el estado de transición específico (costo energético).
4. Derivación del Flujo: La densidad de flujo neto ($J_B$) se deriva como la diferencia entre los flujos hacia adelante y hacia atrás, utilizando las probabilidades calculadas de las configuraciones de transferencia.

Contribuciones Clave y Resultados
La contribución principal es la derivación de una nueva ecuación de tasa (Ec. 15) para el flujo interfacial de partículas de carga neutra:
$$J_B = \nu_0 \rho_T (1 - x_A)(1 - x_B) \exp\left(-\frac{\epsilon_T}{k_B T}\right) \left( \exp\left(\frac{\mu_S^A}{k_B T}\right) - \exp\left(\frac{\mu_S^B}{k_B T}\right) \right)$$

Características clave de este resultado incluyen:

• Dependencia Explícita de los Potenciales Químicos: A diferencia de la ecuación BV, que depende de la diferencia en los potenciales químicos ($\Delta \mu_S$), esta ley de tasa depende explícitamente de los potenciales químicos individuales ($\mu_S^A$ y $\mu_S^B$) de las dos fases.
• Dependencia de la Fracción de Vacantes: El flujo está modulado explícitamente por las fracciones de vacantes $(1-x_A)$ y $(1-x_B)$, reflejando la probabilidad estadística de encontrar un sitio vacío en un conjunto ergódico.
• Balance Detallado: La ecuación satisface naturalmente el principio de balance detallado, desvaneciéndose cuando $\mu_S^A = \mu_S^B$.
• Explicación de la Ralentización Cinética: El artículo demuestra que cerca de un punto crítico de brecha de miscibilidad o espinodal, la susceptibilidad del soluto ($\chi_S = dx/d\mu_S$) diverge. Dado que el tiempo característico $T$ es inversamente proporcional al producto de la densidad de flujo de intercambio y la susceptibilidad ($T^{-1} \propto J_{eq}^\mu / \chi_S$), la divergencia de $\chi_S$ conduce a una divergencia en $T$. Esto proporciona una explicación teórica para la ralentización drástica observada experimentalmente en los tiempos de carga cerca de puntos críticos.
• Asimetría: La ley derivada predice una asimetría en el tiempo característico $T$ a ambos lados del punto crítico, una característica observada experimentalmente pero que no es reproducida por la simétrica ecuación BV.

Significancia
El artículo afirma que su ley de tasa derivada ofrece una descripción físicamente más consistente de la transferencia interfacial en sistemas donde las fases voluminosas se equilibran rápidamente en comparación con la interfaz. Al tratar los sitios de transferencia como acoplados a sistemas canónicos ergódicos en lugar de sitios fijos, el modelo da cuenta de las fluctuaciones de las energías de estado fundamental y las ocupaciones inherentes a los soluciones sólidas con interacciones. El trabajo sugiere que el fallo de la ecuación de Butler-Volmer para explicar las anomalías cinéticas en los hidruros metálicos proviene de su negligencia de la dependencia explícita de las composiciones de equilibrio y la naturaleza estadística de la disponibilidad de vacantes. Los resultados proporcionan una base teórica para comprender la cinética de carga/descarga en materiales de almacenamiento de energía, particularmente cerca de transformaciones de fase, sin requerir ajustes ad hoc del coeficiente de transferencia o la energía de la barrera.