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

Mittels Molekulardynamiksimulationen mit neuronalen Netzwerkpotentialen zeigt diese Studie, dass die schnelle Ionenmigration in MAPbI$_3$ durch eine konzertierte kollektive Bewegung von MA-Interstitien und ladungsabhängigen I-bezogenen Defekten angetrieben wird, während MA-Leerstellen unbeweglich bleiben, wodurch das konventionelle Verständnis des ionischen Transports in hybriden Perowskiten revidiert wird.

Das Wesentliche

Stellen Sie sich eine Solarzelle vor, die aus einem speziellen Kristall namens MAPbI3 besteht. Betrachten Sie diesen Kristall nicht als starren Steinblock, sondern als einen weichen, quetschbaren Schwamm, der aus winzigen Bausteinen besteht. In diesem Schwamm gibt es zwei Haupttypen von Bausteinen: schwere Metallbausteine (Blei und Jod) und leichtere, organische „Molekül"-Bausteine (genannt MA, die wie kleine Methylammonium-Moleküle sind).

Das Problem ist, dass dieser Schwamm nicht perfekt ist. Manchmal fehlen Bausteine (was Leerstellen erzeugt), und manchmal werden zusätzliche Bausteine dort hineingedrückt, wo sie nicht hingehören (was Zwischengitteratome erzeugt). Wenn sich diese „Defekte" herumzubewegen beginnen, können sie dazu führen, dass die Solarzelle im Laufe der Zeit zusammenbricht.

Lange Zeit haben Wissenschaftler versucht herauszufinden, wie genau sich diese Defekte bewegen und wie schnell sie sind. Die Zahlen, die sie aus Experimenten erhielten, waren völlig unterschiedlich, wie eine Gruppe von Menschen, die die Geschwindigkeit eines Autos schätzen und Antworten von „Schrittempo" bis „überschall" erhalten.

Diese Arbeit verwendet eine hochintelligere Computersimulation (angetrieben durch Künstliche Intelligenz), um diese Defekte in Echtzeit zu beobachten, wie eine Hochgeschwindigkeitskamera, die einen Tanzboden filmt. Hier ist das Ergebnis, einfach erklärt:

1. Der „Geist" gegen den „Schwertransporter"

In diesem Kristall sind die Jod-Defekte (die Halogenid-Ionen) wie Geister. Sie sind leicht und flink. Ob ein Jod-Atom fehlt (eine Leerstelle) oder ein zusätzliches hineingedrückt wurde (ein Zwischengitteratom), es rast sehr leicht herum. Die Energie, die benötigt wird, um sie in Bewegung zu setzen, ist sehr gering, wie das Schieben eines Einkaufswagens auf einem glatten Boden.

2. Der überraschende Tänzer (das MA-Molekül)

Die große Überraschung in dieser Arbeit betrifft die MA-Moleküle. Diese sind viel größer und schwerer als die Jod-Atome. Man könnte erwarten, dass sie langsam, schwerfällig und schwer zu bewegen sind – wie der Versuch, ein Flügelklavier durch einen Raum zu schieben.

• Der alte Glaube: Wissenschaftler glaubten, dass diese großen Moleküle feststecken oder sich sehr langsam bewegen.
• Die neue Entdeckung: Die Simulation zeigte, dass die MA-Zwischengitteratome (die zusätzlichen Moleküle) tatsächlich genauso schnell sind wie die Jod-Geister!

Wie ist das möglich?
Die Arbeit erklärt, dass diese großen Moleküle nicht allein bewegen. Sie bewegen sich in einer Gruppenumarmung. Stellen Sie sich drei Personen auf einem Tanzboden vor. Anstatt dass eine Person versucht, sich hindurchzudrängen, drehen und verschieben sie sich alle gemeinsam in einer koordinierten, „konzertierten" Bewegung. Einer macht einen Schritt vorwärts, die anderen drehen sich, um Platz zu machen, und plötzlich hat sich die ganze Gruppe verschoben. Diese Teamarbeit ermöglicht es den schweren MA-Molekülen, fast genauso schnell herumzurennen wie die winzigen Jod-Atome.

3. Der, der stillsteht

Es gibt eine Ausnahme: MA-Leerstellen (Löcher, in denen ein MA-Molekül fehlt). Die Simulation zeigte, dass diese Löcher im Wesentlichen unbeweglich sind. Selbst als die Temperatur in der Simulation stark erhöht wurde, bewegten sich diese Löcher nicht. Es ist, als wäre das Loch am Boden festgeklebt. Dies deutet darauf hin, dass, wenn man MA-Bewegung in einer Solarzelle sieht, es wahrscheinlich die zusätzlichen Moleküle sind, die sich bewegen, und nicht die leeren Stellen.

4. Warum die Zahlen verwirrend waren

Die Arbeit legt nahe, dass der Grund, warum frühere Experimente so unterschiedliche Antworten lieferten (einige sagten, es sei langsam, andere sagten, es sei schnell), darin liegt, dass sie unterschiedliche Dinge gemessen haben.

