计算化学公社

 找回密码 Forget password
 注册 Register
Views: 6882|回复 Reply: 12
打印 Print 上一主题 Last thread 下一主题 Next thread

[CP2K] hse06+admm dzvp配合dftd3矫正计算band报错求助

[复制链接 Copy URL]

148

帖子

0

威望

888

eV
积分
1036

Level 4 (黑子)

各位大佬,最近本人尝试用杂化泛函配合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的经验也烦请指点一二,万分感谢!



hse06-admm-band报错.png (86.87 KB, 下载次数 Times of downloads: 28)

hse06-admm-band报错.png

收敛信息.png (73.28 KB, 下载次数 Times of downloads: 24)

收敛信息.png

3809

帖子

3

威望

1万

eV
积分
20334

Level 6 (一方通行)

围观吃瓜群众

2#
发表于 Post on 2021-8-19 14:04:01 | 只看该作者 Only view this author
这里应该只是个Wanring,这个cutoff应该指的是杂化泛函计算中,Coulomb算符有个截断半径,而不是网格的截断,程序建议这个截断半径不要超过晶格最短边的一半。
日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。不做培*,不接代*,不接*发谢谢。

148

帖子

0

威望

888

eV
积分
1036

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2021-8-19 14:07:07 | 只看该作者 Only view this author
卡开发发 发表于 2021-8-19 14:04
这里应该只是个Wanring,这个cutoff应该指的是杂化泛函计算中,Coulomb算符有个截断半径,而不是网格的截断 ...

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

3809

帖子

3

威望

1万

eV
积分
20334

Level 6 (一方通行)

围观吃瓜群众

4#
发表于 Post on 2021-8-19 18:52:52 | 只看该作者 Only view this author
Aridea 发表于 2021-8-19 14:07
大大,官网pbe0有cutoff r的设置不能忽略;而hse06没这个设置所以这个warning可以忽略的,现在最关心scf ...

按道理HSE06也可以做truncated Coulomb,不过我不知道他们程序咋实现的,当时只是看了他们文章。他们程序有些地方乱的很。报错我没看到他详细的提示,暂时判断不出来。
日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。不做培*,不接代*,不接*发谢谢。

148

帖子

0

威望

888

eV
积分
1036

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2021-8-19 19:24:23 | 只看该作者 Only view this author
本帖最后由 Aridea 于 2021-8-19 19:58 编辑

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

84

帖子

0

威望

626

eV
积分
710

Level 4 (黑子)

6#
发表于 Post on 2021-8-19 22:30:16 来自手机 | 只看该作者 Only view this author
我也是换了好多组合,都不行,收敛都不收敛。

84

帖子

0

威望

626

eV
积分
710

Level 4 (黑子)

7#
发表于 Post on 2021-8-20 00:06:11 来自手机 | 只看该作者 Only view this author
我感觉cp2k是不支持HSE06的K点band structure计算,我去google group搜了一下,也有人问过类似的问题。

148

帖子

0

威望

888

eV
积分
1036

Level 4 (黑子)

8#
 楼主 Author| 发表于 Post on 2021-8-20 09:35:27 | 只看该作者 Only view this author
leeru 发表于 2021-8-20 00:06
我感觉cp2k是不支持HSE06的K点band structure计算,我去google group搜了一下,也有人问过类似的问题。

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

84

帖子

0

威望

626

eV
积分
710

Level 4 (黑子)

9#
发表于 Post on 2021-8-23 12:44:21 来自手机 | 只看该作者 Only view this author
采用

84

帖子

0

威望

626

eV
积分
710

Level 4 (黑子)

10#
发表于 Post on 2021-8-23 12:46:11 来自手机 | 只看该作者 Only view this author
我的情况是采用gamma点hse06 ADMM可以收敛,但一旦加入K sampling就报错

148

帖子

0

威望

888

eV
积分
1036

Level 4 (黑子)

11#
 楼主 Author| 发表于 Post on 2021-8-23 19:25:07 | 只看该作者 Only view this author
嘿嘿,多谢兄弟我也是这样的

3809

帖子

3

威望

1万

eV
积分
20334

Level 6 (一方通行)

围观吃瓜群众

12#
发表于 Post on 2021-9-4 11:58:53 | 只看该作者 Only view this author
Aridea 发表于 2021-8-23 19:25
嘿嘿,多谢兄弟我也是这样的

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

一般程序如果没经过设计,对于杂化泛函也不能读取密度指定KPATH,因为杂化泛函没办法使用Harris-Foulkes泛函来通过非自洽得到高对称点上的本征态,所以需要手动指定K网格和权重,以及零权重的高对称点(实际这方面的操作很多程序在体验上还有优化的余地,不过他们不去做罢了)。
日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。不做培*,不接代*,不接*发谢谢。

148

帖子

0

威望

888

eV
积分
1036

Level 4 (黑子)

13#
 楼主 Author| 发表于 Post on 2021-9-5 23:55:05 | 只看该作者 Only view this author
卡开发发 发表于 2021-9-4 11:58
要是cp2k杂化泛函只能算gamma的那我无能为力了。

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

原来如此,涨姿势了

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-18 18:32 , Processed in 0.181559 second(s), 23 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list