| 查看: 4390 | 回复: 19 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
zs806750658
|
[交流]
LAMMPS输出的库伦能为什么是正值，而MS的discover算出的库伦能是负值 已有4人参与
|
|||
| 相同的模型，相同的系宗，相同的力场，相同的参数设置，为什么用Material Studio的discover计算的单点能中库伦能是-4万多，而Lammps跑完输出的能量中库伦能是正5万多,看过手册发现两个软件计算库伦能的函数是一样的，为什么相差如此之大呢。而其他的范德华能、键能等相差不是很大，哪位高手能给予指点，不胜感激 |
回复此楼» 猜你喜欢
面上合作单位盖章
已经有8人回复
-
面上函评意见出来了，像什么等级？
已经有15人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有14人回复
-
小白求助 投论文要求的highlights应该如何写
已经有4人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有12人回复
-
怎么看青基中了没有啊
已经有6人回复
-
国社科又开始会评了，不知道这次命运如何
已经有11人回复
-
有多少人是今天查系统知道结果的？
已经有19人回复
-
基金未中，这种答复是模板吗？
已经有7人回复
-
哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识
已经有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
5