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

Utilizando simulações de dinâmica molecular com potenciais de redes neurais, este estudo revela que a migração iônica rápida em MAPbI$_3$ é impulsionada pelo movimento coletivo concertado de intersticiais de MA e defeitos relacionados ao I dependentes de carga, enquanto as vacâncias de MA permanecem imóveis, revisando assim a compreensão convencional do transporte iônico em perovskitas híbridas.

O essencial

Imagine uma célula solar feita de um cristal especial chamado MAPbI3. Pense neste cristal não como um bloco rígido de pedra, mas como uma esponja macia e maleável feita de pequenos blocos de construção. Dentro desta esponja, há dois tipos principais de blocos: blocos de metal pesado (Chumbo e Iodo) e blocos orgânicos mais leves, chamados "moléculas" (chamados MA, que são como pequenas moléculas de metilamônio).

O problema é que esta esponja não é perfeita. Às vezes, blocos desaparecem (criando vacâncias), e às vezes blocos extras são espremidos onde não deveriam estar (criando intersticiais). Quando esses "defeitos" começam a se mover, podem fazer com que a célula solar se degrade ao longo do tempo.

Há muito tempo, cientistas têm tentado descobrir exatamente como esses defeitos se movem e com que velocidade. Os números obtidos em experimentos variavam muito, como um grupo de pessoas adivinhando a velocidade de um carro e obtendo respostas que iam desde "passo de caminhada" até "supersônico".

Este artigo utiliza uma simulação computacional superinteligente (impulsionada por Inteligência Artificial) para observar esses defeitos se movendo em tempo real, como uma câmera de alta velocidade filmando uma pista de dança. Eis o que descobriram, explicado de forma simples:

1. O "Fantasma" vs. O "Transportador Pesado"

Neste cristal, os defeitos de Iodo (os íons haleto) são como fantasmas. São leves e ágeis. Seja um átomo de Iodo faltando (uma vacância) ou um extra espremido (um intersticial), ele se move muito facilmente. A energia necessária para fazê-los mover é muito baixa, como empurrar um carrinho de compras em um piso liso.

2. O Dançarino Surpreendente (A Molécula MA)

A grande surpresa neste artigo envolve as moléculas MA. Elas são muito maiores e mais pesadas do que os átomos de Iodo. Você poderia esperar que fossem lentas, pesadas e difíceis de mover — como tentar empurrar um piano de cauda através de um cômodo.

• A Crença Antiga: Cientistas pensavam que essas grandes moléculas estavam presas ou se moviam muito lentamente.
• A Nova Descoberta: A simulação mostrou que os intersticiais MA (as moléculas extras) são na verdade tão rápidos quanto os fantasmas de Iodo!

Como isso é possível?
O artigo explica que essas grandes moléculas não se movem sozinhas. Elas se movem em um abraço coletivo. Imagine três pessoas em uma pista de dança. Em vez de uma pessoa tentar espremer-se para passar, todas giram e se deslocam juntas em um movimento coordenado e "concertado". Uma dá um passo à frente, as outras giram para abrir espaço e, de repente, todo o grupo se deslocou. Esse trabalho em equipe permite que as pesadas moléculas MA se movam quase tão rápido quanto os minúsculos átomos de Iodo.

3. Aquele Que Fica Parado

Há uma exceção: Vacâncias MA (buracos onde uma molécula MA está faltando). A simulação mostrou que esses buracos são essencialmente imóveis. Mesmo quando a temperatura foi aumentada bastante na simulação, esses buracos não se moveram. É como se o buraco estivesse colado ao chão. Isso sugere que, se você ver o MA se movendo em uma célula solar, provavelmente são as moléculas extras se movendo, e não os espaços vazios.

4. Por Que Os Números Estavam Confusos

O artigo sugere que a razão pela qual experimentos passados deram respostas tão diferentes (alguns dizendo que é lento, outros dizendo que é rápido) é porque estavam medindo coisas diferentes.