• Die schnelle Bewegung (Energiebarriere von 0,15–0,20 eV) ist das, was tief im Inneren des Kristalls passiert (Volumendiffusion), worauf sich diese Studie konzentrierte.
• Die langsamere Bewegung, die in anderen Studien berichtet wurde, findet möglicherweise an den Rändern der Kristallkörner oder an den Grenzen zwischen ihnen statt, wo Dinge stecken bleiben und sich anders bewegen.

Das große Ganze

Diese Studie schreibt das Regelbuch darüber um, wie wir diese Materialien verstehen. Sie sagt uns, dass:

1. Teamarbeit zählt: Selbst große, schwere Moleküle können sich schnell bewegen, wenn sie sich gemeinsam in einem koordinierten Tanz bewegen.
2. Ladung spielt keine große Rolle: Im Gegensatz zu den Jod-Defekten, deren Geschwindigkeit sich je nach ihrer elektrischen Ladung ändert, bewegen sich die MA-Moleküle mit der gleichen Geschwindigkeit, egal ob sie geladen oder neutral sind.
3. Der „langsame" MA ist ein Mythos: Die Idee, dass der organische Teil des Kristalls eine langsame, träge Engstelle ist, ist falsch; er ist tatsächlich ziemlich agil, wenn er als Team bewegt wird.

Indem man versteht, dass diese Defekte so beweglich sind und sich auf bestimmte Weise bewegen, können Wissenschaftler nun bessere Wege entwickeln, um diese Defekte zu „passivieren" (zu stopfen) oder zu verhindern, dass sie sich bewegen, was dazu beitragen sollte, Solarzellen und Leuchten viel länger haltbar zu machen.

Technische Zusammenfassung

Technische Zusammenfassung: Ein vereinheitlichtes mikroskopisches Bild der Kationen- und Anionenmigration in MAPbI₃

Problemstellung
Die Defektpassivierung und die Einschränkung der Defektmobilität sind entscheidend für die Verlängerung der Lebensdauer von Halogenid-Perowskit-Bauelementen. Während experimentelle Studien den Ionentransport im Prototyp-Material MAPbI₃ (Methylammoniumbleijodid) umfassend gemessen haben, zeigen die Ergebnisse eine erhebliche Streuung der berichteten Migrationsbarrieren (im Bereich von 0,10 bis 0,94 eV). Darüber hinaus bleibt die spezifische Natur und Zusammensetzung der wandernden Defekte – ob es sich um Leerstellen, Zwischengitteratome oder ausgedehnte Defekte handelt – experimentell schwer zu identifizieren. Herkömmliche rechnerische Ansätze, die auf der Nudged-Elastic-Band-Methode (NEB) basieren, stoßen bei diesen weichen Materialien aufgrund des Vorhandenseins mehrerer konkurrierender Pfade mit niedriger Energie und der niedrigen Symmetrie des Methylammonium-Kations (MA) an ihre Grenzen. Zudem fehlt es Standard-Molekulardynamik-Simulationen (MD) aufgrund von Einschränkungen klassischer Kraftfelder oft an der für eine zuverlässige Extraktion von Barrieren erforderlichen Genauigkeit.

Methodik
Um diese Herausforderungen zu bewältigen, setzten die Autoren Molekulardynamik-Simulationen (MD) ein, die durch auf Dichtefunktionaltheorie (DFT)-Daten trainierte maschinell erlernte Kraftfelder (MLFFs) angetrieben wurden.

• Entwicklung des Kraftfelds: Die Studie nutzte das Allegro-Modell, ein lokales E(3)-äquivariantes Graph-Neural-Network-Potential (NNP), implementiert über das NequIP-Paket.
• Trainingsdaten: Die NNPs wurden auf DFT-Daten trainiert, die mit dem SCAN+rVV10-Funktional (einschließlich van-der-Waals-Wechselwirkungen) im VASP-Paket generiert wurden. Die Datensätze wurden mittels eines „On-the-Fly"-Sampling-Ansatzes aus Kurzzeit-MD-Läufen erstellt, um die Einbeziehung vielfältiger Strukturen sicherzustellen.
• Defektmodellierung: Die Simulationen modellierten Punktdefekte in 2×2×2-Supergittern (96 Atome) und Produktionsläufe in 6×6×6-Supergittern (2592 Atome). Die Studie betrachtete Jodid- (I) und MA-Leerstellen sowie Zwischengitteratome in ihren unter intrinsischen Bedingungen vorherrschenden Ladungszuständen: $I^-_I$, $V^+_I$, $I^0_{MA}$, $I^+_{MA}$, $V^0_{MA}$ und $V^-_{MA}$. Für neutrale und geladene MA-Defekte wurden separate NNPs trainiert, indem DFT-Daten mit Fermi-Niveaus nahe den Valenz- oder Leitungsbändern ausgewählt wurden.
• Validierung: Die Genauigkeit der NNPs wurde durch den Vergleich von Migrationsbarrieren validiert, die mittels Climbing-Image-NEB (CI-NEB) berechnet wurden, gegen DFT-Referenzwerte, wobei Abweichungen von weniger als 0,11 eV gezeigt wurden. Kraftvergleiche ergaben einen Bestimmtheitsmaß ($R^2$) von 0,99.
• Simulationsprotokoll: Drei unabhängige 2-Nanosekunden-MD-Simulationen wurden bei fünf Temperaturen zwischen 500 K und 600 K durchgeführt. Diffusionskoeffizienten wurden unter Verwendung der Einstein-Beziehung extrahiert, wobei alle Atome einer bestimmten Spezies berücksichtigt wurden, um Kick-out-Prozesse statt Tracer-Diffusion zu erfassen.

