四十年硅经典原子间势的热学、缺陷与力学性能基准评测：面向摩擦学、磨损与切削

《Computational Materials Science》：Four decades of classical interatomic potentials for silicon: Benchmarking thermal, defect and mechanical performance for tribology, wear and machining

　　

七种经典硅原子间势（classical interatomic potentials）在热响应与缺陷响应方面接受基准评测。不同硅势在切削过程中预测出不同的应力适应（stress-accommodation）路径。Erhart–Albe与Tersoff倾向于强烈局域化的亚表面（subsurface）变形。Stillinger–Weber促进切削过程中的瞬态位错（transient dislocation）活动。电子显微术（electron microscopy）揭示非晶化（amorphisation）、层错（stacking faults）与部分位错（partial dislocations）。

该研究面向硅在半导体、微机电系统（micro-electromechanical systems, MEMS）、红外光学和精密工程中的纳米级表面与低亚表面损伤需求，指出单点金刚石车削、精密磨削、晶圆划片和抛光中，硅去除受脆性断裂与高度受限塑性变形竞争控制；分子动力学（molecular dynamics, MD）可解析纳米—皮秒尺度变形，但以往多固定一种势或只比较单一性质，缺少把势公式—基础性质—缺陷机制—应力适应—加工亚表面结构联系起来的机制基准，因此难以解释同一接触载荷下不同势何以预测根本不同的亚表面变形。研究人员构建分级基准：先筛选七种硅势的热响应、相内聚能、广义层错能（generalized stacking-fault energy, GSFE）和拉伸行为，再将五种代表性势推进至位错核、位错线/核能、应力驱动迁移率、层错扩展和纳米切削，并以磨削单晶硅的透射电子显微镜（transmission electron microscopy, TEM）观察作为制造尺度机制参照。结论是没有任何单一势在热学、弹性、缺陷和加工性质上全面优越，切削响应由势能面编码的应力适应路径而非单一峰值应力或单一能量描述符决定，势选择应依据所关心的机制和观测量。
研究方法概括不超过250字：采用LAMMPS与OVITO，七势为SW、Tersoff、Erhart–Albe ABOP、Kumagai（Mod）、Screened ABOP、Starikov ADP和Pun–Mishin modified Tersoff；热评估用64,000原子NPT加热300—4000 K、250 ps；相内聚能约8,000原子；GSFE用76,800原子；拉伸用直径20.68 nm、长48.98 nm的?100?纳米线；五种势用于位错偶极、迁移率和层错扩展；切削含约663,000个Si原子与约496,000个C原子刀具，?25°前角、10°后角、250 m s^?1、50 nm；实验样本为(011)取向、直径200 mm单晶硅片，减薄至50 μm，在VG401 MKII精密磨床上以去离子水冷却进行条件1火花光磨和条件2无火花光磨，Ar⁺离子抛光后以Titan Themis G2 300在300 kV下观察。
结果：3.1 NPT加热热评估：通过连续加热轨迹比较密度、热膨胀系数（coefficient of thermal expansion, CTE）和表观热失稳温度（T_inst），得到T_inst约2309—3427 K，均显著高于实验平衡熔点（T_m）1685 K，说明直接快速加热值应视为超热/热失稳指标而非T_m。3.2 相稳定性：通过零压弛豫和内聚能计算，所有势均保持SiI最稳定并再现SiI>bct5>Si-II的零压能量顺序；SiI内聚能4.336—4.6299 eV/atom，SW低估约6.4%，高压相差异更大。3.3 GSFE剖面：通过{111}滑移集?112?和shuffle集?110?剪切能景计算，滑移集γ_usf为1.918—2.948 J m^?2，Tersoff最低、SW最高；shuffle集Tersoff最高2.281 J m^?2、SW最低0.841 J m^?2；滑移集与shuffle集排序相反，γ_usf一致不保证完整剖面一致。3.4 位错核结构：通过{110}、{111} shuffle和{111} glide位错偶极弛豫及应力迁移前后比较，ABOP与Tersoff初始核较局域且glide集可见Shockley部分位错解离，SW迁移后核显著无序和扩展，Kumagai重构明显，Screened ABOP视觉紧凑但拟合核能极低。3.5 位错线能量与拟合弹性参数：通过位错偶极间距依赖拟合，拟合剪切模量G约51.8—76.0 GPa，实验约68 GPa，Kumagai 65.6 GPa最接近；E_c差异大，Screened ABOP的{110}核能0.1109 eV ?^?1，较ABOP的0.6398 eV ?^?1低约5.8倍。3.6 位错迁移率：通过300 K、3 GPa剪切下位移轨迹，SW在三种构型均极高迁移，Tersoff呈强晶体学选择性（{110}/glide受限而shuffle可持续），ABOP受限，Kumagai间歇脱钉，Screened ABOP快速初始位移后停滞。3.7 层错扩展：通过600 K、3 GPa下预存glide集层错带扩展，Kumagai较SW产生更快更宽的部分位错分离，说明动态层错扩展不能仅由γ_usf推断。3.8 硅纳米线拉伸性能：通过?100?纳米线300 K、5×10⁸ s^?1单轴拉伸，ABOP给出最高极限拉伸强度17.39 GPa，Tersoff 15.59 GPa且断裂应变0.314最大，SW与ADP较低；杨氏模量65.06—99.93 GPa，ABOP最高，Tersoff、Pun–Mishin、ADP等较低。3.9 硅纳米切削：通过IDS/DXA分析、切削力分布和区域平均应力，ABOP与Tersoff形成紧凑转变/无序区、无持续DXA可检位错且力高而窄；SW形成较弥散区并唯一产生显著瞬态a/2?110?位错爆发和Shockley部分位错，力低而宽；Kumagai无持续位错但转变区渐宽不规则，处中间；Screened ABOP形成最扩展、可深入数纳米的转变/无序区并仅晚期出现约35 ?短位错，力中等而不规则；von Mises应力约6—10 GPa、Kumagai最低，Tresca约3.4—5.7 GPa、Tersoff较高，ABOP与Screened ABOP压应力幅值最大。4 磨削诱发亚表面变形的实验TEM基准：通过条件1和条件2磨削截面BF-TEM、DF-TEM和HRTEM观察，磨削硅近表面出现几十纳米非晶α-Si层及非晶、畸变晶体、{111}层错和部分/完美位错共存结构，Burgers矢量与1/6?112?、1/3?111?和1/2?110?型缺陷一致，表明加工亚表面变形由非晶化、层错和位错介导塑性竞争共存而非单一机制主导。
讨论总结：该研究强调差异最根本体现在严重局域应力的适应路径，内聚能、杨氏模量、γ_usf、峰值von Mises应力或拉伸强度等单一描述符均不能独立解释切削响应；ABOP与Tersoff代表强局域适应，SW代表叠加结构无序的瞬态位错介导塑性且其高迁移率提示定量位错动力学需谨慎，Screened ABOP代表伴随异常低拟合核能的广泛无序介导适应，Kumagai为兼具较真实弹性和较宽结构无序的中间响应。TEM支持非晶化、{111}层错和部分位错共存的混合机制，但磨削为重复磨粒接触累积损伤，而MD为理想单程约200 ps切削，比较应保持机制性而非一一验证。局限包括单一晶体学/刀具/速度/加载历史、高应变率、经典势未显式表征电子与掺杂效应、直接加热T_inst非T_m、位错能拟合采用固定泊松比0.3和r₀=2 ?；未来应将同一基准用于机器学习势，但更高表示能力仍需机制特定验证。
研究结论翻译：本研究证明原子间势选择根本影响硅加工原子模拟所预测的变形机制。基准中表观热失稳温度跨2309—3427 K，第一{111}滑移集Shockley部分位移的γ_usf为1.918—2.948 J m^?2，杨氏模量为65.06—99.93 GPa，拟合{110}位错核能为0.1109—0.6398 eV ?^?1，表明与任一单一体性质一致不足以确立势对高度受限接触变形的可迁移性。核心发现是加工响应不由峰值应力单独决定，而由势能面编码的应力适应路径决定。ABOP（Erhart–Albe）和经典Tersoff倾向强局域适应，形成较紧凑初始缺陷结构并受限持续位错传播；Stillinger–Weber更易激活瞬态位错介导塑性；Screened ABOP促进广泛无序介导适应并伴异常低位错核能；Kumagai势呈中间响应，兼具较真实弹性行为和动态加载下增加的结构无序。值得注意的是，Stillinger–Weber预测最高滑移集不稳定层错能2.948 J m^?2却表现极高应力驱动位错迁移率，说明静态不稳定层错垒不足以预测位错动力学，位错核结构、局部键重排和完整缺陷能景也必须考虑。实验TEM观察为这些原子趋势提供独立制造尺度参照：磨削硅表现共存非晶化、{111}层错和部分位错活动，表明加工诱发亚表面变形来自竞争和共存机制而非单一主导过程；比较仍应定性，因为磨削通过重复磨粒接触累积损伤，而MD模拟为理想单程纳米切削。总体结果支持机制特定选势策略：无单一经验势在热、弹性、缺陷和加工性质上全面优越；本基准内ABOP和Tersoff倾向强局域亚表面变形，Stillinger–Weber有助于探索位错介导塑性，Screened ABOP突出异常低拟合位错核能与无序介导应力适应的后果，Kumagai提供兼具较真实弹性和更广泛结构无序的中间响应。