• O movimento rápido (barreira de energia de 0,15–0,20 eV) é o que acontece no interior profundo do cristal (difusão volumétrica), que é o foco deste estudo.
• O movimento mais lento relatado em outros estudos pode estar ocorrendo nas bordas dos grãos do cristal ou nas fronteiras entre eles, onde as coisas ficam presas e se movem de maneira diferente.

O Quadro Geral

Este estudo reescreve o manual de regras sobre como entendemos esses materiais. Ele nos diz que:

1. Trabalho em equipe importa: Mesmo moléculas grandes e pesadas podem se mover rápido se se moverem juntas em uma dança coordenada.
2. Carga não importa muito: Ao contrário dos defeitos de Iodo, cuja velocidade muda dependendo de sua carga elétrica, as moléculas MA se movem na mesma velocidade, estejam elas carregadas ou neutras.
3. O "MA lento" é um mito: A ideia de que a parte orgânica do cristal é um gargalo lento e pesado está errada; na verdade, é bastante ágil quando se move em equipe.

Ao entender que esses defeitos são tão móveis e se movem de maneiras específicas, os cientistas agora podem projetar melhores formas de "passivar" (entupir) esses defeitos ou impedir que se movam, o que deve ajudar a fazer células solares e luzes durarem muito mais tempo.

Resumo técnico

Resumo Técnico: Uma Imagem Microscópica Unificada da Migração de Cátions e Ânions em MAPbI3

Enunciação do Problema
A passivação de defeitos e a restrição da mobilidade de defeitos são críticas para estender a vida útil de dispositivos de perovskitas de haleto. Embora estudos experimentais tenham medido extensivamente o transporte iônico no material protótipo MAPbI$_3$ (iodeto de chumbo e metilamônio), os resultados exibem uma dispersão significativa nas barreiras de migração reportadas (variando de 0,10 a 0,94 eV). Além disso, a natureza específica e a composição dos defeitos migrantes — sejam eles vacâncias, intersticiais ou defeitos estendidos — permanecem difíceis de identificar experimentalmente. Abordagens computacionais convencionais, que dependem do método da Banda Elástica Empurrada (NEB), enfrentam dificuldades nestes materiais moles devido à existência de múltiplos caminhos de baixa energia concorrentes e à baixa simetria do cátion metilamônio (MA). Adicionalmente, simulações padrão de Dinâmica Molecular (DM) frequentemente carecem da precisão necessária para a extração confiável de barreiras, devido a limitações nos campos de força clássicos.

Metodologia
Para abordar esses desafios, os autores empregaram simulações de Dinâmica Molecular (DM) impulsionadas por Campos de Força Aprendidos por Máquina (MLFFs) treinados em dados de Teoria do Funcional da Densidade (DFT).

• Desenvolvimento do Campo de Força: O estudo utilizou o modelo Allegro, um potencial de rede neural (NNP) local equivariante a E(3) baseado em grafos, implementado através do pacote NequIP.
• Dados de Treinamento: Os NNPs foram treinados em dados DFT gerados usando o funcional SCAN+rVV10 (incluindo interações de van der Waals) dentro do pacote VASP. Os conjuntos de dados foram construídos utilizando uma abordagem de amostragem "on-the-fly" a partir de execuções de DM de curto prazo, garantindo a inclusão de estruturas diversas.
• Modelagem de Defeitos: As simulações modelaram defeitos pontuais em supercélulas 2$\times$2$\times$2 (96 átomos) e execuções de produção em supercélulas 6$\times$6$\times$6 (2592 átomos). O estudo considerou vacâncias e intersticiais de iodeto (I) e MA em seus estados de carga mais proeminentes sob condições intrínsecas: $I^-_I$, $V^+_I$, $I^0_{MA}$, $I^+_{MA}$, $V^0_{MA}$ e $V^-_{MA}$. NNPs separados foram treinados para defeitos de MA neutros e carregados, selecionando dados DFT com níveis de Fermi próximos às bandas de valência ou condução.
• Validação: A precisão dos NNPs foi validada comparando barreiras de migração calculadas via NEB de Imagem de Escalada (CI-NEB) com benchmarks DFT, mostrando desvios inferiores a 0,11 eV. Comparações de força resultaram em um coeficiente de determinação ($R^2$) de 0,99.
• Protocolo de Simulação: Três simulações de DM independentes de 2 nanosegundos foram realizadas em cinco temperaturas entre 500 K e 600 K. Coeficientes de difusão foram extraídos usando a relação de Einstein, tratando todos os átomos de uma espécie específica para contabilizar processos de "kick-out" em vez de difusão de traçador.

