本文介绍利用Materials Studio建模后,转换成LAMMPS的data文件。

1. 建模。在MS中建立模型,如下图所示。这里建的是一个CH4分子在石墨烯上。

 2. 赋予力场。模型导入LAMMPS中拟采用pcff力场参数。故在MS中,先用Forcite模块中的力场模块进行力场赋值。具体方式为:打开Forcite模块中的Calculation —>点击Energy选项卡 —> Forcefield中选择pcff —>点击旁边的more,打开对话框如下图所示。取消勾选“Calculate automatically”,依次点击箭头所示是的两个Calculate(分别赋力场参数和电荷)。完成后,关闭对话框。

3. 导出。点击程序菜单栏中的“File”—>"Export",如下图所示。输出类型设置为.car,点击保存。

 

 4. 转换。利用lammps中自带的msi2lmp.exe进行转换,生成data文件。具体方式如下:

(1)找到msi2lmp.exe所在文件夹。如果没有配置环境变量,则需要将上一步生成的文件全部拷贝至这个文件夹中,并在此文件夹中进行后续操作。如果已经配置了环境变量,则在原来那个文件夹中操作即可。

(2)在文件夹的地址栏中,输入cmd,打开命令行,如下图所示。我的已经配置了环境变量,所以我没有把文件拷入msi2lmp.exe所在的文件夹中。

(3)输入命令:msi2lmp Graphene_CH4 -class 2 -frc pcff -i >data.dat

其中:Graphene_CH4 为MS导出的.car文件,-class后的2表示力场的类型,具体的可参考文件夹中的README文件。这里直接复制过来。

   -- -class  (or -c)
        # is the class of forcefield to use (I  or 1 = Class I e.g., CVFF)
                                            (O  or 0 = OPLS-AA)
                                            (II or 2 = Class II e.g., CFFx, pcff, COMPASS)

data.dat文件(名字可以随便取)记录了转换过程,如果文件最后出现“Normal program termination”,则表示转换成功。生成的data文件即包含了模型中各原子的力场信息参数。在LAMMPS的in文件中,直接用read_data命令即可调用。

END

Logo

旨在为数千万中国开发者提供一个无缝且高效的云端环境,以支持学习、使用和贡献开源项目。

更多推荐