计算化学公社

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

[CP2K] 求助cp2k优化激发态结构出错

[复制链接 Copy URL]

93

帖子

0

威望

557

eV
积分
650

Level 4 (黑子)

这个是报错类型Data structures of source and target full matrix are not equivalent
我是先按照了使用CP2K结合Multiwfn对周期性体系模拟UV-Vis光谱和分析电子激发态 - 思想家公社的门口:量子化学·分子模拟·二次元 (sobereva.com)找到了
&EXCITED_STATES T
     STATE 42                    
    &END EXCITED_STATES

要算的第几个激发态的结构,然后运行报错,请问是输入文件那部分有问题呢?

未命名图片.png (36.2 KB, 下载次数 Times of downloads: 13)

未命名图片.png

1.out

137.4 KB, 下载次数 Times of downloads: 4

hqcl-90.inp

12.24 KB, 下载次数 Times of downloads: 28

5万

帖子

99

威望

5万

eV
积分
112356

管理员

公社社长

2#
发表于 Post on 2022-9-21 05:28:46 | 只看该作者 Only view this author
先尝试用常规的TDDFT而非sTDA

北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

205

帖子

3

威望

1894

eV
积分
2159

Level 5 (御坂)

3#
发表于 Post on 2022-9-21 09:23:40 | 只看该作者 Only view this author
大概率是缺少以下参数和AUX基组,顺便请问一下社长使用普通泛函(PBE)可以计算激发态么,还是需要杂化泛函?杂化泛函耗时基本4-5倍提升了,着实算不太动
    &AUXILIARY_DENSITY_MATRIX_METHOD
      METHOD BASIS_PROJECTION
      ADMM_PURIFICATION_METHOD NONE
      EXCH_SCALING_MODEL NONE
      EXCH_CORRECTION_FUNC PBEX
    &END AUXILIARY_DENSITY_MATRIX_METHOD

93

帖子

0

威望

557

eV
积分
650

Level 4 (黑子)

4#
 楼主 Author| 发表于 Post on 2022-9-21 11:14:11 | 只看该作者 Only view this author
sobereva 发表于 2022-9-21 05:28
先尝试用常规的TDDFT而非sTDA

去掉sTDA还是会报一样的错误

93

帖子

0

威望

557

eV
积分
650

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2022-9-21 15:40:42 | 只看该作者 Only view this author
Eudaimonia 发表于 2022-9-21 09:23
大概率是缺少以下参数和AUX基组,顺便请问一下社长使用普通泛函(PBE)可以计算激发态么,还是需要杂化泛函? ...

我想应该不是这个问题,我在group上看到有个人的输入没有这些但是还是正常跑不会报这个错https://groups.google.com/g/cp2k/c/ZIKOl2Dq2T4/m/5rE00PfwBwAJ

5万

帖子

99

威望

5万

eV
积分
112356

管理员

公社社长

6#
发表于 Post on 2022-9-21 22:31:50 | 只看该作者 Only view this author
Eudaimonia 发表于 2022-9-21 09:23
大概率是缺少以下参数和AUX基组,顺便请问一下社长使用普通泛函(PBE)可以计算激发态么,还是需要杂化泛函? ...

普通泛函可以,例如下文提到的文章
使用CP2K结合Multiwfn对周期性体系模拟UV-Vis光谱和考察电子激发态
http://sobereva.com/634http://bbs.keinsci.com/thread-28006-1-1.html

不过用来发表的目的,但凡杂化泛函算得动都最好用杂化泛函

至于激发态优化,合适的情况用纯泛函来节约时间倒是未尝不可。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

205

帖子

3

威望

1894

eV
积分
2159

Level 5 (御坂)

7#
发表于 Post on 2022-9-22 15:17:15 | 只看该作者 Only view this author
sobereva 发表于 2022-9-21 22:31
普通泛函可以,例如下文提到的文章
使用CP2K结合Multiwfn对周期性体系模拟UV-Vis光谱和考察电子激发态
...

了解了,谢谢sob老师

10

帖子

0

威望

1829

eV
积分
1839

Level 5 (御坂)

8#
发表于 Post on 2022-12-6 10:34:10 | 只看该作者 Only view this author
为了节省空间,输入文件中有ADDED_MOS选项,会造成数据不匹配。可以把那行注释掉试一试

268

帖子

0

威望

2288

eV
积分
2556

Level 5 (御坂)

