| 查看: 2151 | 回复: 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 ---------------------------------------------------------------------------------------------------------------- |
回复此楼» 猜你喜欢
2026国自然函评费到账
已经有17人回复
-
让我中一个面上吧！
已经有15人回复
-
2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战
已经有5人回复
-
跳槽后在研项目怎么办？
已经有10人回复
-
什么时候开奖？
已经有10人回复
-
今天放榜吗？
已经有16人回复
-
93BebMhtakh前后11位开头都是大写
已经有8人回复
-
范进中举一文的中心思想
已经有4人回复
-
只有每年这种时候来逛逛小木虫
已经有26人回复
-
filecode，4个jtjc了
已经有16人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
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
- 专业: 水文地质学(含地热地质学)
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
- 专业: 水文地质学(含地热地质学)
5楼2014-05-23 14:57:33
40