24小时热门版块排行榜    

查看: 1994  |  回复: 13
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

chunwang

铜虫 (小有名气)

[交流] 【求助】【求助】配位数的计算

用分子动力学模拟得到了径向分布函数,但不知道在程序里怎样实现配位数的计算。计算公式知道,就是不知道怎么操作,还请各位高手指点指点!

[ Last edited by lei0736 on 2009-11-24 at 19:58 ]
回复此楼

» 收录本帖的淘帖专辑推荐

分子模拟 硫系玻璃

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

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

老虎大王

木虫 (著名写手)

★ ★ ★
lei0736(金币+3,VIP+0):谢谢 4-27 08:58
你算不出配位数,那是怎么算出RDF的呢?用MD算RDF,第一步是找出不同距离内的配位数,第二步是归一化得到RDF。

如果你只有一个RDF曲线,然后再通过积分来反算配位数,可能会有一些计算误差的。
9楼2009-04-24 21:12:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

老虎大王

木虫 (著名写手)

★ ★ ★ ★ ★
lei0736(金币+5,VIP+0):谢谢 4-27 11:30
或者,如果你有构型,并且用现成的程序输出了RDF,你不想修改人家的程序来输出配位数,那么建议你用搜索构型的方法求配位数,而不必根据RDF去积分求配位数(你不清楚人家具体的归一化方法,积分容易出错)。求配位数的方法是:
1。根据RDF或其它知识,选择一个配位距离(比如RDF第一谷的位置,设为r0)
2. 将任一原子设为中心原子,搜索与其距离在r0以内的原子个数(注意周期性边界条件,对多组分体系可以分类)
3. 将第2步遍及所有中心原子
4. 将第2和第3步遍及所有不同时刻的构型。
5. 求平均。(包括中心原子平均和构型(即时间)平均)。
这就是配位数的算法,简单得很。若要求配位数随配位距离的变化,再选不同的距离,多一层循环就是了。

[ Last edited by 老虎大王 on 2009-4-27 at 10:24 ]
11楼2009-04-27 10:22:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 chunwang 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见