PAOFLOW: an automated suite for ab initio electronic, transport, and topological properties of materials

現代の材料科学は、強力ではあるものの、しばしば扱いにくいツールに依存しています。それは、コンピュータ上で原子をシミュレーションすることで、物質がどのように振る舞うかを予測する能力です。科学者たちは、電気の流れから、特定の鉱物がなぜ超伝導体になり得るのかに至るまで、あらゆることを理解するためにこれらのシミュレーションを使用します。しかし、最も正確なシミュレーションは非常に低速です。それらは、材料の電子を記述する数学的空間内のあらゆる一点に対して、複雑な方程式を解くことをコンピュータに要求します。材料の特性を明確に把握するために、研究者はしばしば数百万のこれらの点をサンプリングする必要がありますが、このプロセスは、最速のスーパーコンピュータを用いても数日または数週間かかることがあります。このボトルネックにより、人工知能（AI）が新材料を発見するために必要な、大規模で一貫したデータベースを構築することが困難になってきました。もしコンピュータが一つの材料を分析するのに時間がかかりすぎれば、次の画期的な発見を見つけるために必要な何千もの候補をスクリーニングすることができなくなります。

研究チームは、精度を損なうことなくこの速度の問題を解決する、PAOFLOW 3.0と呼ばれる新しいソフトウェア・スイートを開発しました。このプログラムは、標準的なコンピュータ・シミュレーションによる重厚で詳細な結果を取り込み、より軽量でコンパクトな形式へと変換する「翻訳者」として機能します。この新しい形式は、材料の正確な電子挙動を保持しながら、科学者が特性をほぼ瞬時に計算することを可能にします。この変換を自動化することで、このソフトウェアは、AIのために膨大な量の高品質なデータを生成し、以前は研究コストが高すぎて研究できなかった複雑な物理現象を探索することを可能にします。

この成果の核心は、電子の動きの数学をソフトウェアがどのように扱うかにあります。標準的なシミュレーションは、結晶全体に広がる波として電子を記述します。これは正確ですが、計算コストが非常に高い波の操作を必要とします。PAOFLOWは、これらの波の記述を取り込み、局在化した原子中心の軌道へと投影します。これは、風景の高解像度写真を、詳細なランドマークの地図へと変換することに似ています。地図はナビゲートや測定がはるかに容易になりますが、地形に関する不可欠な情報はすべて保持されています。この変換によって、「タイトバインディング（強結合）」モデル、すなわち、選択されたエネルギー範囲において数学的に正確な、材料の電子構造の簡略化された表現が作成されます。このモデルが一度構築されれば、ソフトウェアは電子バンドを補間（つまり、元のデータポイント間の隙間を埋めること）し、任意の密度のグリッド上で材料の挙動の滑らかで連続的な図を生成できます。これにより、フェルミ面の形状（電子が材料内をどのように移動するかを決定するもの）など、極めて高い精度を必要とする特性の計算が可能になります。

ソフトウェアの新バージョンであるPAOFLOW 3.0は、これらのモデルを用いて科学者ができることを大幅に拡張しています。現在、世界で最も広く使用されている2つのシミュレーション・コードからの計算をサポートしており、研究者は異なる計算ワークフロー間を、解析手法を変えることなく切り替えることができます。また、遷移金属酸化物のような、強く相互作用する電子を持つ材料を正確に記述するために必要な補正である「ハバードUおよびV」パラメータを算出するための自己整合的な手法を導入しています。以前は、これらのパラメータは実験データに基づいて推測またはフィッティングされることが多かったのですが、この新しいアプローチは、これらを第一原理から直接計算するため、推測の余地を排除します。さらに、このソフトウェアは「環境依存型」モデルを生成します。これは、原子間の結合を記述する数学的ルールが、局所的な環境の混み具合に応じて変化することを意味します。これは、原子配置が均一ではない表面、界面、または材料のねじれた層を研究する上で極めて重要です。

