VASP计算有限温度下的晶格常数

一般VASP算得的晶格常数是在0k下的，现在如果要算300k下的晶格常数，用VASP可否实现。一种近似是认为二者一样。反正是在常温常压下，距离熔点很远。一般来说，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中就直接输出自由能呀!为何还通过计算声子频率再得到自由能呢?怎样得到的呢?答：VASP输出的自由能只是电子部分的，非原子的自由能。看看VASP在comput.mat.sci.上面的一篇文章就知道了vasp的理论基础。G.KresseandJ.Furthmüller,Efficiencyofab-initiototalenergycalculationsformetalsandsemiconductorsusingaplane-wavebasisset,Comput.Mat.Sci.6,15-50(1996).如果希望采用准谐近似来计算热学性质，可以参看这样的一篇文章（尽管是对金属单质体系的）PHYSICALREVIEWB,VOLUME65,064302，2002AbinitiocalculationofthethermalpropertiesofCu:PerformanceoftheLDAandGGA