大家好，我是小马老师。
本文介绍lammps模拟中pcff势设置方法。
cvff、pcff是ms文件转换为lammps data文件最常用的两种势。
相比于cvff势，pcff势参数更多，但在设置方式上和cvff势过程是一样的。
本文以沥青材料为例，介绍pcff势设置方法。
(1) 建模并转换为data文件
在ms中使用AC模块建立沥青模型，使用forcite模块设置pcff力场，导出为asphalt.car。
使用msi2lmp转换为data文件：
msi2lmp asphalt -class II -frc pcff -i >data.asphalt
转换之后得到asphalt.data，用记事本打开asphalt.data，找到势函数部分：
Pair Coeffs # lj/class2/coul/long
1 0.0640000000 4.0100000000 # cp
2 0.0710000000 4.0270000000 # s
………………
16 0.0130000000 1.0980000000 # ho
Bond Coeffs # class2
1 1.4170 470.8361 -627.6179 1327.6345 # cp-cp
………………
38 0.9650 532.5062 -1282.9050 2004.7658 # oh-ho
Angle Coeffs # class2
1 118.9000 61.0226 -34.9931 0.0000 # cp-cp-cp
2 120.0500 44.7148 -22.7352 0.0000 # cp-cp-c1
…………
119 108.1900 53.1250 -8.5016 0.0000 # cp-oh-ho
Dihedral Coeffs # class2
1 8.3667 0.0000 1.1932 0.0000 0.0000 0.0000# cp-cp-cp-cp
2 0.0000 0.0000 3.9661 0.0000 0.0000 0.0000# cp-cp-cp-hc
因篇幅有限，只列出部分势参数。
pcff在lammps中对应的势函数多为class2类型。
pcff势对应的lj势为9-6势，如果需要修正参数的话，需要注意参数的转换。
(2) in文件的写法
data文件中势函数部分已经把全部的参数列出，在in文件中只需写明势的类型即可，势的类型就是data中各种势“#”后面的名称。
pcff势的设置语句要放到read_data命令的前面。
pcff默认的pair势有长程库仑力，因此需要设置kspace_style。
in文件pcff势具体设置为：
pair_style lj/class2/coul/long 10
bond_style class2
angle_style class2
dihedral_style class2
improper_style class2
kspace_style pppm 1e-4
read_data asphalt.data
读取文件之后，如果不需要替换参数，直接就可以进行弛豫计算。
公众号定期分享lammps命令和模拟技巧，感谢分享和关注^_^。
END
公众号案例代码以及学习交流，请加QQ群：754749935
点击关注微信公众号：lammps加油站
