| 查看: 3987 | 回复: 19 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
zs806750658
|
[交流]
LAMMPS输出的库伦能为什么是正值，而MS的discover算出的库伦能是负值 已有4人参与
|
|||
| 相同的模型，相同的系宗，相同的力场，相同的参数设置，为什么用Material Studio的discover计算的单点能中库伦能是-4万多，而Lammps跑完输出的能量中库伦能是正5万多,看过手册发现两个软件计算库伦能的函数是一样的，为什么相差如此之大呢。而其他的范德华能、键能等相差不是很大，哪位高手能给予指点，不胜感激 |
回复此楼» 猜你喜欢
333求调剂
已经有10人回复
-
307求调剂
已经有3人回复
-
考研化学308分求调剂
已经有4人回复
-
298求调剂
已经有7人回复
-
263求调剂
已经有5人回复
-
276求调剂。有半年电池和半年高分子实习经历
已经有6人回复
-
求调剂材料学硕080500，总分289分
已经有21人回复
-
北科281学硕材料求调剂
已经有8人回复
-
08工学调剂
已经有9人回复
-
291 求调剂
已经有4人回复
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
你好！ 我也遇到类似问题，但根据我的测试，我倾向于怀疑是lammps中special_bonds的设置问题。 默认情况下，lammps是不计算1-2,1-3,1-4原子的pairwise energy的。所以我怀疑这个设置MS和lammps不一样，但因为我查不到MS是如何设置的，所以我对比了不同的情况，最后发现，如果在lammps中不计算1-2,1-3，但计算1-4，所得到的LJ能和静电能（Ecoul+Elong）和MS中是基本一致的。 简单来说，就是在lammmps中加一句 special_bonds lj/coul 0 0 1 但这仅限于我的case，不过我建议楼主试一下。 |
16楼2015-12-16 18:20:59
lsloneil
专家顾问 (正式写手)
-
专家经验: +49 - 模拟EPI: 2
- 应助: 148 (高中生)
- 金币: 8640.9
- 散金: 507
- 红花: 30
- 帖子: 505
- 在线: 598.4小时
- 虫号: 556471
- 注册: 2008-05-11
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 分子模拟
★ ★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
月只蓝: 金币+3, 感谢指导 2015-06-02 13:05:32
小木虫: 金币+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 |
2楼2015-06-01 06:08:18
zs806750658 
铜虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 725.8
- 散金: 12
- 红花: 11
- 帖子: 462
- 在线: 497.3小时
- 虫号: 1405619
- 注册: 2011-09-17
- 性别: GG
- 专业: 矿物学(含矿物物理学）
3楼2015-06-02 04:25:18
lsloneil
专家顾问 (正式写手)
-
专家经验: +49 - 模拟EPI: 2
- 应助: 148 (高中生)
- 金币: 8640.9
- 散金: 507
- 红花: 30
- 帖子: 505
- 在线: 598.4小时
- 虫号: 556471
- 注册: 2008-05-11
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 分子模拟
4楼2015-06-02 09:34:05
