24小时热门版块排行榜    

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

Waxwane

新虫 (初入文坛)

[求助] 关于Ewald sum算法的相关问题 已有1人参与

自己作死地用C编写Monte Carlo的代码,模拟的具体问题是粘土矿物对Ar的吸附。
然后需要加上Ewald Sum算法的问题算矿物羟基的长程力。
但是这个矿物是三斜晶系的,书上给的Ewald Sum是立方格子。
我查过LAMMPS里面,前几年有个Ewald/n方法专门算三斜格子,现在和Ewald本身合并了。
我决定用三斜的格子算,因为矿物的周期性是三斜的,用立方重复不了。但是实验室其他人不搞算法,他们说关于三斜粘土矿物的文献,没有听过把格子划分为三斜格子的说法。
我认为算短程力的格子可以是立方的,但长程力不行,那些文献里没有写清楚。
不知道我的观点正不正确,因为我只能找到立方格子的源码和伪代码,没见过三斜的代码。
还有,如果有编程简单一点的,被大家比较公认的算长程力的算法,求推荐。我已经要被Ewald sum坑死了。
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

DDXZS

铜虫 (小有名气)

您好,我是做粘土矿物吸附的,刚刚开始接触分子动力学,想向您请教一下有什么推荐的书籍或者好的学习方法,希望能得到您的帮助,万分感谢您!
6楼2020-05-25 20:55:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Waxwane 的主题更新
信息提示
请填处理意见