| 查看: 1038 | 回复: 9 | |||
| 当前主题已经存档。 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
liuyalihd铜虫 (小有名气)
|
[交流]
【求助】Ewald长程力中Fourier空间形状问题?
|
||
|
我看某文献1中计算Ewald长程力时,分成实空间和Fourier空间. 实空间中其实还是在截断半径内计算两个粒子之间的作用. Fourier空间是在半球形空间(z轴0到nz, y轴从-ny到ny, x轴从-nx到nx)内计算. 但在gromacs4.0.5中,我看了Ewald.c的代码,发现它没有分成实空间和Fourier空间两部分计算.它似乎仅在Fourier空间中计算. 最让我奇怪的是它的Fourier空间不是半球形,而是在半正方体(z轴0到nz, y轴从-ny到ny, x轴从-nx到nx)中计算. 请问为何gromacs4.0.5中计算Ewald的Fourier空间形状不是传统的球形呢? gromacs4.0.5中计算Ewald是咋回事?为何没有实空间的计算呢? [ Last edited by zdhlover on 2009-11-15 at 02:29 ] |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有9人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有5人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有7人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有6人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有6人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有5人回复
liuyalihd
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 117.7
- 散金: 577
- 帖子: 59
- 在线: 8.7小时
- 虫号: 531310
- 注册: 2008-03-23
- 专业: 计算机体系结构
|
致 zixi: 关于您的第一条：非常感谢，我再找找real空间的计算。 关于您的第二条： 1）原理部分： 前面是我没有讲清楚。球形空间是Ewald求和法的周期作用区域，这在很多书中都提到了（如：分子模拟-陈正隆，9.2节 图9-5）。 2）代码部分：另外，我给出的代码不是gromacs的。但是我看在gromacs的ewald.c（src/mdlib/ewald.c）中有函数do_ewald()，我看代码非常像是计算Ewald中的Fourier计算。但是它没有if判断镜像是否在球形空间中，我认为应该是球形才计算，在正方体空间计算不符合原理呀！？ 3）正方体空间与球形空间相比，去掉if语句的开销不提，球形空间减少了一半的计算量。 |
7楼2009-09-04 12:37:02
bay__gulf
金虫 (著名写手)
刘苏州
- 模拟EPI: 8
- 应助: 9 (幼儿园)
- 贵宾: 4.85
- 金币: 2332.8
- 红花: 1
- 帖子: 1344
- 在线: 271小时
- 虫号: 592012
- 注册: 2008-09-03
- 专业: 理论和计算化学
- 管辖: 分子模拟
2楼2009-09-01 20:52:07
liuyalihd
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 117.7
- 散金: 577
- 帖子: 59
- 在线: 8.7小时
- 虫号: 531310
- 注册: 2008-03-23
- 专业: 计算机体系结构
3楼2009-09-01 22:22:09
yongleli
木虫 (正式写手)
- 模拟EPI: 1
- 应助: 93 (初中生)
- 金币: 3769.1
- 散金: 369
- 红花: 24
- 帖子: 764
- 在线: 101.3小时
- 虫号: 303595
- 注册: 2006-12-03
- 专业: 原子和分子物理
4楼2009-09-02 15:01:01
