A Unified microscopic picture of cation and anion migration in MAPbI$_3$

Mediante simulaciones de dinámica molecular con potenciales de redes neuronales, este estudio revela que la migración iónica rápida en MAPbI$_3$ está impulsada por un movimiento colectivo concertado de intersticiales de MA y defectos relacionados con I dependientes de la carga, mientras que las vacantes de MA permanecen inmóviles, redefiniendo así la comprensión convencional del transporte iónico en perovskitas híbridas.

Lo esencial

Imagina una celda solar hecha de un cristal especial llamado MAPbI3. Piensa en este cristal no como un bloque rígido de piedra, sino como una esponja suave y elástica hecha de pequeños bloques de construcción. Dentro de esta esponja, hay dos tipos principales de bloques: bloques de metales pesados (Plomo y Yodo) y bloques orgánicos más ligeros, llamados "moléculas" (denominados MA, que son como pequeñas moléculas de metilamonio).

El problema es que esta esponja no es perfecta. A veces faltan bloques (creando vacantes) y, a veces, bloques extra se aprietan donde no pertenecen (creando intersticiales). Cuando estos "defectos" comienzan a moverse, pueden hacer que la celda solar se degrade con el tiempo.

Durante mucho tiempo, los científicos han intentado averiguar exactamente cómo se mueven estos defectos y a qué velocidad lo hacen. Los números que obtuvieron de los experimentos variaban enormemente, como un grupo de personas adivinando la velocidad de un coche y obteniendo respuestas que iban desde "ritmo de caminata" hasta "supersónica".

Este artículo utiliza una simulación informática superinteligente (potenciada por Inteligencia Artificial) para observar cómo se mueven estos defectos en tiempo real, como una cámara de alta velocidad filmando una pista de baile. Esto es lo que encontraron, explicado de forma sencilla:

1. El "Fantasma" vs. El "Transportista Pesado"

En este cristal, los defectos de Yodo (los iones haluro) son como fantasmas. Son ligeros y rápidos. Ya sea que falte un átomo de Yodo (una vacante) o que haya uno extra apretado (un intersticial), se desplazan con gran facilidad. La energía necesaria para ponerlos en movimiento es muy baja, como empujar un carrito de la compra sobre un suelo liso.

2. El Bailarín Sorprendente (La Molécula MA)

La gran sorpresa de este artículo involucra a las moléculas MA. Son mucho más grandes y pesadas que los átomos de Yodo. Podrías esperar que fueran lentas, torpes y difíciles de mover, como intentar empujar un piano de cola a través de una habitación.

• La Vieja Creencia: Los científicos pensaban que estas grandes moléculas estaban atascadas o se movían muy lentamente.
• El Nuevo Descubrimiento: La simulación mostró que los intersticiales MA (las moléculas extra) son en realidad tan rápidos como los fantasmas de Yodo.

¿Cómo es esto posible?
El artículo explica que estas grandes moléculas no se mueven solas. Se mueven en un abrazo grupal. Imagina a tres personas en una pista de baile. En lugar de que una persona intente apretarse para pasar, todas rotan y se desplazan juntas en un movimiento coordinado y "concertado". Uno da un paso adelante, los demás rotan para hacer espacio, y de repente todo el grupo se ha desplazado. Este trabajo en equipo permite que las pesadas moléculas MA se desplacen casi tan rápido como los diminutos átomos de Yodo.

3. El Que Se Queda Fijo

Hay una excepción: las Vacantes MA (huecos donde falta una molécula MA). La simulación mostró que estos huecos son esencialmente inmóviles. Incluso cuando la temperatura se elevó mucho en la simulación, estos huecos no se movieron. Es como si el hueco estuviera pegado al suelo. Esto sugiere que si ves movimiento de MA en una celda solar, es probable que sean las moléculas extra las que se mueven, no los espacios vacíos.

4. Por Qué Los Números Estaban Confusos

El artículo sugiere que la razón por la que los experimentos pasados dieron respuestas tan diferentes (algunos diciendo que es lento, otros que es rápido) es porque estaban midiendo cosas distintas.

