24小时热门版块排行榜

返回列表

胡博洋的洋

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 1157.8
帖子: 93
在线: 123.9小时
虫号: 2351552
注册: 2013-03-16
性别: GG
专业: 理论和计算化学

[求助] 经过MS模型导入到LAMMPS，在运算时出现的错误求助已有2人参与

我的模型是从MS软件中导出来的，力场设置的是cvff，用LAMMPS中的msi2lmp工具导出data文件。data文件我也认真看了，没有错误。但是我在提交LAMMPS运算时出现了如下情况：
------------------------------------------------------------------------------------------------------
LAMMPS (22 Mar 2013)
Scanning data file ...
  25 = max angles/atom
  88 = max dihedrals/atom
Reading data file ...
  orthogonal box = (0 0 0) to (10.4507 6.0874 4.752)
  4 by 2 by 2 MPI processor grid
  28 atoms
  116 angles
  272 dihedrals
Finding 1-2 1-3 1-4 neighbors ...
  0 = max # of 1-2 neighbors
  0 = max # of 1-3 neighbors
  0 = max # of 1-4 neighbors
  1 = max # of special neighbors
4 atoms in group 1
16 atoms in group 2
4 atoms in group 3
4 atoms in group 4
Replicating atoms ...
  orthogonal box = (0 0 0) to (31.3521 18.2622 14.256)
  4 by 2 by 2 MPI processor grid
  756 atoms
  3132 angles
  7344 dihedrals
Finding 1-2 1-3 1-4 neighbors ...
  0 = max # of 1-2 neighbors
  0 = max # of 1-3 neighbors
  0 = max # of 1-4 neighbors
  1 = max # of special neighbors
WARNING: Resetting reneighboring criteria during minimization (min.cpp:173)
Ewald initialization ...
  G vector (1/distance) = 0.284516
  estimated absolute RMS force accuracy = 0.000147646
  estimated relative force accuracy = 1.02535e-05
  KSpace vectors: actual max1d max3d = 731 10 4630
Setting up minimization ...
Memory usage per processor = 36.3555 Mbytes
Step Temp E_pair E_mol TotEng Press
   0          0 -29439.954 139654.04 110214.08 -1054599.1
   43          0 -3.625706e+16 88135.753 -3.625706e+16 -2.3722708e+18
Loop time of 1.14906 on 16 procs for 43 steps with 756 atoms

Minimization stats:
  Stopping criterion = linesearch alpha is zero
  Energy initial, next-to-last, final =
      110214.082772 -3.62570600766e+16 -3.62570600766e+16
  Force two-norm initial, final = 5892.7 1.29073e+31
  Force max component initial, final = 355.236 8.16331e+30
  Final line search alpha, max atom move = 9.9195e-74 8.0976e-43
  Iterations, force evaluations = 43 424

Pair  time (%) = 0.652155 (56.7553)
Bond  time (%) = 0.0531141 (4.62238)
Kspce time (%) = 0.35794 (31.1505)
Neigh time (%) = 0.00358602 (0.312082)
Comm  time (%) = 0.0253774 (2.20853)
Outpt time (%) = 0 (0)
Other time (%) = 0.0568924 (4.9512)

Nlocal: 47.25 ave 56 max 41 min
Histogram: 2 2 1 1 2 4 2 1 0 1
Nghost: 4239.88 ave 4298 max 4175 min
Histogram: 3 2 3 0 0 0 0 0 3 5
Neighs: 49139.7 ave 58021 max 42819 min
Histogram: 2 3 0 2 1 6 0 1 0 1

Total # of neighbors = 786235
Ave neighs/atom = 1039.99
Ave special neighs/atom = 0
Neighbor list builds = 2
Dangerous builds = 0
Ewald initialization ...
  G vector (1/distance) = 0.284516
  estimated absolute RMS force accuracy = 0.000147646
  estimated relative force accuracy = 1.02535e-05
  KSpace vectors: actual max1d max3d = 731 10 4630
Setting up run ...
Memory usage per processor = 36.5856 Mbytes
Step Temp E_pair E_mol TotEng Press Volume
   43          0 -3.625706e+16 88135.753 -3.625706e+16 -2.3722708e+18 8162.3914
APPLICATION TERMINATED WITH THE EXIT STRING: Segmentation fault (signal 11)
------------------------------------------------------------------------------------------------------
当我把minimize 命令去掉后，出现如下问题：
------------------------------------------------------------------------------------------------------
LAMMPS (22 Mar 2013)
Scanning data file ...
  25 = max angles/atom
  88 = max dihedrals/atom
Reading data file ...
  orthogonal box = (0 0 0) to (10.4507 6.0874 4.752)
  4 by 2 by 2 MPI processor grid
  28 atoms
  116 angles
  272 dihedrals
Finding 1-2 1-3 1-4 neighbors ...
  0 = max # of 1-2 neighbors
  0 = max # of 1-3 neighbors
  0 = max # of 1-4 neighbors
  1 = max # of special neighbors
4 atoms in group 1
16 atoms in group 2
4 atoms in group 3
4 atoms in group 4
Replicating atoms ...
  orthogonal box = (0 0 0) to (31.3521 18.2622 14.256)
  4 by 2 by 2 MPI processor grid
  756 atoms
  3132 angles
  7344 dihedrals
Finding 1-2 1-3 1-4 neighbors ...
  0 = max # of 1-2 neighbors
  0 = max # of 1-3 neighbors
  0 = max # of 1-4 neighbors
  1 = max # of special neighbors
