如何用Vasp计算声子谱？

发信人: valenhou (Octave), 信区: Gaussian
标 题: Re: [求助]如何用Vasp计算声子谱？
发信站: BBS 大话西游站 (Wed Nov 30 09:50:44 2005), 转信

还是用pwscf或abinit好一些。
目前VASP能计算声子谱的都是采用一种直接的方法：

构造超原胞，把原子移动一下，计算原胞中所有原子所受的力（这个根据体系的周期性，要多移动几个原子），然后根据这个力构造力常数矩阵。
而且一般情况下对LO-TO的split不能计算出（只有在计算了Born effective charge之后，
进一步考虑了 non-analyticity term，才能计算出）

（PHONON就是做这事的，先构造超原胞，根据体系的对称性，告诉要移动原胞中哪些原子，
下面就是每移动一个原子，然后计算原子所受的力，最后把所有的力的收集起来。把这些力
输入到PHONON中，然后就可以计算声子谱了。如前所提到的，可以再计算Born effectiv
charge，把这个也输入到PHONON，这样再计算时，就可以考虑 non-analyticity term，得到 LO-TO的splitting。）

象pwscf和abinit就不同了，采用的是一种linear response theory 的方法（或者称为
density perturbation functional theory，DFPT），直接计算出原子的移动而导致
的势场变化，再进一步构造出动力学矩阵。这种方法在计算谱时，Born effective charge
（对极性的材料）和声子谱都能计算出。

【 在 Yggdrasil (世界之树) 的大作中提到: 】
不要惦记用castep算最好，奇慢无比，而且设置不好极易出虚频
我算一个4原子的单胞居然用了一个多月……虽然机器是挫了点……

【 在 O50249 (icp) 的大作中提到: 】
: MS好像就能计算声子谱

--

--
_o \o_ __| \ / |__ o _ o/ \o/
__|- __/ \__/o \o | o/ o/__ /\ /| |
> > / \ ( \ /o\ / ) | (\ / | < \ / \

※ 来源:・BBS 大话西游站 wforum.zixia.net・[FROM: 202.120.224.18]