| 查看: 4959 | 回复: 22 | ||
胡博洋的洋金虫 (小有名气)
|
[求助]
经过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（即不考虑某些作用角）时却可以正常运行，计算出结果，如果我想考虑全部的键角就会错误，所以在此恳求各位前辈，帮我看看可以怎么解决感激不尽~！ |
回复此楼» 猜你喜欢
308求调剂
已经有4人回复
-
NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗
已经有14人回复
-
材料与化工一志愿南昌大学327求调剂推荐
已经有6人回复
-
化学调剂0703
已经有7人回复
-
327求调剂
已经有11人回复
-
调剂
已经有8人回复
-
梁成伟老师课题组欢迎你的加入
已经有7人回复
-
伙伴们，祝我生日快乐吧
已经有24人回复
-
中科院材料273求调剂
已经有3人回复
-
材料工程专硕274一志愿211求调剂
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 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人回复
7楼2014-05-04 11:02:15
10楼2014-09-22 17:14:47
139564@chen
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 2170.6
- 红花: 1
- 帖子: 134
- 在线: 227.3小时
- 虫号: 1607763
- 注册: 2012-02-09
- 性别: GG
- 专业: 理论和计算化学
2楼2013-12-29 16:55:42
胡博洋的洋
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1157.8
- 帖子: 93
- 在线: 123.9小时
- 虫号: 2351552
- 注册: 2013-03-16
- 性别: GG
- 专业: 理论和计算化学
3楼2013-12-29 20:16:47
ashioulou
铜虫 (初入文坛)
- 应助: 5 (幼儿园)
- 金币: 167.9
- 帖子: 23
- 在线: 35.6小时
- 虫号: 2536624
- 注册: 2013-07-08
- 性别: GG
- 专业: 聚合物共混与复合材料
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
胡博洋的洋: 金币+10, ★★★很有帮助, 谢谢~，最近较忙，一直没去试 2014-01-02 10:03:51
月只蓝: 金币+1, 感谢指导。 2014-01-02 10:42:06
胡博洋的洋: 金币+10, ★★★很有帮助, 谢谢~，最近较忙，一直没去试 2014-01-02 10:03:51
月只蓝: 金币+1, 感谢指导。 2014-01-02 10:42:06
|
??????N????}???????????????^???????????y?\????? ?????????N??in?????????^???????????N???} ????????Q???k?????? 1. ?S??minimize???????y??????С?????\????C 2. ??_???timestep?{С??c????Εr?g |
4楼2013-12-30 10:00:42
5楼2014-01-02 09:31:30
wdwangcn
金虫 (著名写手)
- 应助: 26 (小学生)
- 金币: 1443.6
- 散金: 2430
- 红花: 8
- 帖子: 1044
- 在线: 783.5小时
- 虫号: 674141
- 注册: 2008-12-13
- 性别: GG
- 专业: 微/纳机械系统
6楼2014-05-02 18:27:33
8楼2014-09-22 10:43:34
9楼2014-09-22 11:05:25
