| 查看: 2016 | 回复: 4 | ||
fwforest金虫 (小有名气)
|
[求助]
100金币求助msi2lmp导出LAMMPS数据文件的问题
|
|
我想用LAMMPS和CLAYFF力场做粘土矿物的分子模拟。在MS中首先建好结构，并根据CLAYFF力场中原子的type分配好，然后导出.car和.mdf文件，但是我用msi2lmp转换data file时遇到了问题，报的错误如下： -------------------------------------------------------------------------------------------------- Running msi2lmp v3.9.1 / 08-Oct-2013 ... Forcefield: Class I Forcefield file name: ../frc_files/clayff.frc Reading car file: quartz_alpha.car Reading mdf file: quartz_alpha.mdf Building internal coordinate lists Reading forcefield file Get force field parameters for this system Unable to find equivalent type for ob Unable to find equivalent type for st Unable to find bond data for ob st ---------------------------------------------------------------------------------------------------- 但是实际上，我将自己导出的.car, .mdf文件和msi2lmp自带的例子做了一些比较，没看出来什么问题。但是例子可以运行，我导出的文件就运行不了，期待大神帮我解决一下，非常感谢！ 另外，如果模拟的时候要实现CLAYFF力场和CVFF力场的耦合该怎么做呢？ 我导出的.car文件： ----------------------------------------------------------------------------------------------------- !BIOSYM archive 3 PBC=ON Materials Studio Generated CAR File !DATE Wed May 21 21:42:07 2014 PBC 4.9100 4.9100 5.4020 90.0000 90.0000 120.0000 (P1) O1 -0.280434837 1.762907387 0.629330000 XXXX 1 ob O 0.000 O2 1.068492788 3.127863797 4.230663333 XXXX 1 ob O 0.000 O3 -0.788057951 3.613591183 2.429996667 XXXX 1 ob O 0.000 O4 1.666940000 0.638590000 4.772670000 XXXX 1 ob O 0.000 O5 2.174563113 2.489273797 1.171336667 XXXX 1 ob O 0.000 O6 3.523490739 1.124317387 2.972003333 XXXX 1 ob O 0.000 Si1 1.180315677 2.044366722 0.000000000 XXXX 1 st Si 0.000 Si2 2.549364547 0.000000000 3.601333333 XXXX 1 st Si 0.000 Si3 -1.274682273 2.207814461 1.800666667 XXXX 1 st Si 0.000 end end --------------------------------------------------------------------------------------------------------- 导出的.mdf文件： !BIOSYM molecular_data 4 !Date: Wed May 21 21:42:07 2014 Materials Studio Generated MDF file #topology @column 1 element @column 2 atom_type @column 3 charge_group @column 4 isotope @column 5 formal_charge @column 6 charge @column 7 switching_atom @column 8 oop_flag @column 9 chirality_flag @column 10 occupancy @column 11 xray_temp_factor @column 12 connections @molecule quartz_alpha XXXX_1:O1 O ob ? 0 0 0.0000 0 0 8 1.0000 0.0000 Si1 Si3 XXXX_1:O2 O ob ? 0 0 0.0000 0 0 8 1.0000 0.0000 Si2 Si1 XXXX_1:O3 O ob ? 0 0 0.0000 0 0 8 1.0000 0.0000 Si3 Si2 XXXX_1:O4 O ob ? 0 0 0.0000 0 0 8 1.0000 0.0000 Si1 Si2 XXXX_1:O5 O ob ? 0 0 0.0000 0 0 8 1.0000 0.0000 Si3 Si1 XXXX_1:O6 O ob ? 0 0 0.0000 0 0 8 1.0000 0.0000 Si2 Si3 XXXX_1:Si1 Si st ? 0 0 0.0000 0 0 8 1.0000 0.0000 O1 O4 O5 O2 XXXX_1:Si2 Si st ? 0 0 0.0000 0 0 8 1.0000 0.0000 O2 O6 O4 O3 XXXX_1:Si3 Si st ? 0 0 0.0000 0 0 8 1.0000 0.0000 O3 O5 O6 O1 ! #symmetry @periodicity 3 xyz @group (P1) #end ---------------------------------------------------------------------------------------------------------------- |
回复此楼» 猜你喜欢
焦虑
已经有8人回复
-
308求调剂
已经有4人回复
-
NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗
已经有14人回复
-
材料与化工一志愿南昌大学327求调剂推荐
已经有6人回复
-
化学调剂0703
已经有7人回复
-
327求调剂
已经有11人回复
-
调剂
已经有8人回复
-
梁成伟老师课题组欢迎你的加入
已经有7人回复
-
伙伴们，祝我生日快乐吧
已经有24人回复
-
中科院材料273求调剂
已经有3人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
fwforest
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 3619.4
- 散金: 780
- 红花: 2
- 帖子: 270
- 在线: 399.4小时
- 虫号: 846321
- 注册: 2009-09-12
- 性别: GG
- 专业: 水文地质学(含地热地质学)
2楼2014-05-22 16:09:15
hooge
木虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 1790.1
- 散金: 236
- 红花: 1
- 帖子: 975
- 在线: 135.1小时
- 虫号: 89636
- 注册: 2005-08-29
- 专业: 金属材料的力学行为
3楼2014-05-22 16:16:00
fwforest
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 3619.4
- 散金: 780
- 红花: 2
- 帖子: 270
- 在线: 399.4小时
- 虫号: 846321
- 注册: 2009-09-12
- 性别: GG
- 专业: 水文地质学(含地热地质学)
4楼2014-05-22 16:56:19
fwforest
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 3619.4
- 散金: 780
- 红花: 2
- 帖子: 270
- 在线: 399.4小时
- 虫号: 846321
- 注册: 2009-09-12
- 性别: GG
- 专业: 水文地质学(含地热地质学)
5楼2014-05-23 14:57:33
5