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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » LAMMPS/DL_POLY » 为什么两个碰撞后能量发散了?
51/1返回列表
查看: 1626  |  回复: 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 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 01:34 by ZPa0EcMwuECS
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +5 ASdOkHsho7FD 2026-08-28 6/300 2026-08-30 01:03 by ZPa0EcMwuECS
[硕博家园] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +6 ASdOkHsho7FD 2026-08-28 10/500 2026-08-30 00:53 by ZPa0EcMwuECS
[论文投稿] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 ASdOkHsho7FD 2026-08-28 5/250 2026-08-30 00:41 by ZPa0EcMwuECS
[考研] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +7 ASdOkHsho7FD 2026-08-28 10/500 2026-08-30 00:41 by ZPa0EcMwuECS
[基金申请] 面上意见出来了 +8 黄鸟于飞Chao 2026-08-29 16/800 2026-08-29 20:45 by 黄鸟于飞Chao
[基金申请] 为什么到现在没收到通知？ +4 tannykie 2026-08-29 4/200 2026-08-29 17:14 by lmz0216
[教师之家] 导师吐槽：我怎么摊上了这么个极品研究生！ +9 苏东坡二世 2026-08-23 9/450 2026-08-29 14:41 by hustersqt
[基金申请] 麻烦专家们看看评委们的意见（F口面上） +5 gdd2018 2026-08-28 10/500 2026-08-29 14:10 by Jacob678
[基金申请] 为什么资助数各大高校都创新高，自己申请怎么就这么难 +12 Kittylucky 2026-08-27 13/650 2026-08-29 00:04 by superceng
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜？ 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国际合作可查了，中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 今天务委会开完了，明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 牛来！米来！面来！ +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 在坚冰还盖着北海的时候，我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
信息提示
关闭
请填处理意见
关闭