castep计算声子谱内存不足

菩城

各位老师好，我的体系有60个原子，我想计算声子谱，用castep跑不动，请问vasp是否可以呢，这么大体系是否只能用超算了呢@卡开发发

卡开发发

这个问题说的很笼统，你得说清你的体系、你用的具体方法、服务器的配置等。

CASTEP做声子计算有两种不同的方式：线性响应 or 有限位移。一般有限位移的方法就是数值导数的方法，这个计算需要设置一个cutoff的方法，cutoff的大小决定你系统采用超胞的大小；而线性响应则是通过DFPT，这个过程不必建立超胞，而是转到倒空间求各个q点的声子谱，因此两者内存用量还是有差别。但CASTEP的线性响应局限性很多，如不能使用DFT-D2以及必须使用模守恒赝势等。

VASP这个程序自身线性响应方法只能做q=Gamma的声子谱，得借助phonopy来做，也得建立超胞，计算量也不会小很多。

除了上述主流程序，也可以考虑QE做线性响应，目前新的版本貌似已经支持vdW-DF这类泛函通过线性响应计算声子谱。超胞法的话可以考虑SIESTA，自洽计算的话单节点大概两三百原子应该没啥压力，缺点是结构优化可能有点捉鸡。

菩城

卡开发发 发表于 2016-10-8 23:26
这个问题说的很笼统，你得说清你的体系、你用的具体方法、服务器的配置等。

CASTEP做声子计算有两种不同 ...

好详细，非常感谢老师！按照您说的我试试。