24小时热门版块排行榜

返回列表

混乱的流动

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 306.6
散金: 7
帖子: 62
在线: 57.4小时
虫号: 2497733
注册: 2013-06-05
性别: GG
专业: 多相流热物理学

[求助] Ewald summation求助已有1人参与

请问谁有Ewald summation 的程序段，我想用一下。实现以下目的：在三维计算区域内有1000个球体，带电量为q，其中y-z方向是周期性边界条件，x方向没有周期性。用Ewald 应该能算出库伦静电作用势或者静电作用力吧？

回复此楼

» 猜你喜欢

申博求助已经有4人回复
27年申博自荐-能源/环境方向已经有3人回复
闲聊已经有3人回复
亲密关系已经有13人回复
27秋申博自荐已经有4人回复

1楼 2015-11-09 14:44:16

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

lsloneil

专家顾问 (正式写手)

专家经验: +49
模拟EPI: 2
应助: 148 (高中生)
金币: 8640.9
散金: 507
红花: 30
帖子: 505
在线: 598.4小时
虫号: 556471
注册: 2008-05-11
性别: GG
专业: 理论和计算化学
管辖: 分子模拟

【答案】应助回帖

感谢参与，应助指数 +1

如果使用lammps，可以使用kspace_modify slab这个命令，详细说明参见
http://lammps.sandia.gov/doc/kspace_modify.html

不过说明一下，Ewald summation适用的是电中性体系，如果你的体系非中性，贸然使用Ewald summation会带来误差，而且体系总电荷距离0越远，误差越大。

赞一下

回复此楼

2楼2015-11-09 16:09:41

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

混乱的流动

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 306.6
散金: 7
帖子: 62
在线: 57.4小时
虫号: 2497733
注册: 2013-06-05
性别: GG
专业: 多相流热物理学

引用回帖:

2楼: Originally posted by lsloneil at 2015-11-09 16:09:41
如果使用lammps，可以使用kspace_modify slab这个命令，详细说明参见
http://lammps.sandia.gov/doc/kspace_modify.html

不过说明一下，Ewald summation适用的是电中性体系，如果你的体系非中性，贸然使用Ewald ...

也就是正负电荷量差不多才行？为什么有这样的特性？

发自小木虫IOS客户端

赞一下

回复此楼

3楼2015-11-09 20:02:28

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

lsloneil

专家顾问 (正式写手)

专家经验: +49
模拟EPI: 2
应助: 148 (高中生)
金币: 8640.9
散金: 507
红花: 30
帖子: 505
在线: 598.4小时
虫号: 556471
注册: 2008-05-11
性别: GG
专业: 理论和计算化学
管辖: 分子模拟

引用回帖:

3楼: Originally posted by 混乱的流动 at 2015-11-09 00:02:28
也就是正负电荷量差不多才行？为什么有这样的特性？
...

是的。因为如果总电荷不为0的话，Ewald summation在k空间求和那一部分的级数会发散。通常为了避免这种问题，分子模拟的程序在处理非电中性系统的Ewald summation是都会给系统加上一个均匀分布的背景电荷密度来中和这部分电荷。如果你的体系净电荷不大（比如说只有-2 e），那么这个背景电荷的影响是比较小的，反之，如果你的体系净电荷很大，那么Ewald summation带来的误差就会很大。

赞一下

回复此楼

4楼2015-11-09 23:24:06

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