转载:https://www.cnblogs.com/zll-52011/p/10960905.html
1.从美国矿物数据库下载矿物CIF(有晶格)
2.晶胞CIF导入MS
3.选择display style中的ball and stick视图
4.加键add bonds,看看有没有错。此步骤前根据需要是否补H
5.手动指定每个原子力场类型
edit-atom selection-比如选好H(先选不用区分类型的元素)
edit-edit sets-new(定义一个集)
View-Explorers-Properties Explorer-Filter选atom-改Forcefield Type(先点击那个原子集)
其他元素一样,对于要区分类型的,如O,自己判断ob与oh(与H连),atom selection时一个个按ctrl选择好,再edit sets.
6.所有键加好,原子类型定义好后,扩胞
Build-symmetry-supercell
7.这里MS的优化模块,无论时GULP还是Forcite模块力场经过我不懈的努力还是没法用,我们直接导出文件
File- Export-选择car.cor格式(会导出两个文件)
8.在服务器账号里创建一个文件夹kao(文件夹名称一定要和你导出来的两个文件名字一样,区分大小写),将kao.car和kao.cor放进去
9.将frc_files放在kao文件夹同一个目录
10.已编译msi2lmp.exe(珠海服务器)
11.进入工作目录kao
cd kao/
msi2lmp.exe kao -class I -frc clayff -ignore > out.log
生成后kao文件夹会有data生成,可检查一下data力场对不对,不对的话检查一下frc_files里面的clayff.frc参数准不准。
12.生成的data是没有电荷信息的,可以在VMD里设置并微调为0,或者在in文件中设置,也需要平衡电荷。