求助，使用CP2K计算C60的激发态时，振子强度值过小

不会扣篮的后卫

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


同时，我想比较CP2K和NWChem的激发态计算结果。于是用NWChem程序，在相同的计算级别下跑相同的结构。
发现在NWChem的输出文件中，振子强度值也是过小。

请问各位老师，C60过小的振子强度是不是不合理，请问该怎么修改输入文件来获得合理的振子强度呢?


还有一个问题是关于激发能的值。如下图所示，

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

sobereva

你可以说很小，但如果说过小，必须给定一个参考标准，怎么判断的过小，有什么证据认为这是不合理的

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

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

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


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

snljty2

别的我不懂，不过孤立的C60分子属于Ih点群，对称性极高，有一堆不可约表示，按照群论，应该会有好多好多跃迁是对称性禁阻的，振子强度严格为0。

不会扣篮的后卫

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

这不 ...

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

&SUBSYS

    &CELL

         PERIODIC NONE  

&POISSON

  PERIODIC NONE

  PSOLVER WAVELET

提交任务想要跑出纯泛函PBE的波函数当初猜，但是出现了这个报错。


结合CP2K官网的手册，



我尝试把PSOLVER WAVELET换成PSOLVER IMPLICIT。会出现新的报错，


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

不会扣篮的后卫

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

谢谢老师的解答!

sobereva

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

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

不会扣篮的后卫

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

谢谢卢老师！