| 查看: 3822 | 回复: 19 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
zs806750658
|
[交流]
LAMMPS输出的库伦能为什么是正值，而MS的discover算出的库伦能是负值 已有4人参与
|
|||
| 相同的模型，相同的系宗，相同的力场，相同的参数设置，为什么用Material Studio的discover计算的单点能中库伦能是-4万多，而Lammps跑完输出的能量中库伦能是正5万多,看过手册发现两个软件计算库伦能的函数是一样的，为什么相差如此之大呢。而其他的范德华能、键能等相差不是很大，哪位高手能给予指点，不胜感激 |
回复此楼» 猜你喜欢
全日制（定向）博士
已经有5人回复
-
假如你的研究生提出不合理要求
已经有10人回复
-
萌生出自己或许不适合搞科研的想法，现在跑or等等看？
已经有4人回复
-
Materials Today Chemistry审稿周期
已经有4人回复
-
参与限项
已经有3人回复
-
实验室接单子
已经有4人回复
-
对氯苯硼酸纯化
已经有3人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有12人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
我现在想了下，ecoul为什么是正值也许可以从以下方式理解。正如我前面所述的NaCl电荷点阵分布，因为最近的正负电荷其实在力场中已经有成键的参数bond coeffs来描述，所以这部分吸引力在计算时就不计在内，后面自然是从次近邻的电荷之间的静电力算起，次近邻是同性电荷所以结果是能量为正值。 如果把原子电荷的描述从full改成charge的话，那么最近邻的电荷之间就是异性电荷，能量为负值，我测试过个例子，把atom style从full改成charge计算出来的能量的确为负值。 但是这种理解方式却不能解释elong有时候比ecoul绝对值大，至少NaCl应该不能解释，有可能一些电荷的排布方式会造成ecoul绝对值小于elong。 |
15楼2015-11-16 08:55:01
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