[交流] 有会用CVM集团变分法计算构型熵的吗？

最近看Tamm在2015发表Acta Materialia 99（2015) 307-312上关于NiCrCo合金短程有序熵的计算，采用CVM方法计算fcc的熵。                           S=KB{(z-1)*sum(c(lnc-1)-z/2*sum((y[ij](ln(y[ij])-1)+(z/2-1)}         y[ij] 分别是： Ni–Ni        0.293333333 Ni–Cr        0.423333333 Ni–Co        0.283333333 Cr-Ni        0.423333333 Cr–Cr        0.193333333 Cr–Co        0.386666667 Co-Ni        0.283333333 Co–Cr        0.386666667 Co–Co        0.33 我理解他的公式是 ：fcc  z=12 配位数  c1=c2=c3=1/3     S=KB{(z-1)*[c1lnc1-c1+c2lnc2-c2+c3lnc3-c3]-z/2*[(y11lny11-y11+y12lny12-y12+y13lny13-y13)+(y21lny21-y21+y22lny22-y22+y23lny23-y23)+(y31lny31-y31+y32lny32-y32+y33lny33-y33)]+(z/2-1)}                                        可是为什么得不到他的结果 s=0.944    Tamm 论文地址： http://www.sciencedirect.com/science/article/pii/S1359645415005790?via%3Dihub    有懂CVM的帮帮忙啊。 S.jpg

回复此楼