计算化学公社

标题: castep计算声子谱内存不足 [打印本页]

作者
Author: 菩城 时间: 2016-10-8 16:04
标题: castep计算声子谱内存不足
各位老师好，我的体系有60个原子，我想计算声子谱，用castep跑不动，请问vasp是否可以呢，这么大体系是否只能用超算了呢@卡开发发

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

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

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

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

作者
Author: 菩城 时间: 2016-10-9 19:44

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

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

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

欢迎光临计算化学公社 (http://ccc.keinsci.com/)