4. 拉伸压缩模拟的原子模型构建

分子动力学模拟拉伸与压缩，第一步不是设置应变速率，而是把“被拉伸的物体”在原子尺度上真实地搭出来。初始构型决定了位错形核位置、自由体积分布、链段取向和界面结合状态，这些因素往往比势函数选择更直接地影响屈服强度、弹性模量和断裂应变的模拟结果。本节围绕四类典型结构——晶体、非晶、分子链和界面——说明建模对象、生成方法与初始构型的检查要点。

一、晶体结构：从晶格常数到取向控制

晶体建模的核心是“晶格+取向+缺陷”。以面心立方金属铜为例，先确定晶格常数 a≈0.3615 nm，再按目标尺寸换算重复单元数。若模拟盒为 10 nm×10 nm×10 nm，则每个方向约需 27～28 个晶胞。构建时通常使用 LAMMPS 的 lattice fcc 3.615 配合 region box block 0 27 0 27 0 27 和 create_box、create_atoms 完成。

拉伸方向必须与晶体取向绑定。沿 [100] 拉伸与沿 [111] 拉伸，滑移系开动难度不同，屈服强度可相差数倍。因此建模时要用 orient 参数或旋转矩阵把目标方向对齐到加载轴。例如：

lattice fcc 3.615 orient x 1 0 0 orient y 0 1 0 orient z 0 0 1
region box block 0 27 0 27 0 27 units lattice
create_box 1 box
create_atoms 1 box

若研究含缺陷晶体，可在 create_atoms 后删除特定区域原子形成空位，或插入额外半原子面构造刃位错。初始构型生成后，必须做能量最小化（minimize）和短时弛豫（NVT 或 NPT，约 20～50 ps），让原子回到受力平衡位置，否则初始重叠原子会产生巨大斥力，导致模拟“爆炸”。

二、非晶结构：熔融淬火与自由体积控制

非晶金属或非晶聚合物的建模思路是“先无序、后弛豫”。常用方法是熔融-淬火法：先把晶体加热到熔点以上（如 Cu 加热到 2000 K）保温 100 ps 使其熔化，再以 10¹¹～10¹² K/s 的速率快速冷却到目标温度（如 300 K），得到无长程有序的非晶态。冷却速率越慢，非晶越致密，自由体积越少，模拟屈服强度通常越高。

非晶初始构型的关键检查量是径向分布函数 g(r) 和密度。若 g(r) 第二峰分裂不明显，说明非晶化不充分；若密度明显低于实验值（如 Cu 非晶约 8.0 g/cm³），说明存在过多孔洞，需要延长弛豫或调整淬火速率。非晶拉伸时没有明显位错滑移，塑性变形以剪切转变区（STZ）和自由体积聚集为主，因此初始自由体积分布直接决定剪切带的萌生位置。

三、分子链：聚合物的链长、缠结与取向

聚合物拉伸压缩模拟不能只放一堆单体，必须构建具有合理链长和缠结的分子链。常用工具是 Materials Studio、LAMMPS 的 molecule 模板或 Python 脚本（如 Polymatic）。以聚乙烯为例，单链重复单元为 –CH₂–CH₂–，若目标分子量 10 万，则每条链约含 3500 个碳原子。模拟盒中通常放置 20～50 条链，总原子数控制在 5 万～20 万，兼顾精度与算力。

链的初始构型常用“自避随机行走”生成，再经过能量最小化和 NPT 弛豫使体系达到平衡密度。判断平衡的指标包括：回旋半径 Rg 趋于稳定、端到端距离分布不再漂移、体系密度波动小于 2%。拉伸时，链段沿加载方向取向，取向度可用 Herman 取向因子 S 描述，S 从 0 逐渐趋近 1，对应应力-应变曲线上的取向硬化段。

四、界面：晶界、异质界面与结合能

界面建模要同时处理“几何匹配”和“原子结合”。晶界常用重合位置点阵（CSL）模型，例如 Cu 的 Σ5(310) 对称倾斜晶界，先分别构建两个取向晶粒，再按 CSL 关系拼接，删除重叠原子并弛豫。异质界面（如金属/陶瓷、金属/聚合物）则需考虑晶格失配度 δ=(a₂−a₁)/a₁，失配度大于 5% 时通常需要引入失配位错或选择半共格界面。

界面结合能用 E_b=(E_total−E_A−E_B)/A 计算，其中 E_total 为界面体系总能量，E_A、E_B 为两侧独立体系能量，A 为界面面积。结合能为负且绝对值越大，界面越稳定。拉伸时，界面往往是裂纹萌生点，因此初始构型中界面附近的原子应额外弛豫，避免因拼接产生的残余应力干扰后续力学响应。

五、初始构型的验证与常见问题

模型搭好后，不能直接进入拉伸。建议按以下顺序验证：第一，检查原子是否重叠，最近邻距离不应小于平衡键长的 0.7 倍；第二，运行 10 ps NVT，观察温度、势能是否稳定；第三，计算密度并与实验值或文献值对比，偏差应控制在 3% 以内；第四，对晶体检查 g(r) 是否出现尖锐峰，对非晶检查第二峰是否分裂，对聚合物检查 Rg 是否收敛。

常见问题有三类。一是“初始应力过大”，多因拼接界面或链段重叠，解决方法是延长弛豫或分步升温退火。二是“非晶化不充分”，表现为 g(r) 与晶体相似，需降低淬火速率或提高熔融温度。三是“分子链穿过周期性边界后缠结异常”，需检查链的连续性，必要时用 fix shake 或调整盒子尺寸。只有初始构型通过验证，后续的应力-应变曲线才具有物理意义。