单个水分子势能计算问题

gongrehk

sobereva 发表于 2021-7-19 21:52
是

好的 谢谢老师 非常感谢！

gongrehk

k64_cc 发表于 2021-7-15 20:53
你把体系分成两部分A和B，总势能就可以被分解为U(A), U(B) 和 U(AB)。相互作用能就是相互作用能，不要创 ...

前辈能详细说一下为什么需要改用Cut-off吗？最近碰到了这个问题。
PME的短程是原电荷(ri)和screening charge产生的电势是吗？这部分衰减很快，可以视为短程部分。
倒易空间内的是cancelling charge产生的，是长程部分。
但是这个过程也是计算了原电荷(ri)处的静电势能啊，对于一个分子来说加和这部分不就是该分子的静电势能了吗？为什么能量分解处会有困扰呢?真的很好奇这个问题。

k64_cc

gongrehk 发表于 2021-8-23 18:30
前辈能详细说一下为什么需要改用Cut-off吗？最近碰到了这个问题。
PME的短程是原电荷(ri)和screening ch ...

Ewald sum可以理解为
point charge - gaussian charge + gaussian charge ，
其中前两项加一起是Ewald sum的短程部分，最后一项是长程部分。能量分解要的是point charge的贡献，不需要 - gaussian charge部分，所以Ewald sum短程部分直接拿来用不很恰当。如果只想要point charge部分，那最简单的方法显然是rerun了。

gongrehk

k64_cc 发表于 2021-8-23 19:39
Ewald sum可以理解为
point charge - gaussian charge + gaussian charge ，
其中前两项加一起是Ewald  ...

- gaussian charge + gaussian charge的结果不就相抵消了吗，那这样的话短程部分加上长程部分的结果不还是point charge的贡献吗?

k64_cc

gongrehk 发表于 2021-8-24 17:45
- gaussian charge + gaussian charge的结果不就相抵消了吗，那这样的话短程部分加上长程部分的结果不还 ...

你要是能把长程部分给分解了，可以考虑拿出来单发一个JCP。

gongrehk

k64_cc 发表于 2021-8-25 14:03
你要是能把长程部分给分解了，可以考虑拿出来单发一个JCP。

没明白 不过还是谢谢您  这个问题好像确实也没什么意义，我先跳过这个吧