应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » LAMMPS/DL_POLY » 为什么两个碰撞后能量发散了?
51/1返回列表
查看: 1584  |  回复: 4
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

mathsphy

新虫 (小有名气)

[求助] 为什么两个碰撞后能量发散了? 已有2人参与

大家好,我初学lammps,今天做两个铁原子碰撞时程序提示有原子丢失,所以我就使用了可变时间步长,结果发现它系统的能量不太对,所以想请教一下.
我这个系统就只建了一个5*5*5的晶胞,然后放了三个原子,坐标分别是(0,0,0),(0,3,3),(0,3,0),晶格单位.
第一个原子是用作参考的,第二个原子碰撞第三个原子.系综和设置如下:
fix                                1 all nve
fix                                2 all dt/reset 1 NULL 0.001 0.01
velocity                1 set 1000 0 0
velocity                2 set 0 0 -1000
velocity                3 set 0 0 0

输出的信息如下,其中最右边三列分别代表第1,2,3原子的动能:
Step Dt Time Temp TotEng PotEng KinEng 4 5 6
       0       1e-005            0 1.8351477e+008    47442.293            0    47442.293    23721.146    23721.146            0
      10       1e-005       0.0001 1.8351477e+008    47442.293            0    47442.293    23721.146    23721.146            0
      20       1e-005       0.0002 1.8351477e+008    47442.293            0    47442.293    23721.146    23721.146            0
      30       1e-005       0.0003 1.8351477e+008    47442.293            0    47442.293    23721.146    23721.146            0
      40       1e-005       0.0004 1.8351477e+008    47442.293            0    47442.293    23721.146    23721.146            0
      50       1e-005       0.0005 1.8351477e+008    47442.293            0    47442.293    23721.146    23721.146            0
      60       1e-005       0.0006 1.8351477e+008    47442.293            0    47442.293    23721.146    23721.146            0
      70       1e-005       0.0007 1.8351477e+008    47442.293            0    47442.293    23721.146    23721.146            0
      80       1e-005       0.0008 1.8351477e+008    47442.293            0    47442.293    23721.146    23721.146            0
      90       1e-005       0.0009 1.8351477e+008    47442.293            0    47442.293    23721.146    23721.146            0
     100       1e-005        0.001 1.8351477e+008    47442.293            0    47442.293    23721.146    23721.146            0
     110       1e-005       0.0011 1.8351477e+008    47442.293            0    47442.293    23721.146    23721.146            0
     120       1e-005       0.0012 1.8351477e+008    47442.293 -2.2147594e-005    47442.293    23721.146    23721.146 5.3606842e-015
     130 9.9999999e-006       0.0013 1.8351478e+008    47442.293 -0.00057027847    47442.293    23721.146    23721.147 3.4424435e-012
     140 9.9999994e-006       0.0014 1.8351478e+008    47442.293 -0.002693678    47442.295    23721.146    23721.149 7.6714312e-011
     150 9.9999977e-006       0.0015 1.8351481e+008    47442.293 -0.010243151    47442.303    23721.146    23721.157 1.108068e-009
     160 9.9999835e-006 0.0015999999 1.8351504e+008    47442.293 -0.067726385    47442.361    23721.146    23721.215 4.9021806e-008
     170 9.9998924e-006 0.0016999994 1.8351669e+008    47442.293  -0.49445878    47442.787    23721.146    23721.641 2.5772336e-006
     180 9.9997992e-006 0.0017999979 1.8351836e+008    47442.293  -0.92702466     47443.22    23721.146    23722.074 9.0716298e-006
     190 9.9996775e-006 0.0018999954 1.8352057e+008    47442.293    -1.497937    47443.791    23721.146    23722.645 2.3667644e-005
     200 9.9995503e-006 0.0019999916 1.8352292e+008    47442.293   -2.1048515    47444.398    23721.146    23723.252 4.6713654e-005
     210 9.9994425e-006 0.0020999866 1.8352492e+008    47442.293   -2.6212737    47444.914    23721.146    23723.768 7.2437174e-005
     220 9.9993716e-006 0.0021999806 1.8352627e+008    47442.293   -2.9705227    47445.263    23721.146    23724.117 9.3009162e-005
     230 9.9993898e-006 0.0022999743 1.8352597e+008    47442.293   -2.8944083    47445.188    23721.146    23724.041 8.8322947e-005
     240 9.9997456e-006 0.0023999692 1.8351944e+008    47442.287   -1.2123728      47443.5    23721.146    23722.353 1.5352359e-005
     250       1e-005 0.0024999687 1.8347824e+008    47442.269    9.4212164    47432.847    23721.146      23711.7 0.00094058064
     260       1e-005 0.0025999687 1.8330903e+008     47442.22    53.116248    47389.104    23721.146    23667.928  0.029882795
     270       1e-005 0.0026999687 1.8275428e+008    47442.011    196.32033    47245.691    23721.146    23524.134   0.41077332
     280       1e-005 0.0027999687 1.8045156e+008     47440.57    790.17944     46650.39    23721.146    22922.405    6.8394779
     290 8.1905831e-012 0.0028299688 2.7355272e+020 7.0718929e+016    165.60334 7.0718929e+016    23721.146 3.5359435e+016 3.5359493e+016
     300 8.1905831e-012 0.0028299689 2.7355272e+020 7.0718929e+016    6.8461037 7.0718929e+016    23721.146 3.5359435e+016 3.5359493e+016
     310 8.1905831e-012 0.0028299689 2.7355272e+020 7.0718929e+016   -2.9300923 7.0718929e+016    23721.146 3.5359435e+016 3.5359493e+016
     320 8.1905831e-012  0.002829969 2.7355272e+020 7.0718929e+016    -2.556017 7.0718929e+016    23721.146 3.5359435e+016 3.5359493e+016
     330 8.1905831e-012 0.0028299691 2.7355272e+020 7.0718929e+016   -1.4125551 7.0718929e+016    23721.146 3.5359435e+016 3.5359493e+016
     340 8.1905831e-012 0.0028299692 2.7355272e+020 7.0718929e+016  -0.44704995 7.0718929e+016    23721.146 3.5359435e+016 3.5359493e+016
     350 8.1905831e-012 0.0028299693 2.7355272e+020 7.0718929e+016 -0.0086311757 7.0718929e+016    23721.146 3.5359435e+016 3.5359493e+016

