[论文解读] Revisiting the Common Neighbour Analysis and the Centrosymmetry Parameter

Peter Mahler Larsen|arXiv (Cornell University)|Mar 19, 2020

Graph theory and applicationsMathematics被引用 21

一句话总结

本文提出了分子动力学中两种广泛使用的结构分析方法的改进变体：区间共邻分析（i-CNA）和最小权重匹配中心对称参数（MWM-CSP）。i-CNA 通过在一系列阈值范围内搜索，提高了对原子无序情况下的鲁棒性，避免了误分类；MWM-CSP 则利用图匹配方法，确保了中心对称性计算的一致性、连续性和物理直观性——两种方法在不增加参数的前提下，提升了准确性和可靠性。

ABSTRACT

We review two standard methods for structural classification in simulations of crystalline phases, the Common Neighbour Analysis and the Centrosymmetry Parameter. We explore the definitions and implementations of each of their common variants, and investigate their respective failure modes and classification biases. Simple modifications to both methods are proposed, which improve their robustness, interpretability, and applicability. We denote these variants the Interval Common Neighbour Analysis, and the Minimum-Weight Matching Centrosymmetry Parameter.

研究动机与目标

• 解决传统 CNA 和 CSP 方法在原子无序和不同长度尺度下分类晶体结构时的局限性。
• 解决现有 CNA 和 CSP 实现中在多晶或无序体系下存在的分类偏差和不一致结果问题。
• 开发稳健的无参数替代方法，在保持可解释性的同时提高结构识别的准确性。
• 确保在小原子位移下中心对称性测量的连续性和一致性，这是可靠模拟分析的关键要求。

提出的方法

• 提出区间共邻分析（i-CNA），通过在一系列邻居阈值范围内评估结构特征，而非依赖单一固定阈值，从而提高对热或结构涨落的鲁棒性。
• 为每个原子使用第一或第二壳层邻居距离计算局部长度尺度，以定义原子特异的阈值，实现在多相体系中的尺度不变性。
• 用最小权重匹配（MWM）替代 CSP 中的贪心或启发式匹配，在原子完全图上进行匹配，以平方距离作为边权重，寻找相反邻居的最佳配对。
• 确保每个原子恰好有一个相反邻居，并最小化总权重和，从而生成在小扰动下连续、不可微但稳定的 CSP 函数。
• 采用混合策略：首先应用更快的贪心边选择（GES）方法；仅当 GES 无法生成有效匹配时，才调用完整的 MWM 算法，从而提高计算效率。
• 通过包括 Bain 变换、理想晶体结构和多晶 Ru 模拟在内的测试案例验证新方法，以评估分类偏差和一致性。

实验结果

研究问题

• RQ1 如何使共邻分析（CNA）在分子动力学模拟中对原子位移和不同局部长度尺度更具鲁棒性？
• RQ2 现有 CSP 实现失败的原因是什么？如何使中心对称性计算在小原子畸变下保持连续且具有物理意义？
• RQ3 能否开发一种无参数、鲁棒的 CNA 变体，避免因阈值敏感性导致的误分类？
• RQ4 所提出的最小权重匹配 CSP 方法是否能在不同对称性和无序程度下产生一致且直观的结果？
• RQ5 与传统 CNA 和 CSP 相比，新方法在分类准确性、偏差和计算效率方面表现如何？

主要发现

• 区间共邻分析（i-CNA）通过在一系列邻居阈值范围内搜索，而非依赖单一固定阈值，显著提高了在无序或受扰动的局部环境中结构识别的准确率。
• i-CNA 方法成功通过了 Bain 变换测试，表明其无分类偏差，而传统 CNA 变体则存在此类问题。
• 最小权重匹配中心对称参数（MWM-CSP）生成了连续、稳定且直观的 CSP 值，即使在小原子位移下也能正确反映中心对称性的程度。
• 基于贪心的 CSP 方法（GES 和 GVM）无法产生一致结果：GVM 在旋转下不连续，而 GES 将截然不同的结构分配相似的 CSP 值，导致中心对称与非中心对称环境的分离效果差。
• 在多晶 hcp Ru 模拟中，MWM-CSP 在 CSP 分布中产生单一宽峰，而 GVM 产生虚假的次峰，GES 则导致人为狭窄的峰，表明分辨率差且存在过拟合。
• 混合 MWM-CSP 方法在单个 CPU 线程上可实现约 30,000 次 CSP 计算/秒，且默认使用更快的 GES 方法，仅在必要时回退至 MWM，有效平衡了速度与准确性。