| 查看: 4063 | 回复: 19 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
zs806750658
|
[交流]
LAMMPS输出的库伦能为什么是正值，而MS的discover算出的库伦能是负值 已有4人参与
|
|||
| 相同的模型，相同的系宗，相同的力场，相同的参数设置，为什么用Material Studio的discover计算的单点能中库伦能是-4万多，而Lammps跑完输出的能量中库伦能是正5万多,看过手册发现两个软件计算库伦能的函数是一样的，为什么相差如此之大呢。而其他的范德华能、键能等相差不是很大，哪位高手能给予指点，不胜感激 |
回复此楼» 猜你喜欢
人工智能320调剂08工类还有机会吗
已经有17人回复
-
考研英一数一338分
已经有10人回复
-
求助调剂，跨调
已经有15人回复
-
085600材料与化工329分求调剂
已经有20人回复
-
085600材料与化工349分求调剂
已经有15人回复
-
求调剂
已经有13人回复
-
一志愿华南理工大学331分材料求调剂
已经有11人回复
-
271求调剂
已经有40人回复
-
求调剂学校
已经有3人回复
-
332求调剂
已经有15人回复
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
谢谢，你所说的轨迹文件是dump出来的吗，轨迹文件是minimization还是MD计算出来的？轨迹文件如果可以用MS读出charge的话，应该在文件里就可以找到相应的charge值，假设轨迹文件是文本文件的话。 MS里体系原子电荷变化有几个可能： 1. LAMMPS计算过程中电荷变化了，正到负（不大可能），因为电荷没有平衡所以计算过程电荷应该没变。 2. LAMMPS计算给出的电荷在MS读入时所有电荷都变正了，MS读入是忽略了正负号。 3. LAMMPS计算给出的电荷没有读入，MS重新赋予原子电荷（不大可能），因为MS赋予电荷不会全部都是正值。 所以最有可能的就是2. 这个有可能从轨迹文件中得到验证，如果轨迹文件是文本就能找出LAMMPS输出的电荷。 现在我对这个short-range的ecoul一般为正值很不理解，想象下一个向NaCl的电荷点阵，正负电荷是距离最近的，总的能量应该是负值才对，不知你对这个有什么好的理解。谢谢。 |
13楼2015-11-12 13:51:09
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
10