| 查看: 2093 | 回复: 5 | ||||
混乱的流动银虫 (小有名气)
|
[求助]
Ewald summation求助 已有1人参与
|
| 请问谁有Ewald summation 的程序段，我想用一下。实现以下目的：在三维计算区域内有1000个球体，带电量为q，其中y-z方向是周期性边界条件，x方向没有周期性。用Ewald 应该能算出库伦静电作用势或者静电作用力吧？ |
回复此楼» 猜你喜欢
招收生物学/细胞生物学调剂
已经有4人回复
-
一志愿武汉理工，总分321，英一数二，求老师收留。
已经有5人回复
-
324求调剂
已经有8人回复
-
284求调剂
已经有12人回复
-
348求调剂
已经有3人回复
-
329求调剂
已经有7人回复
-
求调剂，一志愿 南京航空航天大学大学 ，080500材料科学与工程学硕
已经有4人回复
-
考研调剂
已经有9人回复
-
304材料求调剂
已经有4人回复
-
一志愿郑大085600，310分求调剂
已经有3人回复
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
混乱的流动
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 306.6
- 散金: 7
- 帖子: 62
- 在线: 57.4小时
- 虫号: 2497733
- 注册: 2013-06-05
- 性别: GG
- 专业: 多相流热物理学
3楼2015-11-09 20:02:28
lsloneil
专家顾问 (正式写手)
-
专家经验: +49 - 模拟EPI: 2
- 应助: 148 (高中生)
- 金币: 8640.9
- 散金: 507
- 红花: 30
- 帖子: 505
- 在线: 598.4小时
- 虫号: 556471
- 注册: 2008-05-11
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 分子模拟
4楼2015-11-09 23:24:06
chuntao118
新虫 (著名写手)
- 应助: 93 (初中生)
- 金币: 2590.9
- 散金: 1402
- 红花: 25
- 帖子: 1611
- 在线: 241.2小时
- 虫号: 1791567
- 注册: 2012-05-02
- 专业: 计算数学与科学工程计算
5楼2015-11-10 08:26:32
lsloneil
专家顾问 (正式写手)
-
专家经验: +49 - 模拟EPI: 2
- 应助: 148 (高中生)
- 金币: 8640.9
- 散金: 507
- 红花: 30
- 帖子: 505
- 在线: 598.4小时
- 虫号: 556471
- 注册: 2008-05-11
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 分子模拟
6楼2015-11-10 09:26:52
