| 查看: 2019 | 回复: 1 | ||
[求助]
希望大佬们帮忙看一下这个invalid atom type in fix shake报错是为什么啊？ 已有1人参与
|
|
在data文件中mass已经都定义出来了原子类型了： 804 atoms 734 bonds 1144 angles 1626 dihedrals 120 impropers 10 atom types 12 bond types 24 angle types 31 dihedral types 3 improper types -0.063168202 29.936831798 xlo xhi -1.070661042 28.929338958 ylo yhi -0.531017684 29.468982316 zlo zhi Masses 1 15.999400 # o* 2 1.007970 # h* 3 1.007970 # h 4 12.011150 # c2 5 12.011150 # c' 6 15.999400 # o' 7 14.006700 # n2 8 1.007970 # hn 9 12.011150 # c1 10 12.011150 # c3 但是in文件中group了水分子也就是type 1 2，为什么fix shake报错？ units real atom_style full boundary p p p #kspace_style pppm 1e-4 #Out of range , can't compute PPPM kspace_style ewald 1e-4 ########################################### pair_style lj/cut/coul/long 10.0 10.0 pair_modify mix arithmetic ########################################### bond_style harmonic angle_style harmonic dihedral_style harmonic improper_style cvff read_data PAM.test thermo_style custom step fmax fnorm pe evdwl emol ecoul epair ebond eangle edihed elong thermo 100 group water type 1 2 group poly type 3 4 5 6 7 8 9 10 neighbor 2.0 bin neigh_modify delay 0 every 1 check yes timestep 0.5 #1.0 min_style sd min_modify dmax 0.2 minimize 1.0e-5 1.0e-6 2000 10000 reset_timestep 0 compute mol poly chunk/atom molecule compute 1 poly com/chunk mol compute 2 poly gyration/chunk mol compute 3 poly msd/chunk mol fix ave1 poly ave/time 1 100 1000 c_1[1] c_1[2] c_1[3] c_2 c_3[1] c_3[2] c_3[3] c_3[4] file msdandgyra.data mode vector fix 1 all nve fix 2 all temp/berendsen 300.0 300.0 10.0 fix 3 water shake 0.0001 20 1000 t 12 报错是这样的： ERROR: Invalid atom type index for fix shake (../fix_shake.cpp:144) Last command: fix 3 water shake 0.0001 20 1000 t 12 求大佬帮忙看一下  |
回复此楼» 猜你喜欢
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
-
2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？
已经有20人回复
-
今天务委会开完了，明天出结果吗
已经有22人回复
-
在坚冰还盖着北海的时候，我看到了怒放的梅花。
已经有9人回复
-
明天应该可查了！？
已经有6人回复
-
有没有大神帮我看看基金代码
已经有10人回复
-
放榜前的不淡定
已经有19人回复
» 本主题相关商家推荐: (我也要在这里推广)
2楼2018-08-09 16:50:04
28