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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » LAMMPS/DL_POLY » lammps计算的能量有误？请大家帮我分析下是为什么
111/2返回列表12下一页
查看: 2185  |  回复: 10
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

gtssongchi

木虫 (小有名气)

[求助] lammps计算的能量有误？请大家帮我分析下是为什么 已有3人参与

大家好，本人lammps初学，找了一个bcc-Fe的结构练练手，算了下体系能量lammps得到的能量如下
Step       Temp   E_pair           E_mol    TotEng            Press
       0            0   -231.65991            0   -231.65991    -91.492453

而之前我用VASP也算过同样体系，得到的能量为
  FREE ENERGIE OF THE ION-ELECTRON SYSTEM (eV)
  ---------------------------------------------------
  free  energy   TOTEN  =      -448.766382 eV

  energy  without entropy=     -448.747734  energy(sigma->0) =     -448.760166



不知道为什么会有这么大的差别，lammps中用了两种势函数计算，得出的结果相差不大。
回复此楼

» 猜你喜欢
1楼 2017-09-20 09:01:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

匿名

用户注销 (正式写手)
本帖仅楼主可见
赞一下
5楼2017-09-20 10:37:10
已阅   申请模拟EPI   回复此楼   编辑   查看我的主页

cl1994861

金虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与，应助指数 +1
gtssongchi: 金币+5, ★★★很有帮助, 感谢 2017-09-20 12:47:05
首先上面的lammps你没说单位制是metal还是real还是其他什么的，单位不一样结果也不一样，还有能量跟温度有关的，你上面的显示的温度是0，当然能量也跟其他温度的不同，望有用
赞一下(1人) 回复此楼
2楼2017-09-20 09:32:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

匿名

用户注销 (正式写手)
本帖仅楼主可见
赞一下
3楼2017-09-20 10:33:31
已阅   申请模拟EPI   回复此楼   编辑   查看我的主页

匿名

用户注销 (正式写手)
★ ★ ★ ★ ★
感谢参与，应助指数 +1
gtssongchi: 金币+5 2017-09-20 12:54:05
本帖仅楼主可见
赞一下
4楼2017-09-20 10:35:44
已阅   申请模拟EPI   回复此楼   编辑   查看我的主页

gtssongchi

木虫 (小有名气)
引用回帖:
2楼: Originally posted by cl1994861 at 2017-09-20 09:32:09
首先上面的lammps你没说单位制是metal还是real还是其他什么的，单位不一样结果也不一样，还有能量跟温度有关的，你上面的显示的温度是0，当然能量也跟其他温度的不同，望有用

你好，单位设置为metal，我在lammps中做的是run 0，也就是直接计算体系能量。而vasp中也是静态计算的能量，所以我当时以为它们的结果是可以比较的。
赞一下 回复此楼
6楼2017-09-20 12:48:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gtssongchi

木虫 (小有名气)
引用回帖:
5楼: Originally posted by xmc057711 at 2017-09-20 10:37:10
所以绝对值比较是没有意义的

感谢，我想可能的确是这个原因，请问“lammps单原子势能为0”是在手册中哪里有介绍呢？
赞一下 回复此楼
7楼2017-09-20 12:53:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

用户注销 (正式写手)
本帖仅楼主可见
赞一下
8楼2017-09-20 13:32:49
已阅   申请模拟EPI   回复此楼   编辑   查看我的主页

匿名

用户注销 (正式写手)
本帖仅楼主可见
赞一下
9楼2017-09-20 13:34:41
已阅   申请模拟EPI   回复此楼   编辑   查看我的主页

abinitio

版主 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与，应助指数 +1
gtssongchi: 金币+5, 谢谢交流 2017-09-20 14:15:30
引用回帖:
5楼: Originally posted by xmc057711 at 2017-09-20 10:37:10
所以绝对值比较是没有意义的

同意5楼意见，绝对值没有意义，经典md和abnitio md二者的机理不同，初始的假设也不同，近似不同。
赞一下 回复此楼
10楼2017-09-20 14:01:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 gtssongchi 的主题更新
111/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 vZfe6xYu34yj 2026-09-14 3/150 2026-09-16 08:51 by SjQobnC04j06
[博后之家] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 vZfe6xYu34yj 2026-09-14 3/150 2026-09-16 08:27 by SjQobnC04j06
[找工作] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-16 05:15 by LH5NkK5Ud4bh
[考研] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +5 6F5UbRU2I5hL 2026-09-14 5/250 2026-09-16 05:03 by LH5NkK5Ud4bh
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 s3fFTmArrBt6 2026-09-14 3/150 2026-09-16 04:23 by DgNGHc3h5tPl
[博后之家] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-14 4/200 2026-09-16 04:11 by DgNGHc3h5tPl
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +4 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-16 03:23 by DgNGHc3h5tPl
[考研] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +6 QUjhNVAcOSff 2026-09-13 6/300 2026-09-16 02:02 by DgNGHc3h5tPl
[考博] 申博发邮件 +4 Lmengk 2026-09-15 4/200 2026-09-16 00:05 by 木心想继续深造
[公派出国] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 Aj1rhIDL5ixY 2026-09-14 3/150 2026-09-15 19:59 by CgyNCDVNhGVg
[找工作] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-15 18:47 by CgyNCDVNhGVg
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +4 23jxep3nCNZb 2026-09-14 4/200 2026-09-15 18:11 by CgyNCDVNhGVg
[考研] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +5 s3fFTmArrBt6 2026-09-14 5/250 2026-09-15 16:59 by CgyNCDVNhGVg
[考博] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +6 s3fFTmArrBt6 2026-09-14 7/350 2026-09-15 16:47 by CgyNCDVNhGVg
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-14 5/250 2026-09-15 16:39 by CgyNCDVNhGVg
[公派出国] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +4 s3fFTmArrBt6 2026-09-14 4/200 2026-09-15 16:26 by CgyNCDVNhGVg
[公派出国] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +7 QUjhNVAcOSff 2026-09-13 9/450 2026-09-15 14:10 by VMDqgnlKimtZ
[教师之家] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 07:50 by C79jjtjAKjEn
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 05:38 by 5BDX0d0WFp7t
[考研] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +5 LwdutQ8HoqWP 2026-09-13 5/250 2026-09-14 22:25 by vZfe6xYu34yj
信息提示
关闭
请填处理意见
关闭