Ewald initialization ...
  G vector (1/distance) = 0.284516
  estimated absolute RMS force accuracy = 0.000147646
  estimated relative force accuracy = 1.02535e-05
  KSpace vectors: actual max1d max3d = 731 10 4630
Setting up run ...
Memory usage per processor = 36.7132 Mbytes
Step Temp E_pair E_mol TotEng Press Volume
   0          0 -29439.954 139654.04 110214.08 -1054599.1 8162.3914
   100 17.156213 -27586.553 82839.726 55254.847 189936.49 7175.6492
ERROR on proc 4: Angle atoms 35 54 45 missing on proc 4 at step 136 (neigh_bond.cpp:168)
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 4
ERROR on proc 1: Angle atoms 511 530 521 missing on proc 1 at step 136 (neigh_bond.cpp:168)
ERROR on proc 5: Angle atoms 539 558 549 missing on proc 5 at step 136 (neigh_bond.cpp:168)
ERROR on proc 6: Angle atoms 203 222 213 missing on proc 6 at step 136 (neigh_bond.cpp:168)
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 1
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 5
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 6
ERROR on proc 0: Angle atoms 7 26 17 missing on proc 0 at step 136 (neigh_bond.cpp:168)
ERROR on proc 9: Angle atoms 632 643 627 missing on proc 9 at step 136 (neigh_bond.cpp:168)
ERROR on proc 11: Angle atoms 716 727 711 missing on proc 11 at step 136 (neigh_bond.cpp:168)
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 0
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 9
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 11
ERROR on proc 3: Angle atoms 679 698 689 missing on proc 3 at step 136 (neigh_bond.cpp:168)
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 3
ERROR on proc 8: Angle atoms 380 391 375 missing on proc 8 at step 136 (neigh_bond.cpp:168)
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 8
ERROR on proc 12: Angle atoms 72 83 67 missing on proc 12 at step 136 (neigh_bond.cpp:168)
ERROR on proc 13: Angle atoms 576 587 571 missing on proc 13 at step 136 (neigh_bond.cpp:168)
ERROR on proc 2: Angle atoms 175 194 185 missing on proc 2 at step 136 (neigh_bond.cpp:168)
ERROR on proc 10: Angle atoms 464 475 459 missing on proc 10 at step 136 (neigh_bond.cpp:168)
ERROR on proc 14: Angle atoms 240 251 235 missing on proc 14 at step 136 (neigh_bond.cpp:168)
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 2
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 12
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 13
ERROR on proc 7: Angle atoms 707 726 717 missing on proc 7 at step 136 (neigh_bond.cpp:168)
ERROR on proc 15: Angle atoms 744 755 739 missing on proc 15 at step 136 (neigh_bond.cpp:168)
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 7
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 10
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 14
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 15
------------------------------------------------------------------------------------------------------
当我删掉一些angles（即不考虑某些作用角）时却可以正常运行，计算出结果，如果我想考虑全部的键角就会错误，所以在此恳求各位前辈，帮我看看可以怎么解决感激不尽~！

回复此楼

» 猜你喜欢

真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗已经有6人回复
请问哪里可以有青B申请的本子可以借鉴一下。已经有3人回复
孩子确诊有中度注意力缺陷已经有14人回复
三甲基碘化亚砜的氧化反应已经有4人回复
请问下大家为什么这个铃木偶联几乎不反应呢已经有5人回复
请问有评职称，把科研教学业绩算分排序的高校吗已经有5人回复
2025冷门绝学什么时候出结果已经有3人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有4人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复
AI论文写作工具：是科研加速器还是学术作弊器？已经有3人回复

» 本主题相关价值贴推荐，对您同样有帮助:

lammps计算中出现原子距离太近现象已经有11人回复
LAMMPS计算无反应，nan 已经有4人回复
怎样把material studio中sio2的模型数据导入lammps中已经有9人回复
lammps计算结果xyz文件中每个数字代表的字符长短。已经有10人回复
LAMMPS并行计算的问题（cpu——time关系）已经有17人回复
用MS建模导入LAMMPS计算，是不是不能用compass力场？已经有9人回复
lammps求助：计算含有不同原子体系的热导率的in文件该怎么修改已经有7人回复
在LAMMPS中怎么建立图示的模型已经有6人回复
【讨论】lammps 用sw势计算硅的熔点怎么这么高（2400k),附in文件已经有13人回复

不得求之

1楼2013-07-16 10:17:57

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

回帖支持 ( 显示支持度最高的前 50 名 )

笑笑愁

金虫 (小有名气)

未知

应助: 23 (小学生)
金币: 709
散金: 10
红花: 4
帖子: 270
在线: 250.7小时
虫号: 1105596
注册: 2010-09-22
性别: GG
专业: 能源化工

【答案】应助回帖

★
月只蓝: 金币+1, 感谢指导！ 2014-05-16 18:31:19

引用回帖:

6楼: Originally posted by wdwangcn at 2014-05-02 18:27:33
请教vmd的方法，能给个链接么...

1、VMD可以从六维等网站下载，然后熟悉一下基本的VMD基本的操作就可以了，不用学的太复杂。
2、MS导出pdb格式，然后导入VMD，然后用topo writelammpsdata命令转化成lammps所需的data文件
3、lammps中使用dump ~ ~ atom ~ dump.lammpstrj命令产生轨迹文件，次文件可以直接导入VMD，并且有动态效果