势函数理论计算：分子力场参数化、拟合策略与交叉验证方法

σ=2.5 Å 的势阱与实验晶格常数的冲突

势函数参数化是连接第一性原理计算与大规模分子动力学模拟的桥梁，也是最容易被低估难度的环节。本项目在开发Cu-Al合金的EAM（嵌入原子法）势函数时，初始基于Lennard-Jones势拟合Cu-Cu相互作用，取σ=2.5 Å、ε=0.4 eV，结果纯Cu的晶格常数预测为3.52 Å，与实验3.615 Å偏差达2.6%。更严重的是，Cu-Al合金的有序化温度（B2结构的形成温度）被高估了200 K。问题的根源在于：Lennard-Jones势是两体势，无法描述金属的嵌入电子密度效应，而EAM势的嵌入能项（F(ρ̄)）对短程排斥的敏感性远超Lennard-Jones的r⁻¹²项。项目组放弃了两体势路径，转向基于DFT数据（VASP计算了Cu、Al、Cu₃Al的12个构型的能量-体积曲线）的EAM势拟合，采用Mishin等人的 fitting procedure（Acta Mater. 49, 2001），最终晶格常数偏差缩至0.2%，有序化温度偏差缩至30 K。这个经历确立了势函数开发的基本原则：两体势只适用于范德华体系，金属/合金必须多体势，且拟合数据必须覆盖DFT+实验双重空间。

困境累积：DFT数据与实验数据的权重博弈

势函数拟合的目标函数通常包含多个子项：E_coh（ cohesive energy）、a₀（晶格常数）、C₁₁/C₁₂/C₄₄（弹性常数）、E_vac（空位形成能）、E_sf（层错能）、以及高对称路径的声子频率。每个子项的权重决定了势函数的”性格”——权重偏向E_coh，势函数预测相图精确但弹性常数偏差大；权重偏向Cij，力学响应精确但缺陷能可能失真。本项目在Cu的EAM势拟合中测试了5种权重方案：

方案E通过微调权重分配，在 cohesive properties 和 defect properties 之间取得了平衡。关键的折中在于：Cij的权重不能太高（金属EAM势对剪切模量的描述天然弱于体积模量），E_vac的权重不能太低（空位迁移是扩散和蠕变的核心）。拟合采用Levenberg-Marquardt非线性最小二乘法，迭代2000步，RMS误差从初始0.45 eV收敛至0.08 eV。

关键抉择：拟合范围与可移植性测试

势函数的可移植性决定了它能否从”训练体系”外推到”应用体系”。本项目在Cu-Al合金势中，训练数据覆盖了5种结构：FCC Cu、BCC Al、B2 CuAl、L1₂ Cu₃Al、以及一个非晶态构型。但应用时遇到了Cu-Al界面（半共格界面，错配度4.3%）的模拟，界面能预测为0.42 J/m²，而实验HRTEM测量值为0.38 J/m²，偏差10%。进一步分析发现，训练数据中缺少错配位错核心的高能量构型，而EAM势在描述位错芯的压缩-剪切耦合响应时出现了非物理的软化。解决策略是：在训练集中加入3个含错配位错的超胞构型（DFT计算，每个约200原子，耗时24小时），重新拟合后界面能偏差缩至3%。这个经历揭示了一个深层规律：势函数的外推能力取决于训练数据的”覆盖度”——不仅要覆盖平衡态，还要覆盖非平衡态的极端构型（位错、晶界、表面重构）。对于通用力场（如ReaxFF、COMB），训练数据可能涉及数千个构型，拟合成本更高，但可移植性更强。

解决验证：三类势函数的开发参数与验证

经过Cu-EAM、Al₂O₃-Buckingham和SiC-ReaxFF三类体系的开发，项目组建立了势函数参数化的规范流程：

验证结果：

– Cu-EAM：a₀偏差0.2%，C11偏差3%，C12偏差2%，C44偏差8%（EAM对剪切模量描述弱），E_vac偏差8%，熔点偏差4%

– Al₂O₃-Buckingham（Born-Mayer）：a₀偏差0.5%，C11偏差5%，表面能偏差12%（离子势对表面极化描述不足），声子谱在800 cm⁻¹以上偏差明显

– SiC-ReaxFF：a₀偏差1.2%，C11偏差15%（ReaxFF精度的代价），但化学反应能垒（SiC氧化）偏差仅5%，证明了ReaxFF在反应动力学上的优势

反思边界：多体势的精度天花板与数据饥渴

当前势函数开发存在一个根本性矛盾：精度越高，数据需求越大，可移植性越窄。EAM势在纯金属上可以达到DFT精度的90%（对于平衡态），但合金体系可能降至70-80%。ReaxFF可以描述化学反应，但弹性常数偏差可能达15-20%。机器学习势（MTP、ACE、SNAP）在理论上可以逼近DFT精度，但需要1000-10000个训练构型，每个构型都需要DFT计算，总成本可能超过百万核时。项目组的策略是：对于结构-力学性质的预测（如弹性、扩散、蠕变），使用EAM/MEAM势；对于化学反应预测（如氧化、腐蚀、催化），使用ReaxFF或ML势；对于高精度局部模拟（如裂纹尖端、位错芯），使用DFT+ML势的混合方案。温度效应（有限温度势）和量子效应（零点能、隧道效应）在当前势函数框架中未纳入，对于轻元素（H、He）在低温下的扩散，需要额外的量子修正。如需势函数参数化或分子动力学模拟服务，请访问科研学术网首页，或返回分子动力学栏目了解力场选型与模拟流程。