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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » LAMMPS/DL_POLY » LAMMPS输出的库伦能为什么是正值，而MS的discover算出的库伦能是负值
52/1返回列表
查看: 4396  |  回复: 19
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

zs806750658

铜虫 (正式写手)

[交流] LAMMPS输出的库伦能为什么是正值，而MS的discover算出的库伦能是负值 已有4人参与

相同的模型，相同的系宗，相同的力场，相同的参数设置，为什么用Material Studio的discover计算的单点能中库伦能是-4万多，而Lammps跑完输出的能量中库伦能是正5万多,看过手册发现两个软件计算库伦能的函数是一样的，为什么相差如此之大呢。而其他的范德华能、键能等相差不是很大，哪位高手能给予指点，不胜感激
回复此楼

» 猜你喜欢
1楼 2015-05-31 06:42:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zs806750658

铜虫 (正式写手)
引用回帖:
16楼: Originally posted by brigsun at 2015-12-16 18:20:59
你好！
我也遇到类似问题，但根据我的测试，我倾向于怀疑是lammps中special_bonds的设置问题。
默认情况下，lammps是不计算1-2,1-3,1-4原子的pairwise energy的。所以我怀疑这个设置MS和lammps不一样，但因为我查 ...

感谢回复，你的方法是正确的，lammps默认不计算1-4原子的pairwise energy, 但是MS里面默认是计算该能量的，至于MS如何设置，需要打开Forcite模块下的力场文件，里面可以设置，见附图。如果想关闭1-4原子的pairwise energy的计算，需要把COU_EXCLUDE_1-4后面的F改成T. 至于Discover模块，可以在相应模块里面的设置中修改。
LAMMPS输出的库伦能为什么是正值，而MS的discover算出的库伦能是负值
forcefield.jpg
赞一下(1人) 回复此楼
19楼2016-06-23 06:02:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 20 个回答

lsloneil

专家顾问 (正式写手)
★ ★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
月只蓝: 金币+3, 感谢指导 2015-06-02 13:05:32
你用的是什么pair_style， 又是怎么把静电能输出的？

如果你使用的相互作用势利用ewald summation或 pppm 来计算静电相互作用，那么静电相互作用被分为短程（实空间）和长程（倒空间）两部分进行计算。在用thermo_style进行输出时，ecoul对应于短程部分，通常是正值，elong对应于长程部分，通常是负值，两者之和才是总静电相互作用能。

类似的帖子
http://simulation.haotui.com/viewthread.php?tid=23854
赞一下(1人) 回复此楼
2楼2015-06-01 06:08:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zs806750658

铜虫 (正式写手)
引用回帖:
2楼: Originally posted by lsloneil at 2015-06-01 06:08:18
你用的是什么pair_style， 又是怎么把静电能输出的？

如果你使用的相互作用势利用ewald summation或 pppm 来计算静电相互作用，那么静电相互作用被分为短程（实空间）和长程（倒空间）两部分进行计算。在用therm ...

非常感谢你的指导，我的pair_style 选择的是lj/class2/coul/long 10.0，采用的是ewald summation计算的静电相互作用。thermol_style的静电输出包括ecoul和elong，其中ecoul=54507.739;elong=-9033.3904，两者之和还是正值45474.3486，但discover计算出的总势能却是负值-41512.690317。理论上来讲，两者的库伦能计算公式一样，而且参数设置一样，可是计算出来的结果却不同
赞一下 回复此楼
3楼2015-06-02 04:25:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lsloneil

专家顾问 (正式写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
3楼: Originally posted by zs806750658 at 2015-06-01 08:25:18
非常感谢你的指导，我的pair_style 选择的是lj/class2/coul/long 10.0，采用的是ewald summation计算的静电相互作用。thermol_style的静电输出包括ecoul和elong，其中ecoul=54507.739;elong=-9033.3904，两者之和还 ...

我也不是非常清楚为什么会出现这样的结果。我不是很清楚Material Studio里关于ewald summation的设置。你的体系是否是电中性，如果不是电中性的话有时可能会有奇怪的结果。

假如可能的话，可以做一个简单的测试：这个体系只包含一个正离子和一个负离子，两者所带电荷量相等，这样整个体系呈电中性。然后你用Material Studio和lammps计算相同一个configruation的静电能，看是否还会出现这样奇怪的现象。
赞一下(1人) 回复此楼
4楼2015-06-02 09:34:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 20 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 为什么国自然不能直接公布 +5 bjdxyxy 2026-08-26 5/250 2026-09-01 17:54 by 小伟大博士
[文学芳草园] 梦想 +5 myrtle 2026-08-26 7/350 2026-09-01 15:18 by myrtle
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 8/400 2026-09-01 12:58 by ssyjh
[基金申请] 基金未中，这种答复是模板吗？ +6 zhaosm1982 2026-08-27 7/350 2026-08-31 21:18 by qdxxmc
[基金申请] 哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识 +6 yuleib84 2026-08-26 7/350 2026-08-31 19:46 by 鱼翔浅底1
[基金申请] 面上意见出来了 +12 黄鸟于飞Chao 2026-08-29 23/1150 2026-08-31 18:57 by 黄鸟于飞Chao
[基金申请] 能否申诉？ +7 echo8914667 2026-08-30 8/400 2026-08-31 17:00 by yihongxu
[基金申请] 中青基了要发朋友圈吗？ +7 349506619 2026-08-28 7/350 2026-08-31 13:39 by 冼亮淀粉酶
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜？ 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？ +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 国际合作可查了，中了面上 (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
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 牛来！米来！面来！ +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
信息提示
关闭
请填处理意见
关闭