3. 原子尺度建模基础与工程案例
原子尺度仿真是连接量子力学与连续介质力学的桥梁，其核心价值在于用可计算的原子间相互作用，预测材料在纳米尺度下的结构演化与性能响应。本节围绕三类主流方法——分子动力学（MD）、蒙特卡洛（MC）与密度泛函理论（DFT）——梳理其基本原理、势函数选择、边界条件设置，并通过典型材料案例说明工程应用中的关键决策。
一、三种方法的基本原理与适用边界
分子动力学基于牛顿第二定律，对每个原子求解运动方程。给定初始位置与速度，通过势函数计算原子受力，再用Verlet或Velocity-Verlet算法逐步积分，时间步长通常取1 fs量级（需小于最高振动频率周期的1/10）。MD天然输出轨迹，可统计温度、压力、径向分布函数、均方位移等，适合研究熔化、扩散、位错运动、辐照损伤等动力学过程。
蒙特卡洛不直接求解运动方程，而是通过随机抽样接受或拒绝原子构型，使系统按玻尔兹曼分布演化。Metropolis准则中，若新构型能量降低则接受，若升高则以exp(-ΔE/kT)概率接受。MC擅长处理平衡态统计问题，如合金有序-无序转变、吸附构型搜索、空位团簇演化，但不提供真实时间尺度。
密度泛函理论从电子密度出发求解Kohn-Sham方程，无需经验势函数，可给出键能、能带、磁矩、反应路径等电子结构信息。DFT计算量随原子数呈O(N³)增长，常规体系限于数百原子，适合验证势函数、计算缺陷形成能、解析界面成键机制。
二、势函数选择：从LJ到MEAM
MD结果的可靠性首先取决于势函数。工程中常见的选择路径如下：
- Lennard-Jones势：仅含排斥与色散项，参数少，适合惰性气体、简单液体及教学验证，但无法描述金属键的方向性。
- 嵌入原子法（EAM）：将能量分为嵌入能和电子密度项，能较好描述fcc/bcc金属的弹性常数、空位形成能和表面能，是Cu、Al、Fe等金属体系的主流选择。
- MEAM/ADP：在EAM基础上引入角度依赖，改善对共价键和层错能的描述，适用于Si、金属间化合物及合金。
- 反应势（ReaxFF）：允许键断裂与形成，适合碳氢体系燃烧、氧化及界面反应，但计算成本显著高于EAM。
- 机器学习势：以DFT数据训练，精度接近DFT而速度接近经验势，适合复杂合金与界面，但需严格验证外推能力。
选择时建议先明确目标属性：若关注弹性与热膨胀，EAM通常足够；若涉及化学键重组，需ReaxFF或机器学习势；若仅做方法验证，LJ即可。势函数确定后，应通过晶格常数、内聚能、空位形成能等基准量进行校验，偏差超过5%时需谨慎使用。
三、边界条件：周期、自由与固定
边界条件直接决定模拟能否代表真实材料。周期边界条件（PBC）通过最小镜像约定消除表面效应，适合模拟块体材料；自由边界条件用于纳米颗粒、薄膜和表面体系；固定边界条件则常用于压痕、拉伸等力学加载，将底部若干层原子冻结以施加约束。
工程中常采用混合边界：例如纳米压痕模型中，底部固定、侧面周期、顶部自由，压头以恒定速度下压。此时需注意PBC方向上的盒子尺寸应大于势函数截断半径的两倍，避免原子与自身镜像相互作用。对于DFT计算，真空层厚度一般取10～15 Å，k点网格需做收敛测试，总能量变化小于1 meV/atom方可接受。
四、典型工程案例：Cu纳米压痕的MD模拟
以单晶Cu纳米压痕为例，说明完整建模流程。使用EAM势（如Cu01.eam.alloy），创建约20×20×15 nm³的fcc超胞，约含50万原子。设置x、y方向为周期边界，z方向底部5层固定，顶部自由。系统先在300 K、NVT系综下弛豫100 ps，使能量收敛。
随后以半径5 nm的虚拟球形压头沿[001]方向下压，下压速度10 m/s，最大深度3 nm。模拟中每1 ps记录一次载荷-位移曲线。典型结果表现为：初始阶段载荷随深度近似线性上升，对应弹性变形；当深度超过约1 nm后，载荷出现锯齿状突降，对应位错形核与滑移。通过共同邻域分析（CNA）可识别层错和孪晶结构，位错密度随压入深度增加而升高。
五、方法选择与参数设置速查
| 方法 | 时间尺度 | 空间尺度 | 关键参数 | 典型输出 |
|---|---|---|---|---|
| MD | ps～μs | nm～百nm | 势函数、Δt、系综、边界 | 轨迹、力学响应、扩散系数 |
| MC | 无真实时间 | nm～μm | 能量模型、温度、步数 | 平衡构型、有序度、占位 |
| DFT | ps以下 | <2 nm | 泛函、k点、截断能、真空层 | 能带、形成能、反应路径 |
实际工程中，三类方法常形成多尺度链条：DFT计算缺陷形成能与迁移势垒，用于拟合或验证EAM势；MD模拟位错与界面演化，提取本构关系；MC则用于处理长时间扩散与相变问题。理解各自的能力边界与参数敏感性，是获得可信原子尺度结论的前提。