• El movimiento rápido (barrera de energía de 0.15–0.20 eV) es lo que ocurre en el interior profundo del cristal (difusión en volumen), que es en lo que se centró este estudio.
• El movimiento más lento reportado en otros estudios podría estar ocurriendo en los bordes de los granos del cristal o en los límites entre ellos, donde las cosas se atascan y se mueven de manera diferente.

El Cuadro General

Este estudio reescribe el reglamento de cómo entendemos estos materiales. Nos dice que:

1. El trabajo en equipo importa: Incluso las moléculas grandes y pesadas pueden moverse rápido si se desplazan juntas en una danza coordinada.
2. La carga no importa mucho: A diferencia de los defectos de Yodo, cuya velocidad cambia según su carga eléctrica, las moléculas MA se mueven a la misma velocidad tanto si están cargadas como si son neutras.
3. El "lento" MA es un mito: La idea de que la parte orgánica del cristal es un cuello de botella lento y perezoso es incorrecta; en realidad es bastante ágil cuando se mueve en equipo.

Al comprender que estos defectos son tan móviles y se mueven de formas específicas, los científicos ahora pueden diseñar mejor formas de "pasivar" (taponar) estos defectos o evitar que se muevan, lo cual debería ayudar a que las celdas solares y las luces duren mucho más.

Resumen técnico

Resumen Técnico: Una Imagen Microscópica Unificada de la Migración de Cationes y Aniones en MAPbI3

Planteamiento del Problema
La pasivación de defectos y la restricción de la movilidad de los defectos son críticas para extender la vida útil de los dispositivos de perovskitas de haluro. Si bien los estudios experimentales han medido extensamente el transporte iónico en el material prototipo MAPbI$_3$ (yoduro de plomo y metilamonio), los resultados exhiben una dispersión significativa en las barreras de migración reportadas (que oscilan entre 0.10 y 0.94 eV). Además, la naturaleza específica y la composición de los defectos migratorios —ya sean vacancias, intersticiales o defectos extendidos— siguen siendo difíciles de identificar experimentalmente. Los enfoques computacionales convencionales que dependen del método de la Banda Elástica Empujada (NEB) tienen dificultades en estos materiales blandos debido a la existencia de múltiples trayectorias de baja energía en competencia y la baja simetría del catión metilamonio (MA). Además, las simulaciones estándar de Dinámica Molecular (DM) a menudo carecen de la precisión necesaria para la extracción fiable de barreras debido a las limitaciones de los campos de fuerza clásicos.

Metodología
Para abordar estos desafíos, los autores emplearon simulaciones de Dinámica Molecular (DM) impulsadas por Campos de Fuerza Aprendidos por Máquina (MLFF) entrenados con datos de la Teoría del Funcional de la Densidad (DFT).

• Desarrollo del Campo de Fuerza: El estudio utilizó el modelo Allegro, un potencial de red neuronal (NNP) local equivariante bajo E(3) basado en redes neuronales de grafos, implementado mediante el paquete NequIP.
• Datos de Entrenamiento: Los NNP fueron entrenados con datos DFT generados utilizando el funcional SCAN+rVV10 (incluyendo interacciones de van der Waals) dentro del paquete VASP. Los conjuntos de datos se construyeron utilizando un enfoque de muestreo "sobre la marcha" a partir de ejecuciones de DM de corto tiempo, asegurando la inclusión de estructuras diversas.
• Modelado de Defectos: Las simulaciones modelaron defectos puntuales en superceldas de 2$\times$2$\times$2 (96 átomos) y ejecuciones de producción en superceldas de 6$\times$6$\times$6 (2592 átomos). El estudio consideró vacancias e intersticiales de yoduro (I) y MA en sus estados de carga más prominentes bajo condiciones intrínsecas: $I^-_I$, $V^+_I$, $I^0_{MA}$, $I^+_{MA}$, $V^0_{MA}$ y $V^-_{MA}$. Se entrenaron NNP separados para defectos de MA neutros y cargados seleccionando datos DFT con niveles de Fermi cercanos a las bandas de valencia o conducción.
• Validación: La precisión de los NNP se validó comparando las barreras de migración calculadas mediante NEB de Imagen de Escalada (CI-NEB) contra referencias DFT, mostrando desviaciones menores a 0.11 eV. Las comparaciones de fuerzas arrojaron un coeficiente de determinación ($R^2$) de 0.99.
• Protocolo de Simulación: Se realizaron tres simulaciones de DM independientes de 2 nanosegundos a cinco temperaturas entre 500 K y 600 K. Los coeficientes de difusión se extrajeron utilizando la relación de Einstein, tratando todos los átomos de una especie específica para dar cuenta de los procesos de expulsión ("kick-out") en lugar de la difusión por trazadores.

