计算化学公社

标题: 求助，使用CP2K计算C60的激发态时，振子强度值过小 [打印本页]

作者
Author: 不会扣篮的后卫 时间: 2023-6-27 13:21
标题: 求助，使用CP2K计算C60的激发态时，振子强度值过小
本帖最后由不会扣篮的后卫于 2023-6-27 13:23 编辑

各位老师好，我用优化好的C60 unit cell结构跑TDDFT的计算，来算C60的20个态的激发能。
在输入文件中使用杂化泛函PBE0,基组使用6-31G*和gamma point。结合CP2K的输出文件的数据，发现振子强度值过小。
(, 下载次数 Times of downloads: 21)

同时，我想比较CP2K和NWChem的激发态计算结果。于是用NWChem程序，在相同的计算级别下跑相同的结构。
发现在NWChem的输出文件中，振子强度值也是过小。
(, 下载次数 Times of downloads: 21)
请问各位老师，C60过小的振子强度是不是不合理，请问该怎么修改输入文件来获得合理的振子强度呢?

还有一个问题是关于激发能的值。如下图所示，
(, 下载次数 Times of downloads: 20)
虽然是相同的分子结构，计算级别，激发能分布（态数1到10为一组，11到15为第二组，16到20为第三组），但是两个程序CP2K和NWChem还是计算出不同激发能。
请问各位老师，两个程序间的激发能值差距大吗？（附上相关的输入和输出文件）

作者
Author: sobereva 时间: 2023-6-28 00:16
你可以说很小，但如果说过小，必须给定一个参考标准，怎么判断的过小，有什么证据认为这是不合理的

这不叫C60 unit cell，这会令别人以为你算的是C60分子晶体，而你只不过算了个孤立的C60而已

记得CP2K的TDDFT是基于TDA-DFT近似的

CP2K有没有读纯泛函波函数当初猜要说清楚，没读则杂化泛函用SCREEN_ON_INITIAL_P T时结果可能不可靠

另外，明明是孤立体系，干嘛不让Multiwfn产生CP2K输入文件时当非周期性计算，还省得用那么大盒子浪费时间。孤立体系杂化泛函计算也根本没必要用库仑截断。

作者
Author: snljty2 时间: 2023-6-28 20:33
别的我不懂，不过孤立的C60分子属于Ih点群，对称性极高，有一堆不可约表示，按照群论，应该会有好多好多跃迁是对称性禁阻的，振子强度严格为0。

作者
Author: 不会扣篮的后卫 时间: 2023-6-29 05:41
本帖最后由不会扣篮的后卫于 2023-6-29 05:44 编辑

sobereva 发表于 2023-6-28 00:16
你可以说很小，但如果说过小，必须给定一个参考标准，怎么判断的过小，有什么证据认为这是不合理的

这不 ...

谢谢卢老师的建议。
还有一个问题需要请教卢老师。
我把CP2K孤立体系的C60的输入文件的周期性参数取消的话,就如下所示

&SUBSYS

&CELL

PERIODIC NONE

&POISSON

PERIODIC NONE

PSOLVER WAVELET

提交任务想要跑出纯泛函PBE的波函数当初猜，但是出现了这个报错。
(, 下载次数 Times of downloads: 24)

结合CP2K官网的手册，
(, 下载次数 Times of downloads: 24)

我尝试把PSOLVER WAVELET换成PSOLVER IMPLICIT。会出现新的报错，
(, 下载次数 Times of downloads: 24)

请问卢老师，这两个报错在计算孤立体系的C60时可以避免吗？（附上相关输入，输出文件）

作者
Author: 不会扣篮的后卫 时间: 2023-6-29 05:45

snljty2 发表于 2023-6-28 20:33
别的我不懂，不过孤立的C60分子属于Ih点群，对称性极高，有一堆不可约表示，按照群论，应该会有好多好多跃 ...

谢谢老师的解答!

作者
Author: sobereva 时间: 2023-6-29 06:28

不会扣篮的后卫发表于 2023-6-29 05:41
谢谢卢老师的建议。
还有一个问题需要请教卢老师。
我把CP2K孤立体系的C60的输入文件的周期性参数取消 ...

wavelet必须用立方盒子，CP2K使用常识。对于当前版本Multiwfn产生CP2K输入文件，你设成非周期性默认就会用wavelet结合立方盒子

作者
Author: 不会扣篮的后卫 时间: 2023-6-30 23:54

sobereva 发表于 2023-6-29 06:28
wavelet必须用立方盒子，CP2K使用常识。对于当前版本Multiwfn产生CP2K输入文件，你设成非周期性默认就会 ...

谢谢卢老师！

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