| 查看: 2082 | 回复: 5 | |||||
混乱的流动银虫 (小有名气)
|
[求助]
Ewald summation求助 已有1人参与
|
| 请问谁有Ewald summation 的程序段，我想用一下。实现以下目的：在三维计算区域内有1000个球体，带电量为q，其中y-z方向是周期性边界条件，x方向没有周期性。用Ewald 应该能算出库伦静电作用势或者静电作用力吧？ |
回复此楼» 猜你喜欢
323求调剂
已经有6人回复
-
一志愿北京化工大学 070300 学硕 336分 求调剂
已经有4人回复
-
352求调剂
已经有3人回复
-
一志愿东华大学化学070300，求调剂
已经有8人回复
-
277材料科学与工程080500求调剂
已经有7人回复
-
317求调剂
已经有18人回复
-
293求调剂
已经有5人回复
-
280分求调剂 一志愿085802
已经有7人回复
-
0854电子信息求调剂
已经有3人回复
-
263求调剂
已经有4人回复
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
