24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 1662  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料专硕322 +6 哈哈哈吼吼吼哈 2026-04-05 6/300 2026-04-06 00:19 by T可可西里T
[考研] 336材料与化工085600求调剂 +6 水星记infp 2026-04-05 6/300 2026-04-05 23:18 by 来看流星雨10
[考研] 294求调剂 +4 Grey_Ey 2026-04-01 5/250 2026-04-05 23:05 by Grey_Ey
[考研] 化学调剂 +15 艾志恒 2026-04-03 16/800 2026-04-05 22:36 by syh9288
[考研] 材料专硕322分 +10 哈哈哈吼吼吼哈 2026-04-04 10/500 2026-04-05 21:22 by 学员8dgXkO
[考研] 求助071001调剂!!! +3 黄守松 2026-04-05 4/200 2026-04-05 21:01 by barlinike
[考研] 296求调剂 +3 汪!?! 2026-04-05 4/200 2026-04-05 20:13 by 啵啵啵0119
[考研] 求调剂 一志愿西南交通大学085701环境工程 282分 +7 多多爱吃汉堡 2026-04-04 7/350 2026-04-05 19:47 by 福农资环_环境基
[考研] 0703化学321分求调剂 +17 三dd. 2026-03-30 18/900 2026-04-05 18:07 by 蓝云思雨
[考研] 一志愿电子科技大学085600材料与化工 329分求调剂 +10 Naiko 2026-04-04 10/500 2026-04-05 09:40 by sam3303
[考研] 353求调剂 +10 MayUxw1 2026-04-03 10/500 2026-04-05 09:23 by 无际的草原
[考研] 344材料与化工调剂 +9 调剂上岸玘 2026-04-03 9/450 2026-04-04 23:10 by happyddm
[考研] 一志愿南农090401,268,求调剂 +5 一木鸟然 2026-04-04 5/250 2026-04-04 17:07 by babysonlkd
[考研] 265求调剂 +17 林深温澜 2026-04-01 20/1000 2026-04-04 01:09 by userper
[考研] 本科985,专业0812分336求调剂 +4 莫莫很行 2026-04-03 4/200 2026-04-03 21:31 by zhq0425
[考研] 一志愿中国科学院大学265求调剂 +9 恬淡ye 2026-03-31 10/500 2026-04-03 11:10 by txp1986
[考研] 309求调剂 +14 呆菇不是戴夫 2026-04-02 14/700 2026-04-03 09:42 by 蓝云思雨
[考研] 326求调剂 +3 9ahye 2026-04-02 4/200 2026-04-03 08:43 by Jaylen.
[考研] 085602化工求调剂(331分) +9 111@127 2026-03-30 9/450 2026-04-02 20:00 by dick_runner
[考研] 11408 321分求调剂 +3 huchun12138 2026-03-30 4/200 2026-04-01 22:48 by guanxin1001
信息提示
请填处理意见