2. 原子尺度建模基础

分子动力学（MD）模拟晶界与位错，第一步不是写势函数，而是把“原子放在哪里”这件事说清楚。原子尺度建模的本质，是用一组带坐标的粒子去逼近一块真实晶体。要做到这一点，必须同时回答四个问题：晶体用什么点阵描述、方向与面如何标记、原子坐标如何生成、边界如何处理。本节围绕这四个问题，给出单晶与双晶模型的工程化构建方法。

一、晶体结构：从点阵到原子基元

晶体结构 = 空间点阵 + 基元。点阵是抽象的周期性格点，基元是附着在每个格点上的原子组合。面心立方（FCC）常用三个晶格常数描述：晶格常数 a，以及基元原子位于 (0,0,0)、(0,½,½)、(½,0,½)、(½,½,0)（分数坐标）。体心立方（BCC）基元为 (0,0,0) 与 (½,½,½)。密排六方（HCP）需额外给出 c/a 比，理想值为 √(8/3)≈1.633。

工程建模中，金属材料最常用的是 FCC（Cu、Al、Ni、Ag、Au）与 BCC（Fe、W、Mo、Ta）。以 Cu 为例，a = 3.615 Å，最近邻距离 d = a/√2 ≈ 2.556 Å。这个数值直接决定势函数截断半径的选取：截断半径通常取 2.5d～3.5d，即约 6.4～8.9 Å。若截断半径小于 2d，第二近邻被截断，弹性常数会明显失真。

二、晶向与晶面：Miller 指数的工程含义

晶向 [uvw] 表示沿点阵矢量的方向，晶面 (hkl) 表示一组平行等间距的原子面。两者关系由晶带定律约束：若 [uvw] 位于 (hkl) 面内，则 hu + kv + lw = 0。面间距 d(hkl) 对 FCC 为 a/√(h²+k²+l²)（h、k、l 全奇或全偶时），对 BCC 需满足 h+k+l 为偶数。

这些量在建模中不是装饰。晶界能、位错芯结构、滑移系选择都依赖晶面与晶向。例如 FCC 的滑移系为 {111}⟨110⟩，共 12 个；构建刃位错时，位错线沿 [110]，滑移面为 (111)，Burgers 矢量 b = a/2[110]，其模长 |b| = a/√2 ≈ 2.556 Å。若把位错线方向写错，得到的将不是刃位错而是螺位错或混合位错。

三、原子坐标的生成：从分数坐标到笛卡尔坐标

建模时先在分数坐标下生成规则点阵，再乘以晶格矩阵转换为笛卡尔坐标。对正交晶系，转换就是逐轴乘 a、b、c；对 HCP 或三斜晶系，需用完整的晶格矩阵。以 FCC 为例，在 nx×ny×nz 个单胞上循环，对每个单胞叠加四个基元原子，即可得到完整单晶。

# 以 Cu FCC 为例，生成 10×10×10 单晶（分数坐标 → 笛卡尔坐标）
a = 3.615            # Å
basis = [(0,0,0), (0,0.5,0.5), (0.5,0,0.5), (0.5,0.5,0)]
atoms = []
for i in range(10):
    for j in range(10):
        for k in range(10):
            for bx, by, bz in basis:
                x = (i + bx) * a
                y = (j + by) * a
                z = (k + bz) * a
                atoms.append((x, y, z))
# 共 10*10*10*4 = 4000 个原子

生成后必须做两件事：一是检查最近邻距离是否等于理论值（Cu 为 2.556 Å），二是统计原子数是否与理论值一致。若最近邻距离偏大或偏小，说明基元或晶格常数写错；若原子数不符，说明循环边界或基元重复出错。这两步是建模阶段最廉价、最有效的自检。

四、周期性边界条件：用有限原子模拟无限晶体

MD 模拟的原子数通常只有 10³～10⁶，而真实晶体宏观上近似无限。周期性边界条件（PBC）的做法是：把模拟盒在三个方向上首尾相接，任一原子从一侧穿出，就从对侧以相同速度穿入。计算原子间作用力时采用最小镜像约定：两原子距离取其在所有周期镜像中的最短距离。

PBC 带来两个直接约束。第一，模拟盒尺寸必须大于两倍截断半径，否则同一原子会与自身的镜像发生非物理相互作用。若截断半径为 8 Å，盒长至少 16 Å，工程上常取 20 Å 以上。第二，PBC 下不能直接模拟单个位错或自由表面，因为位错应力场会与自身镜像叠加。此时需改用固定边界、柔性边界或把位错放在足够大的盒中并扣除镜像应力。

五、双晶模型：从单晶到晶界的工程构建

双晶模型是研究晶界的基本对象。构建思路是：先分别生成两个取向不同的单晶，再按选定晶界面拼接，最后删除重叠原子并做能量最小化。以 Cu 的 Σ5(310)[001] 对称倾斜晶界为例，取向差 θ = 36.87°，两晶粒绕 [001] 轴相对旋转 θ/2 与 −θ/2。

拼接时需注意三点。第一，两晶粒在界面处的原子不能过近，否则会产生巨大排斥力，通常设定最小间距阈值（如 1.5 Å），删除过近原子。第二，界面两侧需留出足够真空或缓冲层，避免 PBC 下界面与自身镜像耦合。第三，拼接后必须做能量最小化（如共轭梯度法），让原子弛豫到局部能量极小，否则初始构型中的高应力会导致模拟初期“爆炸”。

判断双晶模型是否合理，可观察三个量：界面能是否收敛到文献值附近（Cu Σ5(310) 约 0.7～0.9 J/m²）；界面处原子是否形成规则的重合位置点阵（CSL）结构；弛豫后界面是否保持平直而非严重弯曲。若界面能远高于文献值，通常是重叠原子未删除或取向角写错。

六、建模检查清单

把上述内容落到可执行的检查项，单晶与双晶模型在提交计算前应逐项确认：晶格常数与基元是否与目标材料一致；最近邻距离是否等于理论值；原子总数是否等于理论值；模拟盒是否大于两倍截断半径；双晶取向角与晶界面指数是否匹配；重叠原子是否删除；是否完成能量最小化；界面能是否收敛到合理范围。这八项中任何一项出错，后续的位错迁移、晶界迁移或力学响应计算都会失去意义。

原子尺度建模的工程性正体现在这里：它不追求复杂的算法，而要求每一个坐标、每一个方向、每一个边界条件都经得起核对。把单晶与双晶模型建对，是后续所有晶界与位错模拟的前提。