4. 晶界结构表征
在分子动力学中，晶界并不是一条“线”，而是一个具有有限厚度、原子排列偏离理想点阵的过渡区。要对它做定量表征，必须同时回答三个问题：两侧晶粒的取向关系是什么、晶界面在各自晶体坐标系中如何取向、以及哪些原子真正属于晶界。本节围绕重合位置点阵（CSL）、取向差与晶界面指数，以及基于原子堆积与配位数的晶界原子识别方法展开。
一、从取向差到 CSL：晶界几何的数学描述
晶界几何由五个宏观自由度描述：三个描述两晶粒相对旋转的取向差参数，两个描述晶界面法向的参数。取向差通常用旋转轴 [uvw] 与旋转角 θ 表示，记为 θ/[uvw]。例如面心立方金属中常见的 Σ5 晶界，其取向差为 36.87°/[100]，即两晶粒绕公共 [100] 轴相对旋转 36.87°。
当两套点阵绕某轴旋转特定角度后，会出现一部分阵点周期性重合，这些重合阵点构成的新点阵即重合位置点阵（CSL）。CSL 的密度用 Σ 值表征：Σ 表示 CSL 单胞体积与晶体原胞体积之比，也等于重合阵点占原阵点的比例倒数。Σ 越小，重合密度越高，晶界通常越有序、能量越低。Σ1 对应完美晶体，Σ3 对应孪晶界，Σ5、Σ9、Σ11 等为常见高角度晶界。
需要强调：CSL 只描述两晶粒点阵的几何重合关系，并不直接给出晶界面的位置和取向。同一个 Σ 值可以对应多个不同晶界面指数，其原子结构、能量和迁移行为可能差异显著。因此，CSL 必须与晶界面指数联合使用。
二、晶界面指数与五参数表征
晶界面指数用晶界面法向在两侧晶粒坐标系中的 Miller 指数表示。对立方晶系，常用 (hkl) 表示法向。一个完整的晶界描述可写为：
θ/[uvw] + (hkl)₁ // (hkl)₂
其中 (hkl)₁ 和 (hkl)₂ 分别是晶界面法向在晶粒 1 和晶粒 2 坐标系中的指数，二者通过取向差旋转矩阵关联。例如 Σ5 (310)/[001] 对称倾斜晶界，表示两晶粒绕 [001] 轴相对旋转 36.87°，晶界面法向在两侧分别为 (310) 和 (3-10)，属于对称倾斜晶界。若两侧指数不对称，则为非对称倾斜晶界。
在 MD 建模中，通常先构造两个取向不同的晶粒，再按指定晶界面法向切割并拼接，最后通过能量最小化或高温退火消除拼接产生的重叠原子和过近原子。此时需要检查：拼接后是否出现原子间距小于 0.7a₀（a₀ 为晶格常数）的异常原子对，若有则删除或移动，否则会在后续积分中产生巨大斥力导致体系崩溃。
三、利用原子堆积与配位数识别晶界原子
晶界原子识别的核心思想是：晶界区原子排列偏离理想点阵，其局域堆积密度和配位数与完美晶格不同。常用方法包括：
- 配位数法：以截断半径
r_c（通常取第一近邻距离的 1.2 倍）统计每个原子的近邻数。FCC 完美晶体配位数为 12，BCC 为 8，HCP 为 12。晶界原子常出现配位数 10、11、13、14 等异常值。 - 原子堆积密度法：以原子为中心构造 Voronoi 多面体，计算其体积。晶界区原子 Voronoi 体积偏离完美晶格值超过一定阈值（如 5%）时判为晶界原子。
- 公共近邻分析（CNA）：通过近邻原子间的成键拓扑判断局域结构，FCC、BCC、HCP 各有特征 CNA 签名，无法匹配的原子归为“其他”，即晶界或缺陷原子。
- 势能判据：晶界原子势能通常高于完美晶格原子，可结合势能分布直方图辅助判断。
实际分析中，建议先计算完美晶体的配位数和势能基准，再对含晶界体系做同样统计，以基准值为参考设定阈值。例如某 FCC 金属 a₀ = 3.615 Å，第一近邻距离约 2.556 Å，取 r_c = 3.0 Å，完美晶体配位数为 12；若某原子配位数 ≤ 11 或 ≥ 13，则标记为晶界候选原子。再结合 Voronoi 体积偏差 > 5% 做二次筛选，可有效排除热振动引起的误判。
四、实操建议与结果判读
在 OVITO 或自编 Python 脚本中，可按以下顺序操作：先读取 dump 文件中的原子坐标和类型；用 pbc 处理周期性边界，避免界面处近邻统计错误；对每个原子计算截断半径内的近邻列表；统计配位数并输出直方图。若配位数分布出现明显双峰——主峰在 12、次峰在 10~11 或 13~14——说明晶界原子已被有效区分。
进一步，可将识别出的晶界原子投影到晶界面法向，得到晶界厚度。对 Σ5 (310)/[001] 对称倾斜晶界，典型厚度约 0.5~1.0 nm，对应 2~3 个原子层。若厚度异常大（如超过 2 nm），可能说明建模时晶界面取向偏离设定值，或退火温度过高导致晶界预熔化。
r_c 取值不当会直接改变配位数。FCC 金属中若 r_c 取到第二近邻距离（约 a₀），配位数会变成 18，失去区分度。务必先计算径向分布函数第一峰位置，再取 r_c 为第一峰右侧谷底位置。
综上，CSL、取向差和晶界面指数共同构成晶界的五参数几何描述，而配位数、Voronoi 体积、CNA 和势能则是从原子尺度识别晶界原子的有效判据。二者结合，才能把“晶界在哪里、是什么类型、由哪些原子组成”这三个问题同时回答清楚，为后续晶界能计算、迁移行为模拟和力学响应分析提供可靠的结构基础。