| 查看: 2021 | 回复: 5 | ||||
混乱的流动银虫 (小有名气)
|
[求助]
Ewald summation求助 已有1人参与
|
| 请问谁有Ewald summation 的程序段，我想用一下。实现以下目的：在三维计算区域内有1000个球体，带电量为q，其中y-z方向是周期性边界条件，x方向没有周期性。用Ewald 应该能算出库伦静电作用势或者静电作用力吧？ |
回复此楼» 猜你喜欢
最失望的一年
已经有12人回复
-
拟解决的关键科学问题还要不要写
已经有8人回复
-
存款400万可以在学校里躺平吗
已经有29人回复
-
求推荐英文EI期刊
已经有5人回复
-
请教限项目规定
已经有4人回复
-
国自然申请面上模板最新2026版出了吗？
已经有20人回复
-
26申博
已经有3人回复
-
基金委咋了？2026年的指南还没有出来？
已经有10人回复
-
基金申报
已经有6人回复
-
疑惑？
已经有5人回复
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