| 查看: 4063 | 回复: 19 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
zs806750658
|
[交流]
LAMMPS输出的库伦能为什么是正值，而MS的discover算出的库伦能是负值 已有4人参与
|
|||
| 相同的模型，相同的系宗，相同的力场，相同的参数设置，为什么用Material Studio的discover计算的单点能中库伦能是-4万多，而Lammps跑完输出的能量中库伦能是正5万多,看过手册发现两个软件计算库伦能的函数是一样的，为什么相差如此之大呢。而其他的范德华能、键能等相差不是很大，哪位高手能给予指点，不胜感激 |
回复此楼» 猜你喜欢
085600材料与化工329分求调剂
已经有20人回复
-
085600材料与化工349分求调剂
已经有15人回复
-
求调剂
已经有13人回复
-
一志愿华南理工大学331分材料求调剂
已经有11人回复
-
271求调剂
已经有40人回复
-
求调剂学校
已经有3人回复
-
332求调剂
已经有15人回复
-
297求调剂
已经有25人回复
-
本科西工大 324求调剂
已经有5人回复
-
297工科调剂?
已经有7人回复
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
追问下，前面帖子里你说计算结果没有问题是基于什么标准呢？因为我做过一个优化测试，LAMMPS有没有静电力优化出来的结构是不一样的，LAMMPS没有静电力优化出来的结果和MS的结果差不多。就是LAMMPS用力场分配电荷计算时优化出来的结构和MS是不一样的。这样来说LAMMPS和MS两个必然有个优化出来的结构是错的，我倾向认为LAMMPS优化出来的结构不合理。 还有，问下关于电荷的问题，我对LAMMPS对电荷处理不是很清楚，可能问题比较浅显，望见谅。电荷在计算过程中是不变的吧，是不是要设置计算过程电荷重新平衡，电荷才会变化？最后优化结构的电荷我在log.lammps里没有找到，需要在in文件加入参数才能输出吗？谢谢。 |
11楼2015-11-12 09:17:26
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
