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

《VASP计算有限温度下的晶格常数.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、一般VASP算得的晶格常数是在0k下的，现在如果要算300k下的晶格常数，用VASP可否实现。 一种近似是认为二者一样。反正是在常温常压下，距离熔点很远。 一般来说，VASP计算不能加上温度，除非是做分子动力学的计算，来模拟固定在一定温度下稳定性。要计算高温下固体的能量和晶格常数、体积,可以计算出不同晶格常数（也就是体积）固体的声子振动频率，再得到自由能，利用F(V)的关系得到在某一个温度下的平衡体积。（这种办法是一种准谐近似），PRB上面有很多类似的文章。 有人会问： (1)既然涉及到温度都只能是MD了

2、,为何还能计算高温下的声子振动频率? 答：这是因为MD计算也能求解出声子的振动频率，这个在Rev.Mod.Phys.2001的一篇文章就有提到。 (2)得到声子的振动频率在VASP中是设定IBRION=5,但是只能计算Gama点的振动频率,你的意思是计算Gama点的频率? 答：VASP可以任何k点的声子振动频率，不局限于Gamma点。手册上有提到的，要计算其他非Gamm点的频率需利用一个PHONON的程序。 (3)用VASP计算后的OSZICAR中就直接输出自由能呀!为何还通过计算声子频率再得到自由能呢?

3、怎样得到的呢? 答：VASP输出的自由能只是电子部分的，非原子的自由能。看看VASP在comput.mat.sci.上面的一篇文章就知道了vasp的理论基础。G.KresseandJ.Furthmüller,"Efficiencyofab-initiototalenergycalculations formetalsandsemiconductorsusingaplane-wavebasisset",Comput.Mat.Sci.6,15-50(1996). 如果希望采用准谐近似来计算热学性质，可以参看这

4、样的一篇文章（尽管是对金属单质体系的） PHYSICALREVIEWB,VOLUME65,，2002 AbinitiocalculationofthethermalpropertiesofCu:PerformanceoftheLDAandGGA