[论文解读] On Contact Numbers of Finite Lattice Sphere Packings and the Maximal Coordination of Monatomic Crystals

Samuel Reid|arXiv (Cornell University)|Feb 6, 2016

Computational Geometry and Mesh Generation参考文献 7被引用 3

一句话总结

本文提出了一种算法框架，用于计算三维晶格球体堆积中最大接触数，建立了该几何问题与单原子晶体中最大配位数之间的直接联系。主要贡献在于提出了一种系统性方法，以确定相同原子构成的晶体内可能达到的最高原子接触数（化学键数），并通过基于晶格的球体堆积构型，对最多 n 个原子的情况提供了理论边界和计算验证。

ABSTRACT

We algorithmically characterize the maximal contact number problem for finite congruent lattice sphere packings in $\mathbb{R}^d$ and show that in $\mathbb{R}^3$ this problem is equivalent to determining the maximal coordination of a monatomic crystal.

研究动机与目标

• 解决有限等径球体在三维晶格上堆积时最大接触对数（接触数）的开放性问题。
• 在球体堆积的接触数问题与单原子晶体中的最大配位数之间建立正式联系。
• 开发一种计算上可行的算法，用于计算任意给定原子数 n 在晶格基晶体结构中的最大接触数。
• 利用离散几何方法，为单原子晶体中的化学键数提供理论边界，改进现有渐近估计。

提出的方法

• 本文引入一种基于晶格的模型，其中原子被表示为位于三维整数晶格点上的等径球体，且原子间距离至少为球体半径的两倍。
• 定义最大接触数 $ C_3(\tilde{P}_\text{Λ}, n) $ 为在晶格堆积 $ \tilde{P}_\text{Λ} $ 的任意 n 个球体子集中，接触对的最大数量。
• 该算法将候选构型限制在晶格中的有界立方区域内，具体为 $ 0 \leq \lambda_i \leq \lceil n/3 \rceil $，以确保计算可行性。
• 系统枚举候选晶格点中所有大小为 n 的子集，并计算接触图的边数以识别最大值。
• 该方法依赖于来自离散几何定理的理论边界 $ C_3(\tilde{P}_\Lambda, n) < 6n - 3.665\dots n^{2/3} $，用于验证结果。
• 该框架基于单原子球体堆积对应关系，通过接触图将化学键转化为球体接触。

实验结果

研究问题

• RQ1 在三维晶格上排列的 n 个等径球体的有限堆积中，可能实现的最大接触数（化学键数）是多少？
• RQ2 晶格球体堆积中的最大接触数与单原子晶体中的最大配位数之间有何关系？
• RQ3 能否通过算法方法高效计算任意 n 在三维晶格球体堆积中的最大接触数？
• RQ4 此类堆积中限制接触数的几何与组合约束是什么？与一般球体堆积相比有何不同？
• RQ5 八面体构造及其部分变体在生成最优接触构型方面可达到何种程度的适用性？

主要发现

• 三维有限晶格球体堆积的最大接触数上界为 $ 6n - 3.665\dots n^{2/3} $，该边界比一般球体堆积估计值 $ 6n - 0.926n^{2/3} $ 更为紧密。
• 该算法通过在每个方向上大小为 $ \lceil n/3 \rceil $ 的有界晶格区域内探索所有 n 元子集，成功计算出 $ C_3(\tilde{P}_\Lambda, n) $。
• 对于特定的 n 值，如 $ n = \frac{2k^3 + k}{3} $，八面体构造可达到理论最大接触数，证实其在这些情况下的最优性。
• 该方法证实，单原子晶体中的最大配位数恰好对应于有限晶格球体堆积中的最大接触数。
• 该框架为小规模 n 提供了一种系统化方法，用于生成和验证最优接触构型，将先前构造的范围扩展至 n=20。
• 通过将搜索空间限制在每个晶格方向上大小为 $ \lceil n/3 \rceil $ 的立方区域，该算法确保了计算效率，避免了在整个晶格上进行穷举搜索。