我不明白为啥会发散啊.....请问要怎么解决这个问题?...
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
    本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
  • 附件 1 : 1.lammps
    • 2015-11-28 23:12:36, 71.03 K

» 收录本帖的淘帖专辑推荐

分子动力学

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2015-11-28 23:14:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luftf

木虫 (正式写手)

【答案】应助回帖

感谢参与，应助指数 +1
时间步长太大了吧？
或者你把输出步数间隔设置短一点
赞一下(1人) 回复此楼
2楼2015-11-29 20:18:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mathsphy

新虫 (小有名气)
引用回帖:
2楼: Originally posted by luftf at 2015-11-29 20:18:00
时间步长太大了吧？
或者你把输出步数间隔设置短一点

我是用可变时间步长的,后面它的步长都自动调到1e-11了..
赞一下(2人) 回复此楼
3楼2015-11-30 10:08:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feng0603

新虫 (小有名气)

【答案】应助回帖

楼主！原子丢失，是不是与边界条件的设置有关啊？如果用f边界，原子出去了就回不来了！
赞一下 回复此楼
4楼2015-12-10 11:02:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mathsphy

新虫 (小有名气)
引用回帖:
4楼: Originally posted by feng0603 at 2015-12-10 11:02:16
楼主！原子丢失，是不是与边界条件的设置有关啊？如果用f边界，原子出去了就回不来了！

我用的是周期性边界，应该不是晶界的问题吧。。

发自小木虫Android客户端
赞一下 回复此楼
5楼2015-12-10 23:39:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 mathsphy 的主题更新
51/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 Vi50GxzrFcSG 2026-08-07 6/300 2026-08-08 10:19 by 3OOjAIS77qg2
[基金申请] fileCode有新解读？ +8 Tide man 2026-08-08 16/800 2026-08-08 09:59 by Tide man
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 4/200 2026-08-08 08:27 by 6vVgjDL4CnGu
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 3/150 2026-08-08 03:47 by 6vVgjDL4CnGu
[基金申请] 娱乐 +6 Tide man 2026-08-03 6/300 2026-08-07 22:40 by 铁帽子农民
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变，有变的大神么？ +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] 关于代码变化问题，想知道的进来 +9 且听虎啸 2026-08-07 12/600 2026-08-07 21:53 by 理想三旬2025
[基金申请] 固定端突然变了，今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] 听说今天filecode变了 +24 布布和一二 2026-08-06 47/2350 2026-08-07 16:02 by zhiyanjiang
[基金申请] filecode与中标关系的预测 +3 布布和一二 2026-08-07 3/150 2026-08-07 15:09 by gltch
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[论文投稿] 十年后又回来了，论文投稿求助 +4 哈哈114477 2026-08-01 4/200 2026-08-07 14:39 by jgy194592
[基金申请] 大家散了吧，后缀研究没有意义，别浪费时间了，过好目前的每一天，不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] 关于豆爷回答的JTJC与%2F数量 +4 yang182083 2026-08-06 6/300 2026-08-07 12:43 by 虫友是什么虫
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
[基金申请] 8月时间戳变的，举个手。玩一下，释放压力 +9 archvillain 2026-08-04 11/550 2026-08-05 20:06 by wlwhappy
[基金申请] 好消息？这个有何含义？？？ +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
信息提示
关闭
请填处理意见
关闭