Principais Resultados

• Barreiras de Migração e Mobilidade: As simulações revelaram que intersticiais de iodeto ($I^-_I$) e vacâncias ($V^+_I$) são altamente móveis, com energias de ativação ($E_a$) de 0,18 eV e 0,15 eV, respectivamente. Remarkavelmente, intersticiais de MA ($I^0_{MA}$ e $I^+_{MA}$) exibem mobilidade comparável, com energias de ativação de 0,18 eV e 0,20 eV.
• Vacância de MA Imóvel: Em contraste com os intersticiais, as vacâncias de MA ($V^0_{MA}$ e $V^-_{MA}$) foram encontradas como imóveis dentro das escalas de temperatura e tempo simuladas (até 600 K). Isso sugere um limite inferior para sua energia de ativação de aproximadamente 0,4 eV.
• Mecanismos de Migração:
  • Defeitos de Iodeto: A migração ocorre via um mecanismo de "kick-out", onde um átomo de iodeto vizinho gira para preencher uma vacância, ou um intersticial salta entre pontes Pb-I-Pb.
  • Intersticiais de MA: Apesar de sua natureza molecular, os intersticiais de MA migram via um movimento concertado envolvendo três íons de MA vizinhos. Dois íons de MA em uma gaiola Pb-I e um terceiro vizinho giram para alinhar seus eixos C-N, permitindo que o íon central migre para uma gaiola adjacente.
• Dependência do Estado de Carga: A difusão de defeitos relacionados ao iodeto mostra dependência do estado de carga, enquanto a difusão de defeitos relacionados ao MA é virtualmente independente do estado de carga. Isso é atribuído ao fato de que intersticiais e vacâncias de MA criam níveis de impureza muito rasos que são ionizados à temperatura ambiente, levando a uma carga deslocalizada que não impede significativamente o movimento.
• Consistência com Experimentos: As barreiras de ativação calculadas (0,15–0,20 eV) alinham-se com a extremidade inferior da ampla faixa de valores reportados experimentalmente.

Significado e Alegações
O artigo alega fornecer uma imagem microscópica unificada da migração iônica em MAPbI$_3$ que reconcilia as discrepâncias nas energias de ativação experimentais. Os autores argumentam que os valores mais baixos observados em experimentos (0,15–0,20 eV) correspondem à difusão intrínseca do volume dominada por defeitos de iodeto e intersticiais de MA. Os valores mais altos reportados em outros estudos provavelmente refletem contribuições adicionais, como transporte entre grãos ou a formação de defeitos complexos, em vez de barreiras de migração intrínsecas do volume.

Uma contribuição primária é a revisão da compreensão convencional de que a migração de MA é intrinsecamente lenta devido ao seu tamanho molecular. O estudo demonstra que os intersticiais de MA são quase tão móveis quanto os defeitos de iodeto devido a um mecanismo de migração coletivo e concertado. Além disso, o trabalho destaca que a espécie de cátion A (MA) não afeta significativamente o movimento do haleto, uma vez que a mobilidade de defeitos de haleto em MAPbI$_3$ é comparável à de perovskitas totalmente inorgânicas CsPb(I$_x$Br$_{1-x}$)$_3$. Finalmente, o estudo sublinha a importância do movimento molecular coletivo para permitir a migração iônica rápida em perovskitas híbridas.