Originally posted by dxcharlary at 2011-05-29 20:21:14:
说的很好，补充几点：
5用的是差分，6用的是DFPT，感觉好像DFPT更靠谱一点。不过貌似VASP的计算结果没有PWSCF的好，而且只能算gamma点的。
dx,dy,dz是本征矢，也是震动方向，有时做一些强度或者是稳定性计算 ...

Originally posted by dxcharlary at 2011-05-29 20:23:02:
对了，在这些东西的上面会有你要的动力学矩阵。

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    ...

Originally posted by donggaomin at 2011-05-29 21:53:15:
您好！

请问：1f,2f,3f...等是不是类似于能带图，1f表示其中的一条声频色散带在gamma点上的值？2f表示另一条声频色散带在gamma点上的值？...


我们在文章中常看到的“Gamma点A1(1TO)模式声子最低频率的 ...

Originally posted by donggaomin at 2011-05-29 22:57:23:
请问stractor,
     E（1TO），A1（1TO），F(1TO)  我现在的理解是TO表示横波光学模声子，而1TO中的1表示能带一的，不知道我的理解对不对？

另外，E，A1，F 分别表示什么？从数据中怎么判断？


谢谢！

Originally posted by donggaomin at 2011-05-29 22:01:55:
我看了一下，这些东西上面是“SECOND DERIVATIVES (NOT SYMMETRIZED)”，（在IBRION = 6 时得到的）。


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

475878楼: Originally posted by stractor at 2011-05-29 21:18:00
是的，有限差分的可信度取决于原子的位移
因为有限差分计算Hessen矩阵是
Hij=d2E/dxidxj=/dxidxj
上面的d是微分，rjo就是j原子的平衡位置，E是能量。

上式反映出位移xi, xj不能太大。
太小的话位移后构型的能 ...