24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

astringent

铜虫 (著名写手)

应助: 16 (小学生)
金币: 2210.7
散金: 1233
红花: 4
帖子: 1180
在线: 354.8小时
虫号: 900839
注册: 2009-11-12
专业: 理论和计算化学

[交流] 【求助】相关函数归一化问题已有4人参与

最近一直在看自相关函数的程序，有个问题想不通。
自相关函数的定义是A=(a(t)*a(0))/(a(0)*a(0)),当对体系中所有粒子求和时，是先算出每个例子的(a(t)*a(0)) /(a(0)*a(0))然后求和？还是先分别求a(t)*a(0)以及a(0)*a(0)和然后再相除呢？两者有什么区别？请大家指点。

[ Last edited by ghcacj on 2010-12-8 at 09:27 ]

回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

1楼 2010-12-07 05:47:32

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

manutdhk

木虫 (小有名气)

模拟EPI: 1
应助: 5 (幼儿园)
金币: 5939.4
散金: 70
红花: 5
帖子: 94
在线: 97.3小时
虫号: 983544
注册: 2010-03-27
专业: 核技术及其应用

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

引用回帖:

644082楼: Originally posted by astringent at 2010-12-07 05:47:32:
最近一直在看自相关函数的程序，有个问题想不通。
自相关函数的定义是A=(a(t)*a(0))/(a(0)*a(0)),当对体系中所有粒子求和时，是先算出每个例子的(a(t)*a(0)) /(a(0)*a(0))然后求和？还是先分别求a(t)*a(0)以及a( ...

第二种方案才对。。。
分母只是一个归一化作用的常数，如果不计算分母，对最后的结果不会有影响的，只是Amplitude相差一个倍数。

两种方案在什么时候才会得到一致结果呢？当你所有的原子都同一的时候，也就是说，周期性边界条件下的单质材料。只要有一个原子跟其他的不一致，两种方案就不能得到一致的结果。

此外，分母部分完全没有必要放在循环里计算。
所以，我的建议，先双重循环把分子搞定，然后在另一个版块，把分母搞定，最后二者相除。

赞一下