计算化学公社

标题: hse06+admm dzvp配合dftd3矫正计算band报错求助 [打印本页]

作者
Author: Aridea 时间: 2021-8-19 10:46
标题: hse06+admm dzvp配合dftd3矫正计算band报错求助
各位大佬，最近本人尝试用杂化泛函配合admm加速来计算范德华异质结的band，但是一旦读完k path之后就开始报错，报错信息我贴在了下面了，为了解决这个问题我尝试过了切换了OT算法到diagonalization算法，包括他们restart的波函数也同时切换了的，一样报错；我还尝试过不用restart直接计算也不管用；也试过提升CUTOFF从400到450到500（胞长22 A），EPS_DEFAULT E-10到E-12; EPS_SCHWARZ E-8到E-10都不太行；
之前几乎同样的配置用pbe泛函计算没问题的只是band gap过于低估，不知各位大佬能够给点建议或意见，或者有大佬之前算过hse06+admm下的band的经验也烦请指点一二，万分感谢！

作者
Author: 卡开发发 时间: 2021-8-19 14:04
这里应该只是个Wanring，这个cutoff应该指的是杂化泛函计算中，Coulomb算符有个截断半径，而不是网格的截断，程序建议这个截断半径不要超过晶格最短边的一半。

作者
Author: Aridea 时间: 2021-8-19 14:07

卡开发发发表于 2021-8-19 14:04
这里应该只是个Wanring，这个cutoff应该指的是杂化泛函计算中，Coulomb算符有个截断半径，而不是网格的截断 ...

大大，官网pbe0有cutoff r的设置不能忽略；而hse06没这个设置所以这个warning可以忽略的，现在最关心scf收敛后开始band计算咋一直报错，倒腾了好久没搞定

作者
Author: 卡开发发 时间: 2021-8-19 18:52

Aridea 发表于 2021-8-19 14:07
大大，官网pbe0有cutoff r的设置不能忽略；而hse06没这个设置所以这个warning可以忽略的，现在最关心scf ...

按道理HSE06也可以做truncated Coulomb，不过我不知道他们程序咋实现的，当时只是看了他们文章。他们程序有些地方乱的很。报错我没看到他详细的提示，暂时判断不出来。

作者
Author: Aridea 时间: 2021-8-19 19:24
本帖最后由 Aridea 于 2021-8-19 19:58 编辑

嗯嗯，这个cp2k确实挺乱的,好多奇怪的地方，例如我用pbe算band和pdos，fermi能直接差了1-2eV,不知道啥情况，同一个input同一个scf算出来的费米能居然对不上

作者
Author: leeru 时间: 2021-8-19 22:30
我也是换了好多组合，都不行，收敛都不收敛。

作者
Author: leeru 时间: 2021-8-20 00:06
我感觉cp2k是不支持HSE06的K点band structure计算，我去google group搜了一下，也有人问过类似的问题。

作者
Author: Aridea 时间: 2021-8-20 09:35

leeru 发表于 2021-8-20 00:06
我感觉cp2k是不支持HSE06的K点band structure计算，我去google group搜了一下，也有人问过类似的问题。

兄弟，一定要搞个收敛的试试哈，目前我是收敛了读取k path就开始报错退出，如果兄弟的scf收敛了也和我一样的报错那就是google group说的那样了。到时候记得和我说哈

作者
Author: leeru 时间: 2021-8-23 12:44
采用

作者
Author: leeru 时间: 2021-8-23 12:46
我的情况是采用gamma点hse06 ADMM可以收敛，但一旦加入K sampling就报错

作者
Author: Aridea 时间: 2021-8-23 19:25
嘿嘿，多谢兄弟我也是这样的

作者
Author: 卡开发发 时间: 2021-9-4 11:58

Aridea 发表于 2021-8-23 19:25
嘿嘿，多谢兄弟我也是这样的

要是cp2k杂化泛函只能算gamma的那我无能为力了。

一般程序如果没经过设计，对于杂化泛函也不能读取密度指定KPATH，因为杂化泛函没办法使用Harris-Foulkes泛函来通过非自洽得到高对称点上的本征态，所以需要手动指定K网格和权重，以及零权重的高对称点（实际这方面的操作很多程序在体验上还有优化的余地，不过他们不去做罢了）。

作者
Author: Aridea 时间: 2021-9-5 23:55

卡开发发发表于 2021-9-4 11:58
要是cp2k杂化泛函只能算gamma的那我无能为力了。

一般程序如果没经过设计，对于杂化泛函也不能读取密 ...

原来如此，涨姿势了

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