[論文レビュー] Crystal Structure Generation Based on Polyhedra using Dual Periodic Graphs

Tomoyasu Yokoyama, K. Ichikawa+1 more|arXiv (Cornell University)|Dec 14, 2023

Machine Learning in Materials Science参考文献 18被引用数 7

ひとこと要約

本稿では、多面体モチーフから導出される二重周期的グラフを用いて結晶構造生成のための新規手法を提案する。周期的グラフにおける閉路ベースを活用し、原子の結合性をモデル化する。この手法により、辺空間およびサイクル空間における代数的関係を用いて位相的制約を符号化することで、FCC、HCP、BCC構造の安定な結晶構造の体系的列挙が可能となる。グラフ理論および線形代数を用いた計算的実装により、FCC、HCP、BCC構造への応用に成功した。

ABSTRACT

Crystal structure design is important for the discovery of new highly functional materials because crystal structure strongly influences material properties. Crystal structures are composed of space-filling polyhedra, which affect material properties such as ionic conductivity and dielectric constant. However, most conventional methods of crystal structure prediction use random structure generation methods that do not take space-filling polyhedra into account, contributing to the inefficiency of materials development. In this work, we propose a crystal structure generation method based on discrete geometric analysis of polyhedra information. In our method, the shape and connectivity of a space-filling polyhedron are represented as a dual periodic graph, and the crystal structure is generated by the standard realization of this graph. We demonstrate that this method can correctly generate face-centered cubic, hexagonal close-packed, and body-centered cubic structures from dual periodic graphs. This work is a first step toward generating undiscovered crystal structures based on the target polyhedra, leading to major advances in materials design in areas including electronics and energy storage.

研究の動機と目的

• 多面体モチーフを用いた周期的グラフ理論に基づく、結晶構造生成の体系的手法の開発。
• 多方向周期的グラフにおける閉路ベースを用いて、結晶内の原子結合性をモデル化すること。
• サイクル空間および辺空間の代数的構造を用いて、結晶構造の位相的制約を形式化すること。
• グラフ表現に基づく、二重周期的構造の自動生成を可能とすること。
• 代表的な密な格子および体心格子構造（FCC、HCP、BCC）に対して、本手法の有効性を実証すること。

提案手法

• 原子結合性を表す多方向周期的グラフを構築し、ノードを原子、エッジを化学結合とする。
• グラフ内の閉路を用いてサイクル基底を定義し、エッジベクトルの整数線形結合として表現する。
• グラフのサイクルにおける位相的一致性条件を符号化する代数的関係（αベクトル）を定式化する。
• 線形代数を用いて、αベクトルから導かれるグラム行列に対する行列逆行列計算およびコレスキー分解により、標準的実現を計算する。
• 基準ノードからの最短パスに沿ってエッジ変位ベクトルを経路的に和算することで、原子座標を生成する。
• pymatgenを用いて最終的な結晶構造を構築し、所望の体積に格子スケーリングを行い、サイト統合により化学的一致性を確保する。

Figure S1: Basis of closed paths of FCC structure

実験結果

リサーチクエスチョン

• RQ1 結晶内の多面体モチーフは、位相的制約を伴う周期的グラフとしてどのように体系的に符号化可能か？
• RQ2 周期的グラフのサイクル空間を表現する代数的枠組みは、構造的一致性を保証するためにどのように構築可能か？
• RQ3 二重周期的グラフ表現は、代表的な密な格子および体心格子系の物理的に妥当な結晶構造を生成可能か？
• RQ4 エッジおよびサイクル基底ベクトルは、結晶構造の幾何的実現とどのように関係するか？
• RQ5 グラム行列およびエッジ変位ベクトルは、グラフ位相から原子座標を再構築する際に果たす役割は何か？

主な発見

• FCC、HCP、BCC格子に対して、統一的なグラフ理論的フレームワークを用いて二重周期的構造を成功裏に生成した。
• グラフのサイクル空間における閉路ベースは、整数ベクトル（α）として符号化され、位相的一致性を強制する。
• 標準的実現は、αベクトルから導かれるグラム行列に対する行列逆行列計算およびコレスキー分解により計算された。
• 原子座標は、基準原点からの経路に沿ってエッジ変位ベクトルを和算することで再構築された。
• pymatgenを用いた構造生成および格子スケーリングによる検証が行われ、サイト統合により化学的一致性が確保された。
• 本手法により、原子位置の事前知識が不要な、位相駆動型の自動的結晶構造生成が可能となった。

Figure S2: Basis of closed paths of HCP structure