各位大神，請(qǐng)教一個(gè)問題。
本人按照http://www.gaoyang168.com/t-3181757-1帖子中valenhou001提供的方法擬合線膨脹系數(shù)。

具體步驟為：
大致可采用這樣的步驟：
1）在實(shí)驗(yàn)或理論計(jì)算出來(lái)的晶格常數(shù)或平衡體積(v0)的左右各取6～7體積點(diǎn)，使得所要計(jì)算的體積點(diǎn)的個(gè)數(shù)在10個(gè)左右，而且這些體積的范圍大致在[(1- x) v0,  (1+x) v0]，x為小于0.1~0.3的正數(shù)。

2）對(duì)這10來(lái)個(gè)體積點(diǎn)分別計(jì)算總能e以及聲子譜，然后采用準(zhǔn)諧近似計(jì)算出每個(gè)體積下在不同溫度下的自由能f(v_i, t);

3）然后將數(shù)據(jù)的排列進(jìn)行調(diào)整為：每個(gè)溫度t_i下，自由能排列成
v_1          f(v_1, t_i)
v_2          f(v_2, t_i)
....
v_10         f(v_10, t_i)
然后對(duì)這些10對(duì)v～f數(shù)據(jù)點(diǎn)用狀態(tài)方程擬合得到平衡體積v0(t_i)。有了不同溫度下的平衡體積,然后利用這些數(shù)據(jù)點(diǎn)t_i~v0(t_i)進(jìn)行簡(jiǎn)單的差分求導(dǎo)就可以得到熱膨脹系數(shù)。

但是當(dāng)進(jìn)行到第三步時(shí)，用狀態(tài)方程擬合各溫度下的平衡體積時(shí)，發(fā)現(xiàn)隨著晶格參數(shù)的增加，結(jié)構(gòu)的自由能呈線性變化，逐漸降低，并不能擬合出各溫度下的平衡體積，請(qǐng)問這是怎么回事？ 返回小木蟲查看更多