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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » LAMMPS/DL_POLY » 为什么两个碰撞后能量发散了?
北京石油化工学院2026年研究生招生接收调剂公告
51/1返回列表
查看: 1382  |  回复: 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的回帖

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的回帖
查看全部 5 个回答

luftf

木虫 (正式写手)

【答案】应助回帖

感谢参与，应助指数 +1
时间步长太大了吧？
或者你把输出步数间隔设置短一点
赞一下(1人) 回复此楼
2楼2015-11-29 20:18:00
已阅   回复此楼   关注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的回帖
查看全部 5 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085404 22408 315分 +4 zhuangyan123 2026-03-31 5/250 2026-03-31 13:20 by zhuangyan123
[考研] 一志愿大连理工大学，机械工程学硕，341 +3 西瓜田的守望者 2026-03-30 3/150 2026-03-31 11:08 by asdfzly
[考研] 262求调剂 +7 ZZ..000 2026-03-30 8/400 2026-03-31 10:05 by cal0306
[考研] 材料科学与工程求调剂 +10 深V宿舍吧 2026-03-29 10/500 2026-03-31 04:08 by fmesaito
[考研] 本科211总分289，08工学真心求调剂 +3 utopiaE 2026-03-30 3/150 2026-03-30 23:42 by ms629
[考研] 085600 材料与化工 329分求调剂 +18 Mr. Z 2026-03-25 19/950 2026-03-30 19:01 by Wang200018
[考研] 293求调剂 +3 末未mm 2026-03-30 5/250 2026-03-30 17:23 by 王保杰33
[考研] 求调剂323材料与化工 +10 1124361 2026-03-24 10/500 2026-03-30 16:26 by 690616278
[考研] 329求调剂，一志愿西北工业大学，材料工程（085601） +5 小小机灵虫 2026-03-29 11/550 2026-03-30 15:02 by Wang200018
[硕博家园] 求调剂 有机化学考研356分 +10 Nadiums 2026-03-25 11/550 2026-03-30 11:56 by yjolah
[考研] 318一志愿吉林大学生物与医药 求调剂 +5 笃行致远. 2026-03-28 5/250 2026-03-30 06:56 by ilovexiaobin
[考研] 一志愿北京工业大学，324分求调剂 +6 零八# 2026-03-28 6/300 2026-03-29 21:20 by nanaliuyun
[考研] 311求调剂 +6 冬十三 2026-03-24 6/300 2026-03-29 20:45 by 无际的草原
[考研] 343求调剂085601 +3 要努力学习x 2026-03-29 3/150 2026-03-29 18:35 by wxiongid
[考研] 085600，专业课化工原理，321分求调剂 +5 大馋小子 2026-03-28 5/250 2026-03-29 08:56 by qingfeng258
[考研] 311（085601）求调剂 +4 liziyeyeye 2026-03-28 4/200 2026-03-28 18:50 by 535743368
[考研] 复试调剂，一志愿南农083200食品科学与工程 +5 XQTJZ 2026-03-26 5/250 2026-03-27 14:49 by 狂炫麦当当
[考研] 085600，材料与化工321分调剂 +4 大馋小子 2026-03-27 6/300 2026-03-27 14:11 by 松花缸1201
[考研] 一志愿郑大085600，310分求调剂 +5 李潇可 2026-03-26 5/250 2026-03-27 11:14 by 不吃魚的貓
[考研] 求调剂 一志愿 本科 北科大 化学 343 +6 13831862839 2026-03-24 7/350 2026-03-26 22:57 by 不吃魚的貓
信息提示
关闭
请填处理意见
关闭