计算化学公社

标题: 请问周期性量化计算软件是否有类似RI加速的功能？ [打印本页]

作者
Author: ghifi37 时间: 2019-12-11 09:37
标题: 请问周期性量化计算软件是否有类似RI加速的功能？
请问周期性量化计算软件是否有类似密度拟合近似（RI）加速的功能？
VASP中没找到，好像6的manual里也没提到过。
CP2K有个名称很像的，可惜是平面波结合高斯波函数的。
要不然100来个原子supercell的relax已经几天几天的弄，实在耗不起。

作者
Author: pyscf 时间: 2019-12-11 10:31
用unit cell优化就行了...

作者
Author: ghifi37 时间: 2019-12-11 10:45

pyscf 发表于 2019-12-11 10:31
用unit cell优化就行了...

很多时候是免不了supercell的，比如计算声子、缺陷、吸附、掺杂等等。

作者
Author: 卡开发发 时间: 2019-12-11 11:10
对于KS方程中，V_H的计算，平面波直接通过Fourier变换做(4π/G^2*ρ(G))，不涉及到多中心积分。因此由基函数和算法本身决定了平面波没有RI这种做法，至少目前没看到有人在这方面特别标新立异。

平面波确实能够系统性控制地使计算达到足够高的精度，基函数尺寸也就很大，这导致计算量也很大。另外，因为LCAO的方法的局域特性，也有相应的一些加速手段，因此在大体系计算中LCAO有很大优势。

100多原子，如果没有足够高的对称性，且不说SCF比较长，构型优化需要的步数也会很多。

作者
Author: ghifi37 时间: 2019-12-11 13:01

卡开发发发表于 2019-12-11 11:10
对于KS方程中，V_H的计算，平面波直接通过Fourier变换做(4π/G^2*ρ(G))，不涉及到多中心积分。因此由基函 ...

好的，谢谢！我去查查相应的LCAO加速手段去！

作者
Author: pyscf 时间: 2019-12-12 01:23
计算声子不需要优化supercell 呵呵

作者
Author: beefly 时间: 2019-12-12 10:53
lda和纯gga的情况下，gaussian-pbc和adf-band都用了ri技术

作者
Author: chrinide 时间: 2019-12-12 12:06
TURBOMOLE执行了O(N)线性标度周期性(RI/DF)DFT代码，号称非常快-非常经济(内存消耗非常小)，目前功能并不完全(没有对称性)，而且贵要钱.
PySCF的PBC模块也执行了DF/RI技术，速度非常快，免费，功能还算齐全，你可以一试.
LSDalton的PBC模块也可以用DF/RI，速度就不知道了，功能很少。
其他的商业软件比如FHI-aims，AMS-BAND支持的比较好，功能比较完善，但是非常贵。
另外，MOLPRO也开始支持周期计算了，不过刚开始，功能很菜。

作者
Author: ghifi37 时间: 2019-12-12 15:16

beefly 发表于 2019-12-12 10:53
lda和纯gga的情况下，gaussian-pbc和adf-band都用了ri技术

谢谢回复，我比较一下！

作者
Author: ghifi37 时间: 2019-12-12 15:20

chrinide 发表于 2019-12-12 12:06
TURBOMOLE执行了O(N)线性标度周期性(RI/DF)DFT代码，号称非常快-非常经济(内存消耗非常小)，目前功能并不完 ...

感谢指导！我之前以为pyscf不好做结构优化，今天在你提醒后，看了下manual，发现在geompt中可以调用geomeTRIC和PyBerny进行结构优化，那挺好的，准备试一下。

作者
Author: wei 时间: 2019-12-12 22:38
FHI-aims 2000欧，可以有5个users

作者
Author: pyscf 时间: 2019-12-12 23:33

chrinide 发表于 2019-12-12 12:06
TURBOMOLE执行了O(N)线性标度周期性(RI/DF)DFT代码，号称非常快-非常经济(内存消耗非常小)，目前功能并不完 ...

PySCF的pbc没有梯度吧？

作者
Author: ghifi37 时间: 2020-2-7 16:40

pyscf 发表于 2019-12-12 23:33
PySCF的pbc没有梯度吧？

PySCF的PBC的确没有梯度，是不是就没法做结构优化了？

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