9#
发表于 Post on 2024-3-1 13:16:48 | 只看该作者 Only view this author
sobereva 发表于 2022-9-21 22:31
普通泛函可以,例如下文提到的文章
使用CP2K结合Multiwfn对周期性体系模拟UV-Vis光谱和考察电子激发态
...

卢老师,您好。我在用CP2K PBE0泛函算周期性体系的激发态优化任务。请问这个报错是杂化泛函算不动激发态优化的情况吗?

5万

帖子

99

威望

5万

eV
积分
112356

管理员

公社社长

10#
发表于 Post on 2024-3-2 06:32:31 | 只看该作者 Only view this author
不会扣篮的后卫 发表于 2024-3-1 13:16
卢老师,您好。我在用CP2K PBE0泛函算周期性体系的激发态优化任务。请问这个报错是杂化泛函算不动激发态 ...

估计是ADMM设置不当,没具体文件没法说
确保用的是Multiwfn最新版本产生输入文件,以及CP2K用的是较新版本,此时肯定是兼容的
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

268

帖子

0

威望

2288

eV
积分
2556

Level 5 (御坂)

11#
发表于 Post on 2024-3-2 23:40:06 | 只看该作者 Only view this author
sobereva 发表于 2024-3-2 06:32
估计是ADMM设置不当,没具体文件没法说
确保用的是Multiwfn最新版本产生输入文件,以及CP2K用的是较新版 ...

卢老师,我附上输入和输出文件,麻烦您帮忙看看。

BN.inp

5.44 KB, 下载次数 Times of downloads: 6

out.log

26.26 KB, 下载次数 Times of downloads: 1

5万

帖子

99

威望

5万

eV
积分
112356

管理员

公社社长

12#
发表于 Post on 2024-3-2 23:54:55 | 只看该作者 Only view this author
不会扣篮的后卫 发表于 2024-3-2 23:40
卢老师,我附上输入和输出文件,麻烦您帮忙看看。

明确提示了
*** WARNING in qs_tddfpt2_fhxc_forces.F:577 :: Forces need symmetric ADMM ***
*** kernel corrections

把输入文件里的#      ADMM_KERNEL_CORRECTION_SYMMETRIC T的注释去掉。本来Multiwfn产生跟受力计算有关的任务用了ADMM时输入文件里这就是开启的,不要自己随便改。

在北京科音CP2K第一性原理计算培训班(http://www.keinsci.com/workshop/KFP_content.html)里我专门强调了

北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

268

帖子

0

威望

2288

eV
积分
2556

Level 5 (御坂)

13#
发表于 Post on 2024-3-3 02:58:23 | 只看该作者 Only view this author
sobereva 发表于 2024-3-2 23:54
明确提示了
*** WARNING in qs_tddfpt2_fhxc_forces.F:577 :: Forces need symmetric ADMM ***
*** k ...

卢老师,我按照您的建议,取消受力的注释。但还是有了新的报错,麻烦您再帮我看看报错原因。附上相关的输入和输出文件。

BN.inp

5.44 KB, 下载次数 Times of downloads: 6

out.log

34.16 KB, 下载次数 Times of downloads: 3

5万

帖子

99

威望

5万

eV
积分
112356

管理员

公社社长

14#
发表于 Post on 2024-3-3 03:58:38 | 只看该作者 Only view this author
不会扣篮的后卫 发表于 2024-3-3 02:58
卢老师,我按照您的建议,取消受力的注释。但还是有了新的报错,麻烦您再帮我看看报错原因。附上相关的输 ...

The linear solver didn't converge! Maximum number of iterations reached.
没收敛

这么小晶胞不考虑k点没意义,而CP2K的TDDFT不支持考虑k点的情况
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

268

帖子

0

威望

2288

eV
积分
2556

Level 5 (御坂)

15#
发表于 Post on 2024-3-3 09:51:16 | 只看该作者 Only view this author
sobereva 发表于 2024-3-3 03:58
The linear solver didn't converge! Maximum number of iterations reached.
没收敛

谢谢卢老师。我按照您的建议,扩大了晶胞。再次提交任务,出现了这个内存不够的报错。

是按照这个submit.slm提交的任务。附上报错slurm.out文件,请问卢老师有什么修改建议避免这个报错吗?

slurm.out

5.14 KB, 下载次数 Times of downloads: 0

本版积分规则 Credits rule

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

GMT+8, 2024-11-24 18:38 , Processed in 0.194527 second(s), 24 queries , Gzip On.

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