文献总结：

1.

虚频的位置是个高对称点吧。说明你的单胞或者对称性选的不对。
你先做点检验：
1.把对称性去掉，然后做2-3次结构优化，然后再算声子谱，看看行不行，
2.计算声子的时候，扩更大的胞，看看情况有没有好转。
如果上面两个工作都不行，那么采用如下方法：
1.在虚频方向扩胞，比如虚频的高对称点是（0.5,0,0），就扩2*1*1的胞。
2.找对应的震动方向，具体方法详见phonopy网站，一般这种现象出现在胞的交界原子的震动上面。
3.如果上面的做法还不行，用对应的超胞做分子动力学，最好用metadynamic。
4.如果上面都不行，那就只能用超胞的原子数做结构预测了，不过这个计算量太大了。

 

1.结构优化的时候在INCAR里面加ISYM=0，你在castep里面建的结构和vasp的POSCAR实际上没有任何关系。
2.如果这样，我也没办法，我的经验是好多虚频是因为长程相互作用计算不足造成的。
一个结构优化好的话，只能保证gamma点处没有虚频。非gamma点处的虚频意味着超胞没有优化好。就像我说的（0.5,0,0)有虚频意味着2*1*1的超胞是有虚频的。实际上phonopy是能够给出超胞虚频的震动方向的，不过有个问题那个过程不是一般的麻烦。

 

弄过，一是phonopy找虚频的震动方向特麻烦，二是就算你按照那个方向移动原子坐标，移动之后很可能虚频更多。
还不如分子动力学方法。其实最好是做结构预测，但是这个计算量很成问题。

 

http://emuch.net/html/201406/7591851.html

 

2.

DFPT 和线性响应是同一种方法， DFPT实际上就是基于DFT 基础上上的线性响应。 对体系给一个微扰， 在微扰很小的情况下， 电子态密度变化和微扰是线性关系， 造成态密度变化，从而得到能量变化，构造力学矩阵。常见的软件有PWscf，abaint，castep  
另外一种叫小位移法（或是frozen phonon method)，就是移动一个原子，计算其他原子对其的影响， 计算力学矩阵，从而得到声子色散。因为要考虑到粒子间的作用，所以要建超胞， 以致超胞外的原子对移动原子的影响可以忽略。常见 软件有vasp+phonopy, siesta

 

3.

个人记得热导似乎要牵扯到非谐声子，也就是ph-ph interaction~暂时没有非常标准的计算模式～
简单的通过声子结果再通过简单的公式得到热导可以参考文献：Phys Rev B, 84, 024302。
如果要比较详细的讨论和比较有物理思想的求解，最好看看文章：Phys Rev B, 80, 125203。