24小时热门版块排行榜    

查看: 1628  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 327求调剂 +6 拾光任染 2026-03-15 11/550 2026-03-15 22:47 by 拾光任染
[考研] 085601材料工程315分求调剂 +3 yang_0104 2026-03-15 3/150 2026-03-15 10:58 by peike
[考研] 297一志愿上交085600求调剂 +5 指尖八千里 2026-03-14 5/250 2026-03-14 17:26 by a不易
[考研] 【0703化学调剂】-一志愿华中师范大学-六级475 +5 Becho359 2026-03-11 5/250 2026-03-14 11:35 by 哦哦123
[考研] 学硕285求调剂 +13 Wisjxn 2026-03-12 46/2300 2026-03-14 10:33 by JourneyLucky
[考研] 一志愿天大化工(085600)调剂总分338 +6 蔡大美女 2026-03-09 6/300 2026-03-14 02:46 by JourneyLucky
[考研] 材料与化工 一志愿山大 321分 求调剂 +7 每天散步 2026-03-09 8/400 2026-03-14 02:18 by JourneyLucky
[考研] 一志愿郑大070303,338分,求调剂 +4 dadawaf 2026-03-10 5/250 2026-03-14 01:20 by lsw010101
[考研] 271求调剂 +10 生如夏花… 2026-03-11 10/500 2026-03-14 00:35 by 卖报员小雨
[考研] 308求调剂 +3 是Lupa啊 2026-03-10 3/150 2026-03-14 00:30 by JourneyLucky
[考研] 327求调剂 +4 Ffff03 2026-03-10 4/200 2026-03-14 00:17 by JourneyLucky
[考研] 材料371求调剂 +9 鳄鱼? 2026-03-11 11/550 2026-03-13 22:53 by JourneyLucky
[考研] 26调剂/材料/英一数二/总分289/已过A区线 +6 步川酷紫123 2026-03-13 6/300 2026-03-13 21:59 by 星空星月
[考研] 0703化学一志愿211 总分320求调剂 +5 玛卡巴卡啊哈 2026-03-11 5/250 2026-03-13 21:40 by JourneyLucky
[考研] (081700)化学工程与技术-298分求调剂 +12 11啦啦啦 2026-03-11 35/1750 2026-03-13 21:25 by JourneyLucky
[考研] 材料工程调剂 +4 咪咪空空 2026-03-11 4/200 2026-03-13 19:57 by JourneyLucky
[考研] 085600材料与化工 309分请求调剂 +7 dtdxzxx 2026-03-12 8/400 2026-03-13 14:43 by jxchenghu
[考研] 277求调剂 +4 anchor17 2026-03-12 4/200 2026-03-13 11:15 by 白夜悠长
[考研] 材料调剂,307分 +13 张泳铭1 2026-03-09 17/850 2026-03-13 11:09 by 薛云鹏
[考研] 纺织、生物、化学、材料相关专业招生了 +4 耶耶业 2026-03-09 7/350 2026-03-12 19:05 by Equinoxhua
信息提示
请填处理意见