一种近似是认为二者一样.反正是在常温常压下,距离熔点很远. 一般来说,VASP计算不能加上温度,除非是做分子动力学的计算,来模拟固定在一定温度下的稳定性.要计算高温下固体的能量和晶格常数、体积,可以计算出不同晶格常数(也就是体积下)固体的声子振动频率,再得到自由能,利用F(V)的关系得到在某一个温度下的 平衡体积.(这种办法是一种准谐近似),PRB上面有很多类似的文章. 有人会问: (1)既然涉及到温度都只能是MD了,为何还能计算高温下的声子振动频率? 答:这是因为MD计算也能求解出声子的振动频率,这个在Rev.Mod.Phys.2001的一篇文章就有提到. (2)得到声子的振动频率在VASP中是设定IBRION=5,但是只能计算Gama点的振动频率,你的意思是计算Gama点的频率? 答:VASP可以任何k点的声子振动频率,不局限于Gamma点.手册上有提到的,要计算其他非Gamm点的频率需利用一个PHONON的程序. (3)用VASP计算后的OSZICAR中就直接输出自由能呀!为何还通过计算声子频率再得到自由能呢?怎样得到的呢? G. Kresse and J. Furthmüller, "Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set", Comput. Mat. Sci. 6, 15-50 (1996). 如果希望采用准谐近似来计算热学性质,可以参看这样的一篇文章(尽管是对金属单质体系的) PHYSICAL REVIEW B, VOLUME 65, 064302, 2002 Ab initio calculation of the thermal properties of Cu: Performance of the LDA and GGA