24小时热门版块排行榜    

查看: 2035  |  回复: 10
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

chenjingzhi

金虫 (小有名气)

[求助] 请教个长程力计算的问题 已有2人参与

在Ewald处理长程力中,长程力被分割成三部分,一部分为收敛较快的短程作用,这部分直接在实空间中采取截断进行计算,第二部分为收敛较慢的作用,采用傅里叶变换在倒易空间中计算,这两部分的最优计算时间均与N的二分之三次方成正比,第三部分则为自身电荷的更正,也是直接计算。与Ewald相比,粒子格点Ewald(PME),按我的理解,似乎只是对Ewald计算中的收敛较慢的部分采用快速傅里叶变换进行处理,使这部分的计算时间缩短为NlogN,但是似乎对收敛较快的那部分没有进行改动,计算时间依然为N的二分之三次方,那这样似乎总的计算时间还是和N的二分之三次方成正比,不知道我的理解错在哪里了呢,希望高手指教。
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

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

赵红霞

铁杆木虫 (著名写手)

引用回帖:
9楼: Originally posted by chenjingzhi at 2014-08-12 11:35:00
傅里叶变换没什么啊,可以参考Daan Frenkel的Understanding molecular simulation,里面对Ewald的讲解已经很详细了,也可以找一下Computer simulation of liquids里面的代码来看看,不过是计算静电势能的,如果还需 ...

请问你是哪个学校或者单位的?
10楼2014-08-12 14:35:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 11 个回答

lsloneil

专家顾问 (正式写手)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
chenjingzhi: 金币+2, 有帮助 2014-08-11 12:40:24
收敛较快的短程作用部分计算复杂度为O(N),也就是说与N呈线性关系,这比k空间部分的计算复杂度(O(NlogN) 或者 O(N^3/2))要快很多,所以说k空间部分的计算速度是起决定性作用的。
2楼2014-08-10 10:24:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fhh2626

木虫 (正式写手)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
chenjingzhi: 金币+2, 有帮助 2014-08-11 12:40:35
因为收敛较快的部分是截断半径计算的,这个半径远小于盒子大小,所以对于每个原子,截断半径里面的原子数都不会随N变化,显然可以看出这部分时间复杂度是O(n)而不是O(n^3/2)

[ 发自手机版 http://muchong.com/3g ]
3楼2014-08-10 11:37:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chenjingzhi

金虫 (小有名气)

但是实空间部分的阶段半径不也依赖高斯分布宽度吗,这个宽度似乎随粒子数变化的
4楼2014-08-11 12:39:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 2JOx3r2CYEgw 2026-08-22 4/200 2026-08-23 19:22 by lXgU0RvNTDC9
[论文投稿] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-21 8/400 2026-08-23 18:58 by lXgU0RvNTDC9
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-21 5/250 2026-08-23 18:55 by lXgU0RvNTDC9
[论文投稿] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 DpSrDtM079iu 2026-08-22 3/150 2026-08-23 15:33 by KM5EcsNQRBPn
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 2JOx3r2CYEgw 2026-08-21 5/250 2026-08-23 14:33 by KM5EcsNQRBPn
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-21 6/300 2026-08-23 14:21 by KM5EcsNQRBPn
[基金申请] 2026国自然函评费到账 +15 羊腰板 2026-08-21 16/800 2026-08-23 10:45 by process2012
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-21 9/450 2026-08-23 10:44 by IXZuIJ2Q7OVy
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 2JOx3r2CYEgw 2026-08-22 6/300 2026-08-23 03:40 by OEbVnUOu01ol
[基金申请] 范进中举一文的中心思想 +3 zjuhero 2026-08-22 4/200 2026-08-22 20:02 by mycsru
[基金申请] 今天基金会出结果吗?20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 科研孤儿太难了 +17 我4大白菜 2026-08-20 18/900 2026-08-21 20:57 by zhangev
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[论文投稿] 投稿咨询 +5 wwm09 2026-08-17 7/350 2026-08-21 10:11 by 期刊论文帮手
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 今天系统多次维护,明天很可能放榜! +10 zju2000 2026-08-16 11/550 2026-08-20 20:02 by cl479861084
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
请填处理意见