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

本研究は、ニューラルネットワークポテンシャルを用いた分子動力学シミュレーションにより、MAPbI$_3$における急速なイオン移動が、MA 格子間原子と電荷依存性の I 関連欠陥の協調的な集団運動によって駆動され、MA 空孔は移動しないことを明らかにし、これによりハイブリッドペロブスカイトにおけるイオン輸送に関する従来の理解を修正するものである。

要約

MAPbI3という特殊な結晶でできた太陽電池を想像してみてください。この結晶を、硬い石のブロックではなく、小さな構成要素でできた柔らかくしなやかなスポンジだと考えてください。このスポンジの中には、2 種類の主なブロックがあります。重い金属ブロック（鉛とヨウ素）と、より軽く有機的な「分子」ブロック（MA と呼ばれ、小さなメチルアミン分子のようなもの）です。

問題は、このスポンジが完璧ではないことです。時にはブロックが欠けて（空孔が生まれ）、時には余分なブロックが居場所のないところに押し込まれて（格子間が生まれます）。これらの「欠陥」が動き回ると、太陽電池は時間とともに劣化してしまう可能性があります。

長らく、科学者たちはこれらの欠陥がどのように移動し、どれだけの速さで動くのかを解明しようと試みてきました。実験から得られた数値はバラバラで、まるで人々が車の速さを推測して、「歩行ペース」から「超音速」まで様々な答えが出たかのようでした。

この論文では、人工知能を駆使した超高性能なコンピュータシミュレーションを用いて、これらの欠陥の動きを高速カメラがダンスフロアを撮影するようにリアルタイムで追跡しました。以下に、その発見を分かりやすく説明します。

1. 「ゴースト」と「重荷」

この結晶において、ヨウ素の欠陥（ハロゲン化物イオン）はゴーストのようです。それらは軽く、素早いです。ヨウ素原子が欠けている場合（空孔）でも、余分なものが押し込まれている場合（格子間）でも、非常に簡単に動き回ります。これらを動かすために必要なエネルギーは非常に低く、滑らかな床でショッピングカートを押すようなものです。

2. 驚きのダンサー（MA 分子）

この論文における最大の驚きは、MA 分子に関わるものです。これらはヨウ素原子よりもはるかに大きく、重いです。そのため、ゆっくりと鈍重で、動かすのが難しいだろうと予想されるかもしれません。まるで部屋中にグランドピアノを押し動かそうとするようなものです。

• 従来の考え： 科学者たちは、これらの大きな分子は固定されているか、非常にゆっくりしか動かないと考えていました。
• 新たな発見： シミュレーションは、MA 格子間（余分な分子）が、実際にはヨウ素のゴーストと同じくらい速いことを示しました。

これはどうして可能なのでしょうか？
この論文は、これらの大きな分子は単独で動くのではなく、「集団のハグ」のように動くことを説明しています。ダンスフロアに 3 人の人がいると想像してください。1 人が無理やりすり抜けようとするのではなく、全員が協調して回転し、移動します。一人が前に一歩踏み出すと、他の人たちがスペースを作るために回転し、突然、グループ全体が移動します。このチームワークにより、重い MA 分子は、小さなヨウ素原子とほぼ同じ速さで動き回ることができるのです。

3. 動かずに留まる者

例外が 1 つあります。それはMA 空孔（MA 分子が欠けている穴）です。シミュレーションは、これらの穴が本質的に移動しないことを示しました。シミュレーションで温度を高くしても、これらの穴は動きませんでした。まるで穴が床に接着されているかのようです。これは、太陽電池内で MA が移動しているのを見た場合、それは空の場所ではなく、おそらく余分な分子が移動していることを示唆しています。

4. なぜ数値が混乱していたのか

この論文は、過去の研究がこれほど異なる答え（一部は遅いとし、一部は速いとする）を出した理由は、異なるものを測定していたからだと示唆しています。

• 速い移動（0.15–0.20 eV のエネルギー障壁）は、結晶の奥深く（バルク拡散）で起こる現象であり、この研究が焦点を当てた部分です。
• 他の研究で報告されたより遅い移動は、結晶粒の端や粒界で起こっている可能性があり、そこでは物が詰まって異なった動き方をしているのかもしれません。

全体像

この研究は、これらの材料の理解に関するルールブックを書き換えるものです。それは以下を伝えています。

1. チームワークが重要である： 大きく重い分子でも、協調したダンスのように一緒に動けば、速く移動できる。
2. 電荷はあまり重要ではない： 電荷によって速度が変化するヨウ素の欠陥とは異なり、MA 分子は電荷を持っていようが中性であろうが、同じ速度で移動する。
3. 「遅い」MA は神話である： 結晶の有機部分が遅く、鈍重なボトルネックであるという考えは誤りであり、実際にはチームとして動くときは非常に機敏である。

これらの欠陥が非常に移動性が高く、特定の方法で動くことを理解することで、科学者たちはこれらの欠陥を「パッシベーション（封止）」したり、その移動を止めたりするより良い方法を設計できるようになり、太陽電池や照明の寿命を大幅に延ばすことが期待されます。

技術概要

技術的概要：MAPbI3 におけるカチオンおよびアニオン移動の統合された微視的描像

