24小时热门版块排行榜    

查看: 2400  |  回复: 17
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

donggaomin

金虫 (正式写手)

[求助] 请问VASP5.2中如何设置INCAR才可以计算(gamma点的)phonon frequency?

请问VASP5.2中如何设置INCAR才可以计算(gamma点的)phonon frequency?

说明书上只说了设置IBRION = 6

另外,怎么查看结果?论坛里提到是动力学矩阵的对角化,不是很明白。我知道动力学矩阵对角化后是频率的平方。

谢谢!

VASP5.2中是不是可以计算phonon frequency?
回复此楼
好好学习
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dxcharlary

专家顾问 (职业作家)

【答案】应助回帖

★ ★
donggaomin(金币+5): 谢谢dxcharlary!能不能举例说明一下关于强度和稳定性的计算? 2011-05-29 21:26:56
zzy870720z(金币+2): 谢谢指教 2011-05-29 23:32:27
引用回帖:
Originally posted by stractor at 2011-05-29 20:00:30:
举例
1 f  =   15.937321 THz   100.137142 2PiTHz  531.611795 cm-1    65.911486 meV
             X         Y         Z           dx          dy          dz
     -0.000113  1.858561  0.000000    ...

说的很好,补充几点:
5用的是差分,6用的是DFPT,感觉好像DFPT更靠谱一点。不过貌似VASP的计算结果没有PWSCF的好,而且只能算gamma点的。
dx,dy,dz是本征矢,也是震动方向,有时做一些强度或者是稳定性计算的时候可能会用到。
7楼2011-05-29 20:21:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dxcharlary

专家顾问 (职业作家)

【答案】应助回帖


donggaomin(金币+5): 2011-05-29 21:34:28
zzy870720z(金币+1): 谢谢补充 2011-05-29 23:32:41
对了,在这些东西的上面会有你要的动力学矩阵。
8楼2011-05-29 20:23:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dxcharlary

专家顾问 (职业作家)

【答案】应助回帖

donggaomin(金币+2): 谢谢dxcharlary!!! 2011-05-30 16:39:59
引用回帖:
Originally posted by donggaomin at 2011-05-29 22:01:55:
我看了一下,这些东西上面是“SECOND DERIVATIVES (NOT SYMMETRIZED)”,(在IBRION = 6 时得到的)。


我现在的理解是“力常数矩阵就是HESSIAN矩阵,矩阵元是能量对位移的二阶倒数。动力学矩阵是力常数矩 ...

那个东西是hessian矩阵,k=0的时候,exp(-kr)=1.且单胞内的原子相互作用要远大于跨元胞之间的相互作用,所以动力学矩阵退化为hessian矩阵除以质量的平方根。
另外不好意思,告诉你错了。。。5,6是差分,7,8是DFPT,5,7是不考虑对称性,6,8是考虑对称性的。
16楼2011-05-30 14:28:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 donggaomin 的主题更新
信息提示
请填处理意见