计算化学公社

标题: 求助cp2k优化激发态结构出错 [打印本页]

作者
Author: wmy 时间: 2022-9-20 21:56
标题: 求助cp2k优化激发态结构出错
这个是报错类型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

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

作者
Author: sobereva 时间: 2022-9-21 05:28
先尝试用常规的TDDFT而非sTDA

作者
Author: Eudaimonia 时间: 2022-9-21 09:23
大概率是缺少以下参数和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

作者
Author: wmy 时间: 2022-9-21 11:14

sobereva 发表于 2022-9-21 05:28
先尝试用常规的TDDFT而非sTDA

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

作者
Author: wmy 时间: 2022-9-21 15:40

Eudaimonia 发表于 2022-9-21 09:23
大概率是缺少以下参数和AUX基组，顺便请问一下社长使用普通泛函(PBE)可以计算激发态么，还是需要杂化泛函？ ...

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

作者
Author: sobereva 时间: 2022-9-21 22:31

Eudaimonia 发表于 2022-9-21 09:23
大概率是缺少以下参数和AUX基组，顺便请问一下社长使用普通泛函(PBE)可以计算激发态么，还是需要杂化泛函？ ...

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

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

至于激发态优化，合适的情况用纯泛函来节约时间倒是未尝不可。

作者
Author: Eudaimonia 时间: 2022-9-22 15:17

sobereva 发表于 2022-9-21 22:31
普通泛函可以，例如下文提到的文章
使用CP2K结合Multiwfn对周期性体系模拟UV-Vis光谱和考察电子激发态
...

了解了，谢谢sob老师

作者
Author: Albert0906 时间: 2022-12-6 10:34
为了节省空间，输入文件中有ADDED_MOS选项，会造成数据不匹配。可以把那行注释掉试一试

作者
Author: 不会扣篮的后卫 时间: 2024-3-1 13:16

sobereva 发表于 2022-9-21 22:31
普通泛函可以，例如下文提到的文章
使用CP2K结合Multiwfn对周期性体系模拟UV-Vis光谱和考察电子激发态
...

卢老师，您好。我在用CP2K PBE0泛函算周期性体系的激发态优化任务。请问这个报错是杂化泛函算不动激发态优化的情况吗？
(, 下载次数 Times of downloads: 10)

作者
Author: sobereva 时间: 2024-3-2 06:32

不会扣篮的后卫发表于 2024-3-1 13:16
卢老师，您好。我在用CP2K PBE0泛函算周期性体系的激发态优化任务。请问这个报错是杂化泛函算不动激发态 ...

估计是ADMM设置不当，没具体文件没法说
确保用的是Multiwfn最新版本产生输入文件，以及CP2K用的是较新版本，此时肯定是兼容的

作者
Author: 不会扣篮的后卫 时间: 2024-3-2 23:40

sobereva 发表于 2024-3-2 06:32
估计是ADMM设置不当，没具体文件没法说
确保用的是Multiwfn最新版本产生输入文件，以及CP2K用的是较新版 ...

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

作者
Author: sobereva 时间: 2024-3-2 23:54

不会扣篮的后卫发表于 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）里我专门强调了
(, 下载次数 Times of downloads: 13)

作者
Author: 不会扣篮的后卫 时间: 2024-3-3 02:58

sobereva 发表于 2024-3-2 23:54
明确提示了
*** WARNING in qs_tddfpt2_fhxc_forces.F:577 :: Forces need symmetric ADMM ***
*** k ...

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

作者
Author: sobereva 时间: 2024-3-3 03:58

不会扣篮的后卫发表于 2024-3-3 02:58
卢老师，我按照您的建议，取消受力的注释。但还是有了新的报错，麻烦您再帮我看看报错原因。附上相关的输 ...

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

这么小晶胞不考虑k点没意义，而CP2K的TDDFT不支持考虑k点的情况

作者
Author: 不会扣篮的后卫 时间: 2024-3-3 09:51

sobereva 发表于 2024-3-3 03:58
The linear solver didn't converge! Maximum number of iterations reached.
没收敛

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

是按照这个submit.slm提交的任务。附上报错slurm.out文件，请问卢老师有什么修改建议避免这个报错吗？
(, 下载次数 Times of downloads: 12)

作者
Author: sobereva 时间: 2024-3-3 10:50

不会扣篮的后卫发表于 2024-3-3 09:51
谢谢卢老师。我按照您的建议，扩大了晶胞。再次提交任务，出现了这个内存不够的报错。

是按照这个subm ...

slurm的事我不清楚

作者
Author: 不会扣篮的后卫 时间: 2024-3-3 10:53

sobereva 发表于 2024-3-3 10:50
slurm的事我不清楚

没关系，卢老师。我再去研究一下。

作者
Author: 不会扣篮的后卫 时间: 2024-4-26 11:38
老师，您好。请问您尝试用常规的TDDFT后，解决这个报错了吗？

作者
Author: logzzz 时间: 2024-4-28 11:32
内存不够就是硬伤。解决办法就是换个机器，或者花费大量的时间。

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