| 查看: 2910 | 回复: 0 | |||
laughingjudy新虫 (初入文坛)
|
[求助]
求帮忙看一个TIP4P水分子in文件，提示Invalid atom type in create_atoms mol command
|
|
IN文件如下： # 3d TIP4P water model units real variable T equal 298.16 variable P equal 1 atom_style full dimension 3 boundary p p p molecule H2O h2o.txt lattice fcc 8 orient x 1 0 0 orient y 0 1 0 orient z 0 0 1 region box block 0 4 0 4 0 4 create_box 1 box bond/types 1 angle/types 1 create_atoms 1 box mol H2O 4928458 mass 3 15.9994 mass 4 1.008 bond_style harmonic bond_coeff 1 450 0.9572 angle_style harmonic angle_coeff 1 55 104.52 pair_style lj/cut/tip4p/cut 1 2 1 1 0.15 9 8.5 pair_coeff 3 3 0.1550 3.1536 neighbor 2.0 bin neigh_modify delay 5 every 1 velocity all create $T 102486 dist gaussion timestep 0.01 thermo 1000 fix 1 all shake 0.0001 100 1000 b 1 a 1 fix 2 all npt temp $T $T 100 iso $P $P 1000 drag 0.2 run 1000 h2o.txt如下： #water 3 atoms 2 bonds 1 angles Coords 1 0 0 0 2 0.9572 0 0 3 -0.9572*sin14.52 -0.9572*cos14.52 0 Types 1 1 2 2 3 2 Charges 1 -1.040 2 0.520 3 0.520 Bonds 1 1 1 2 2 1 1 3 Angles 1 1 2 1 3 Shake Flags 1 1 2 1 3 1 Shake Atoms 1 1 2 3 2 1 2 3 3 1 2 3 Shake Bond Types 1 1 1 1 2 1 1 1 3 1 1 1 log文件提示如下： LAMMPS (10 May 2014-ICMS) WARNING: OMP_NUM_THREADS environment is not set. (../comm.cpp:100) using 1 OpenMP thread(s) per MPI task # 3d TIP4P water model units real variable T equal 298.16 variable P equal 1 atom_style full dimension 3 boundary p p p molecule H2O h2o.txt lattice fcc 8 orient x 1 0 0 orient y 0 1 0 orient z 0 0 1 Lattice spacing in x,y,z = 8 8 8 region box block 0 4 0 4 0 4 create_box 1 box bond/types 1 angle/types 1 Created orthogonal box = (0 0 0) to (32 32 32) 1 by 1 by 1 MPI processor grid create_atoms 1 box mol H2O 4928458 ERROR: Invalid atom type in create_atoms mol command (../create_atoms.cpp:166) 不明白是少了语句呢 还是哪里写错了。请各位大神帮忙解答啊  |
回复此楼» 本帖@通知
@御剑江湖
» 猜你喜欢
梦想
已经有8人回复
-
有多少人是今天查系统知道结果的？
已经有19人回复
-
面上合作单位盖章
已经有9人回复
-
面上函评意见出来了，像什么等级？
已经有17人回复
-
为什么国自然不能直接公布
已经有5人回复
-
学科评审组评审是指会评吗？
已经有4人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有14人回复
-
小白求助 投论文要求的highlights应该如何写
已经有4人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有12人回复
-
怎么看青基中了没有啊
已经有6人回复
» 本主题相关商家推荐: (我也要在这里推广)
找到一些相关的精华帖子，希望有用哦~
求助：运行了gromacs，想查看体系的总原子数和水分子数！
已经有4人回复
- 求教给配位水分子加氢原子的方法！
已经有6人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
投票:
20