24小时热门版块排行榜    

查看: 2024  |  回复: 10

chenjingzhi

金虫 (小有名气)

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

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

» 猜你喜欢

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

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

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的回帖

赵红霞

铁杆木虫 (著名写手)

请教楼主,这个长程作用力是自己编写还是采用程序包里的?
最近遇到同样的问题,不知道如何处理
5楼2014-08-11 14:46:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chenjingzhi

金虫 (小有名气)

引用回帖:
5楼: Originally posted by 赵红霞 at 2014-08-11 14:46:11
请教楼主,这个长程作用力是自己编写还是采用程序包里的?
最近遇到同样的问题,不知道如何处理

自己写的
6楼2014-08-11 16:00:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

赵红霞

铁杆木虫 (著名写手)

引用回帖:
6楼: Originally posted by chenjingzhi at 2014-08-11 16:00:09
自己写的...

参考出处呢?我不会傅里叶变换,请楼主赐教
7楼2014-08-12 09:37:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lsloneil

专家顾问 (正式写手)

【答案】应助回帖

引用回帖:
4楼: Originally posted by chenjingzhi at 2014-08-10 16:39:06
但是实空间部分的阶段半径不也依赖高斯分布宽度吗,这个宽度似乎随粒子数变化的

只要使用cutoff的就是短程相互作用,对于短程相互作用计算复杂度就是O(N)。
8楼2014-08-12 10:55:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chenjingzhi

金虫 (小有名气)

引用回帖:
7楼: Originally posted by 赵红霞 at 2014-08-12 09:37:08
参考出处呢?我不会傅里叶变换,请楼主赐教...

傅里叶变换没什么啊,可以参考Daan Frenkel的Understanding molecular simulation,里面对Ewald的讲解已经很详细了,也可以找一下Computer simulation of liquids里面的代码来看看,不过是计算静电势能的,如果还需要什么资料的话,可以再联系
9楼2014-08-12 11:35:00
已阅   回复此楼   关注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的回帖
相关版块跳转 我要订阅楼主 chenjingzhi 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 2026国自然函评费到账 +13 羊腰板 2026-08-21 13/650 2026-08-23 09:25 by lfy8008
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 DpSrDtM079iu 2026-08-22 5/250 2026-08-23 04:07 by OEbVnUOu01ol
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 2JOx3r2CYEgw 2026-08-22 7/350 2026-08-23 03:55 by OEbVnUOu01ol
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 QTy3jDtz1uLt 2026-08-21 9/450 2026-08-23 03:04 by OEbVnUOu01ol
[基金申请] 93BebMhtakh前后11位开头都是大写 +7 且听虎啸 2026-08-17 8/400 2026-08-22 21:55 by 医学老男孩
[基金申请] 范进中举一文的中心思想 +3 zjuhero 2026-08-22 4/200 2026-08-22 20:02 by mycsru
[基金申请] 什么时候开奖? +7 CrisMessi 2026-08-18 8/400 2026-08-22 18:29 by 淀粉搬运工
[基金申请] filecode,4个jtjc了 +13 ziyangfang 2026-08-19 16/800 2026-08-22 17:08 by WH3796
[基金申请] 人气不行了 +8 fansofjerry 2026-08-21 8/400 2026-08-22 16:30 by zyqchem
[基金申请] 放榜前的不淡定 40+4 snowwithsea 2026-08-19 14/700 2026-08-21 23:51 by cratir
[基金申请] 建议基金发布提前给出明确的时间点 +10 kulium 2026-08-21 13/650 2026-08-21 21:58 by alongwaytogo
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +16 医学老男孩 2026-08-20 20/1000 2026-08-21 21:13 by Ldrop2023
[基金申请] 让我中一个面上吧! +11 大萍1987 2026-08-20 13/650 2026-08-21 14:48 by 20120902066
[基金申请] 感觉是下周放榜了 +7 angus9576 2026-08-17 12/600 2026-08-21 13:38 by weiyin
[基金申请] 今天放榜没戏了吧 +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
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
信息提示
请填处理意见