問題提起
欠陥のパッシベーションと欠陥移動性の制限は、ハライドペロブスカイトデバイスの寿命を延伸する上で決定的に重要である。実験的研究は、プロトタイプ材料である MAPbI$_3$（メチルアンモニウム鉛ヨウ化物）におけるイオン輸送を広く測定してきたが、報告された移動障壁の結果には大きなばらつき（0.10 eV から 0.94 eV の範囲）が見られる。さらに、移動する欠陥の具体的な性質と組成（空孔、格子間原子、あるいは拡張欠陥のいずれかであるか）は、実験的に同定することが依然として困難である。従来の計算アプローチは、ヌジッド・エラスティック・バンド（NEB）法に依存しているが、複数の競合する低エネルギー経路の存在と、メチルアンモニウム（MA）カチオンの低対称性のため、このような軟物質においては困難に直面する。加えて、標準的な分子動力学（MD）シミュレーションは、古典的力場の限界により、信頼性の高い障壁抽出に必要な精度を欠くことが多い。

手法
これらの課題に対処するため、著者らは密度汎関数理論（DFT）データに基づいて訓練された機械学習力場（MLFFs）によって駆動される分子動力学（MD）シミュレーションを採用した。

• 力場の開発: 本研究では、NequIP パッケージを介して実装された、局所的な E(3)-共変グラフニューラルネットワークポテンシャル（NNP）であるAllegroモデルを利用した。
• 訓練データ: NNP は、VASP パッケージ内で SCAN+rVV10 汎関数（ファンデルワールス相互作用を含む）を用いて生成された DFT データに基づいて訓練された。データセットは、短い時間の MD 実行からの「オン・ザ・フライ」サンプリング手法を用いて構築され、多様な構造の包含が保証された。
• 欠陥モデリング: シミュレーションでは、2$\times$2$\times$2 超格子（96 原子）における点欠陥と、6$\times$6$\times$6 超格子（2592 原子）における生産実行をモデル化した。本研究では、固有条件下で最も顕著な電荷状態にあるヨウ化物（I）および MA 空孔と格子間原子を考慮した：$I^-_I$、$V^+_I$、$I^0_{MA}$、$I^+_{MA}$、$V^0_{MA}$、および$V^-_{MA}$。中性および帯電した MA 欠陥に対しては、価電子帯または伝導帯付近のフェルミ準位を持つ DFT データを選択することで、個別の NNP を訓練した。
• 検証: NNP の精度は、クライミング・イメージ NEB（CI-NEB）を介して計算された移動障壁を DFT ベンチマークと比較することで検証され、0.11 eV 未満の偏差を示した。力の比較からは、決定係数（$R^2$）が 0.99 となった。
• シミュレーションプロトコル: 500 K から 600 K の間の 5 つの温度において、3 つの独立した 2 ナノ秒の MD シミュレーションが実施された。拡散係数は、トレーサー拡散ではなくキックアウト過程を考慮するために、特定の種に属するすべての原子を扱うアインシュタイン関係を用いて抽出された。

主要な結果

• 移動障壁と移動性: シミュレーションは、ヨウ化物格子間原子（$I^-_I$）および空孔（$V^+_I$）が非常に高い移動性を有し、それぞれ活性化エネルギー（$E_a$）が 0.18 eV および 0.15 eV であることを明らかにした。驚くべきことに、MA 格子間原子（$I^0_{MA}$および$I^+_{MA}$）も、活性化エネルギーがそれぞれ 0.18 eV および 0.20 eV であり、同様の移動性を示す。
• MA 空孔の不動性: 格子間原子とは対照的に、MA 空孔（$V^0_{MA}$および$V^-_{MA}$）は、シミュレーションされた温度および時間スケール（最大 600 K）において不動であることが判明した。これは、それらの活性化エネルギーの下限が約 0.4 eV であることを示唆する。
• 移動メカニズム:
  • ヨウ化物欠陥: 移動は、隣接するヨウ化物原子が回転して空孔を埋めるか、あるいは格子間原子が Pb-I-Pb ブリッジ間をホップする「キックアウト」メカニズムを介して起こる。
  • MA 格子間原子: 分子性であるにもかかわらず、MA 格子間原子は、3 つの隣接する MA イオンを含む協調的運動を介して移動する。Pb-I ケージ内の 2 つの MA イオンと 3 つ目の隣接イオンが回転して C-N 軸を揃え、中央のイオンが隣接するケージへ移動することを可能にする。
• 電荷状態への依存性: ヨウ化物関連欠陥の拡散は電荷状態に依存するが、MA 関連欠陥の拡散は電荷状態に実質的に依存しない。これは、MA 格子間原子および空孔が室温で電離する非常に浅い不純物準位を形成し、運動を著しく妨げない非局在化電荷をもたらすことに起因する。
• 実験との整合性: 計算された活性化障壁（0.15–0.20 eV）は、実験的に報告された広範な値の下限と一致する。

意義と主張
本論文は、実験的な活性化エネルギーの不一致を調和させる、MAPbI$_3$におけるイオン移動の統合された微視的描像を提供すると主張している。著者らは、実験で観測される低い値（0.15–0.20 eV）は、ヨウ化物欠陥および MA 格子間原子によって支配される固有のバルク拡散に対応すると論じている。他の研究で報告された高い値は、固有のバルク移動障壁ではなく、粒間輸送や複合欠陥の形成などの追加的な寄与を反映している可能性が高い。

主要な貢献の一つは、分子サイズのために MA 移動が本質的に遅いという従来の理解の改訂である。本研究は、MA 格子間原子が集合的かつ協調的な移動メカニズムにより、ヨウ化物欠陥とほぼ同等の移動性を有することを示している。さらに、この研究は、A 位カチオン種（MA）がハライド運動に著しく影響を与えないことを強調しており、MAPbI$_3$におけるハライド欠陥の移動性は、無機ペロブスカイト CsPb(I$_x$Br$_{1-x}$)$_3$におけるそれと同等である。最後に、本研究は、ハイブリッドペロブスカイトにおける高速イオン移動を可能にする上で、集合的分子運動の重要性を浮き彫りにしている。