24小时热门版块排行榜    

查看: 993  |  回复: 4
当前主题已经存档。

tang327

新虫 (初入文坛)

[交流] 【求助】怎样计算长程静电作用项

我要计算的力场除了有个LJ相互作用还有长程静电相互作用,请问怎样计算呀?算例中的field文件只有LJ相互作用,不知道怎样把长程静电相互作用加进去,哪位高手会的话,请赐教,万分感谢!!!
回复此楼

» 猜你喜欢

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

tianlangxingaa

铁杆木虫 (著名写手)

★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
lei0736(金币+4,VIP+0):谢谢 1-14 20:16
计算长程相互作用的方法有很多的。

第1个就是1921年ewald提出的ewald summation 方法。 该方法是最精确的方法,但是计算比较慢,scaling 是O(N**2),如果选择参数比较好可以降低到O(N**(3/2))

然后是reaction field method(不过该方法对于homogeneous还可以,对于多组分体系,是不太适合的)

然后就是对ewald方法进行改进:采用快速傅立叶变换技术,使scaling降低到 O(NlogN)
1)particle-particle particle-mesh method
2)Particle mesh ewald methos
3)ENUF method = ewald + NFFT +FFTW (这个是我目前采用的,是我老板开发的)

然后就是multipole based method方法 (这个不太熟悉,不过你可以看看)

1999年wolf提出了一个新的方法,使体系的计算scaliing从 O(NlogN) 降低到了 O(N)

这些就是我所知道的方法。希望对你有帮助

cheers
2楼2010-01-14 16:55:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tang327

新虫 (初入文坛)

谢谢版主了
3楼2010-01-15 09:48:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beefly

专家顾问 (职业作家)

地沟油冶炼专家

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
lei0736(金币+2,VIP+0):谢谢 1-15 21:00
multipole方法:Gaussian的PBC模块使用fast multipole method,O(N);日本有人开发了multipole+有限元,速度比前者快一两个数量级。

wolf的方法已经加入到gulp程序里。和Ewald方法不同,这个方法适用于0D-3D体系。
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
4楼2010-01-15 13:00:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

老虎大王

木虫 (著名写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
lei0736(金币+1,VIP+0):谢谢 1-15 21:00
Ewald方法是最为经典的。
5楼2010-01-15 16:28:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 tang327 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见