Resultados Clave

• Barreras de Migración y Movilidad: Las simulaciones revelaron que los intersticiales de yoduro ($I^-_I$) y las vacancias ($V^+_I$) son altamente móviles, con energías de activación ($E_a$) de 0.18 eV y 0.15 eV, respectivamente. Remarkablemente, los intersticiales de MA ($I^0_{MA}$ e $I^+_{MA}$) exhiben una movilidad comparable con energías de activación de 0.18 eV y 0.20 eV.
• Vacancia de MA Inmóvil: En contraste con los intersticiales, se encontró que las vacancias de MA ($V^0_{MA}$ y $V^-_{MA}$) son inmóviles dentro de las escalas de tiempo y temperatura simuladas (hasta 600 K). Esto sugiere un límite inferior para su energía de activación de aproximadamente 0.4 eV.
• Mecanismos de Migración:
  • Defectos de Yoduro: La migración ocurre mediante un mecanismo de "expulsión" donde un átomo de yoduro vecino gira para llenar una vacancia, o un intersticial salta entre puentes Pb-I-Pb.
  • Intersticiales de MA: A pesar de su naturaleza molecular, los intersticiales de MA migran mediante un movimiento concertado que involucra a tres iones de MA vecinos. Dos iones de MA en una jaula Pb-I y un tercer vecino giran para alinear sus ejes C-N, permitiendo que el ion central migre a una jaula adyacente.
• Dependencia del Estado de Carga: La difusión de defectos relacionados con el yoduro muestra dependencia del estado de carga, whereas la difusión de defectos relacionados con el MA es virtualmente independiente del estado de carga. Esto se atribuye al hecho de que los intersticiales y las vacancias de MA crean niveles de impureza muy superficiales que están ionizados a temperatura ambiente, lo que conduce a una carga deslocalizada que no obstaculiza significativamente el movimiento.
• Consistencia con Experimentos: Las barreras de activación calculadas (0.15–0.20 eV) se alinean con el extremo inferior del amplio rango de valores reportados experimentalmente.

Significado y Afirmaciones
El artículo afirma proporcionar una imagen microscópica unificada de la migración iónica en MAPbI$_3$ que reconcilia las discrepancias en las energías de activación experimentales. Los autores argumentan que los valores más bajos observados en experimentos (0.15–0.20 eV) corresponden a la difusión intrínseca del volumen dominada por defectos de yoduro e intersticiales de MA. Los valores más altos reportados en otros estudios probablemente reflejan contribuciones adicionales, como el transporte entre granos o la formación de defectos complejos, en lugar de barreras de migración intrínsecas del volumen.

Una contribución principal es la revisión de la comprensión convencional de que la migración de MA es intrínsecamente lenta debido a su tamaño molecular. El estudio demuestra que los intersticiales de MA son casi tan móviles como los defectos de yoduro debido a un mecanismo de migración colectivo y concertado. Además, el trabajo destaca que la especie de catión A (MA) no afecta significativamente el movimiento de los haluros, ya que la movilidad de los defectos de haluro en MAPbI$_3$ es comparable a la de las perovskitas totalmente inorgánicas CsPb(I$_x$Br$_{1-x}$)$_3$. Finalmente, el estudio subraya la importancia del movimiento molecular colectivo para permitir una rápida migración iónica en perovskitas híbridas.