Hauptergebnisse

• Migrationsbarrieren und Mobilität: Die Simulationen zeigten, dass Jodid-Zwischengitteratome ($I^-_I$) und Leerstellen ($V^+_I$) hochmobil sind, mit Aktivierungsenergien ($E_a$) von 0,18 eV bzw. 0,15 eV. Bemerkenswerterweise weisen MA-Zwischengitteratome ($I^0_{MA}$ und $I^+_{MA}$) eine vergleichbare Mobilität mit Aktivierungsenergien von 0,18 eV bzw. 0,20 eV auf.
• MA-Leerstellen unbeweglich: Im Gegensatz zu Zwischengitteratomen erwiesen sich MA-Leerstellen ($V^0_{MA}$ und $V^-_{MA}$) innerhalb der simulierten Temperatur- und Zeitskalen (bis 600 K) als unbeweglich. Dies deutet auf eine untere Grenze ihrer Aktivierungsenergie von etwa 0,4 eV hin.
• Migrationsmechanismen:
  • Jodid-Defekte: Die Migration erfolgt über einen „Kick-out"-Mechanismus, bei dem ein benachbartes Jodid-Atom rotiert, um eine Leerstelle zu füllen, oder ein Zwischengitteratom zwischen Pb-I-Pb-Brücken hüpft.
  • MA-Zwischengitteratome: Trotz ihrer molekularen Natur wandern MA-Zwischengitteratome über eine koordinierte Bewegung, die drei benachbarte MA-Ionen umfasst. Zwei MA-Ionen in einem Pb-I-Käfig und ein dritter Nachbar rotieren, um ihre C-N-Achsen auszurichten, wodurch es dem zentralen Ion ermöglicht wird, in einen benachbarten Käfig zu wandern.
• Ladungszustandsabhängigkeit: Die Diffusion jodidbezogener Defekte zeigt eine Abhängigkeit vom Ladungszustand, während die Diffusion MA-bezogener Defekte praktisch unabhängig vom Ladungszustand ist. Dies wird darauf zurückgeführt, dass MA-Zwischengitteratome und -Leerstellen sehr flache Störstellen-Niveaus erzeugen, die bei Raumtemperatur ionisiert sind, was zu einer delokalisierten Ladung führt, die die Bewegung nicht signifikant behindert.
• Konsistenz mit Experimenten: Die berechneten Aktivierungsbarrieren (0,15–0,20 eV) stimmen mit dem unteren Ende der breiten Spanne experimentell berichteter Werte überein.

Bedeutung und Behauptungen
Die Arbeit behauptet, ein vereinheitlichtes mikroskopisches Bild der Ionmigration in MAPbI₃ zu liefern, das die Diskrepanzen in den experimentellen Aktivierungsenergien auflöst. Die Autoren argumentieren, dass die in Experimenten beobachteten niedrigeren Werte (0,15–0,20 eV) der intrinsischen Volumendiffusion entsprechen, die von Jodid-Defekten und MA-Zwischengitteratomen dominiert wird. Die in anderen Studien berichteten höheren Werte spiegeln wahrscheinlich zusätzliche Beiträge wider, wie z. B. Transport zwischen Körnern oder die Bildung komplexer Defekte, anstatt intrinsische Volumen-Migrationsbarrieren.

Ein wesentlicher Beitrag ist die Überarbeitung des konventionellen Verständnisses, dass die MA-Migration aufgrund ihrer molekularen Größe intrinsisch langsam ist. Die Studie zeigt, dass MA-Zwischengitteratome aufgrund eines kollektiven, koordinierten Migrationsmechanismus nahezu so mobil sind wie Jodid-Defekte. Darüber hinaus hebt die Arbeit hervor, dass die A-Kationen-Spezies (MA) die Halogenidbewegung nicht signifikant beeinflusst, da die Mobilität von Halogenid-Defekten in MAPbI₃ mit der in rein anorganischen CsPb(I$_x$Br$_{1-x}$)$_3$-Perowskiten vergleichbar ist. Schließlich unterstreicht die Studie die Bedeutung kollektiver molekularer Bewegung für die Ermöglichung schneller ionischer Migration in hybriden Perowskiten.