lammps：pcff势函数参数设置方法

大家好，我是小马老师。
本文介绍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过程中遇到什么问题，可到公众号或者QQ群留言。
公众号案例代码以及学习交流，请加QQ群：754749935

扫描关注微信公众号：lammps加油站。