有限温度下晶格常数计算方法探讨.pdfVIP

可以使用VASP计算有限温度下的晶格常数

通常VASP计算的晶格常数是在0K下的，现在如果要计算300K下的晶格常数，用VASP是否可以实现。

一种近似方法是认为两者相同。反正是在常温常压下，距离很远。

一般来说，VASP计算不能直接加入温度，除非是进行分子动力学计算，以模拟在特定温度下的稳定

性。要计算高温下固体的能量和晶格常数、体积，可以通过计算不同晶格常数（即不同体积下）固体的声

子振动频率，再得到自由能，利用F(V)的关系得到在某一温度下的平衡体积。（这种方法是一种准谐近

似），PRB上有很多类似的文章。

有人会问:

(1)既然涉及到温度都只能是MD了，为何还能计算高温下的声子振动频率？答：这是因为MD计算

也能求解出声子的振动频率，这个在Rev.Mod.Phys.2001的一篇文章就有提到。

(2)得到声子的振动频率在VASP中是设定IBRION=5，但是只能计算Gama点的振动频率，你的意

思是计算Gama点的频率？答：VASP可以计算任何k点的声子振动频率，不局限于G点。手册

上有提到的，要计算其他非G点的频率需利用一个PHONON的程序。

(3)用VASP计算后的