これらの堅牢なモデルを手にした状態で、ソフトウェアは実世界の応用において直接関連する幅広い物理的特性を計算できます。材料がどのように電気や熱を伝導するか、磁場にどのように反応するか、そして光とどのように相互作用するかを予測できます。例えば、磁場が導体内の移動電子を片側に押し出すホール効果や、温度勾配から電圧を発生させるネルンスト効果を決定できます。また、このソフトウェアは電子の量子的な性質である「スピン」の複雑な物理学も扱い、スピン電流や、電流が正味のスピン蓄積を誘起するラシュバ・エデルシュタイン効果を計算することができます。これらの能力は、電荷ではなく電子のスピンを利用して情報を保存・処理することを目指す分野であるスピントロニクスの発展において不可欠です。

標準的な輸送現象を超えて、ソフトウェアはトポロジカル材料の量子領域にも取り組みます。これらは、電子状態の全体的な配置によって、内部は絶縁体でありながら表面では電気を通すといったユニークな挙動が生じる物質です。PAOFLOWは、保護された形でバンドが交差する特定の場所である「ワイル点」を含む、これらのトポロジカルな特徴を特定できます。また、材料のトポロジカルな性質を分類する整数である「チャーン数」や、隠れた対称性を明らかにする「ミラー・チャーン数」を計算することも可能です。このソフトウェアには、磁場が変化するにつれて材料の電気的特性が周期的に振動する現象である「量子振動」の周波数を特定するために特別に設計されたモジュールが含まれています。これらの振動は、材料のフェルミ面の直接的な指紋を提供し、科学者がシミュレーションを実験的な測定値と直接比較することを可能にします。

ソフトウェアの有用性は、電子と原子格子の振動（フォノン）との相互作用の研究にも及びます。この電子－フォノン結合は、電気が抵抗ゼロで流れる従来の超伝導のメカニズムです。ソフトウェアは、これらの相互作用を極めて高密度なグリッド上に補間できますが、これは標準的な手法では禁止的なコストがかかる作業です。これにより、超伝導転移温度（材料が超伝導体になる点）の精密な計算が可能になります。さらに、ソフトウェアは単一原子のワイヤーのようなナノスケールデバイスを通じた電子の輸送をモデル化し、電流がどのように接合部を流れ、欠陥や不純物によってどのように影響を受けるかを計算できます。

その能力を実証するために、研究者らは単純なシリコンから複雑なねじれ二層グラフェンに至るまで、さまざまな材料にこのソフトウェアを適用しました。二層のカーボンシートが相対的に回転して、数千個の原子を含む巨大な超格子を形成するねじれ二層グラフェンの場合、ソフトウェアは、その材料のフラットバンドや低エネルギー物理学を研究するために使用できる、転送可能なモデルの生成に成功しました。これは、システムのあまりの大きさゆえに、直接的なシミュレーションでは不可能でした。また、ソフトウェアはMoP2のような材料のトポロジカルな特性をモデル化してワイル点の正確な位置を特定し、鉛における電子－フォノン結合をモデル化して、確立された結果を高精度で再現することにも成功しました。

ソフトウェアのアーキテクチャは、パフォーマンスと柔軟性の両方に対応するように設計されています。高レベルのロジックはPythonで記述されており、アクセシビリティと拡張性を高めている一方で、重い計算処理は最適化されたルーチンにオフロードされます。この設計により、スパース代数技術を用いることで、すべての可能な相互作用のフル行列を保持するのではなく、原子間の関連する相互作用のみを保持して、大規模なシステムを扱うことができます。この効率性により、単一のワークステーションで、以前は巨大なスーパーコンピュータを必要とした規模である、数十万の軌道を持つシステムの計算を行うことが可能になりました。

コンパクトなハミルトニアンの構築と解析を自動化することで、PAOFLOW 3.0は、大規模で一貫した材料データベースの作成を長年阻んできた計算上の障壁を取り除きます。これは、同じ基礎となる第一原理データから導出される、電子、光学、輸送、およびトポロジカルな特性を計算するための統一されたフレームワークを提供します。この一貫性は、材料科学のルールを学習するために膨大な量の信頼できるデータを必要とする人工知能モデルの訓練にとって不可欠です。このソフトウェアはオープンソースであり、世界の科学コミュニティがこれを基盤として構築できるように設計されており、多様な研究ワークフローへの統合を保証する幅広い入力フォーマットをサポートしています。新しい量子材料の発見や、より効率的なエネルギー技術の設計へと向かう中で、PAOFLOWのようなツールは、理論的な予測を実用的な現実へと変えるための必要なスピードと精度を提供します。