[论文解读] 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结构。

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）上验证该方法。

提出的方法

• 构建一个表示晶体基元中原子连接性的多向周期图，节点代表原子，边代表化学键。
• 通过图中闭合路径定义环基，以边向量的整数组合表示。
• 制定编码图中环的拓扑一致性条件的代数关系（alpha向量）。
• 应用线性代数，通过Gram矩阵的矩阵求逆与Cholesky分解计算图的标准实现。
• 通过从参考节点沿最短路径的边位移向量累加，生成原子坐标。
• 使用pymatgen构建最终晶体结构，通过晶格缩放至目标体积，并进行原子位置合并以保证化学一致性。

Figure S1: Basis of closed paths of FCC structure

实验结果

研究问题

• RQ1 如何系统性地将晶体中的多面体基元编码为具有拓扑约束的周期图？
• RQ2 何种代数框架可表示周期图的环空间，以确保结构一致性？
• RQ3 双重周期图表示能否生成常见密堆积与体心晶格的物理上合理的晶体结构？
• RQ4 边向量与环基向量如何关联到晶体结构的几何实现？
• RQ5 Gram矩阵与边位移向量在从图拓扑重建原子坐标的过程中起什么作用？

主要发现

• 该方法成功使用统一的图论框架，为FCC、HCP与BCC晶格生成了双重周期结构。
• 图的环空间中的闭合路径基以整数向量（alpha）编码，确保了拓扑一致性。
• 通过Gram矩阵的矩阵求逆与Cholesky分解，计算出图的标准实现。
• 原子坐标通过从参考原点沿路径累加边位移向量重建。
• 通过pymatgen的结构生成与晶格缩放对结果结构进行验证，并通过位置合并确保化学一致性。
• 该方法实现了无需预先知晓原子位置的自动化、拓扑驱动的晶体结构生成。

Figure S2: Basis of closed paths of HCP structure