求助，关于TDDFT输入文件的温度和弥散参数设置

不会扣篮的后卫 · 发表于 Post on 2023-11-20 12:06:46

本帖最后由不会扣篮的后卫于 2023-11-20 21:30 编辑

大家好，我在使用CP2K TDDFT方法来计算金刚烷（Adamantane, C10H16）的孤立体系的激发能。
目前采用了两个不同的杂化泛函，B3LYP和CAM-B3LYP，算出的激发能与实验数据对比如下。

问题1. 感觉计算值和实验值相差的eV值在1到2eV之间，请问大家这个差距合理吗？

问题2. 请问各位老师，同学在CP2K TDDFT加入温度的参数设定，影响计算结果吗（打算在TDDFT计算中模拟和实验值一样的温度，25 ºC）？

问题3. 测试出金刚烷（Adamantane, C10H16）实验值的作者，也在文章中使用高斯09计算该体系。
计算方法如下，

请问大家，弥散函数该如何加到TDDFT的输入文件里呢（附上作者的文章）？

不会扣篮的后卫 · 发表于 Post on 2023-11-24 22:23:29

sobereva 发表于 2023-11-22 19:36
一看SCF收敛情况就知道当前波函数已经很离谱了，导致密度在盒子边缘都有分布，所以有那种提示。B3LYP的SI ...

谢谢卢老师的建议，我去尝试用其他杂化泛函算一算。

sobereva · 发表于 Post on 2023-11-22 19:36:05

不会扣篮的后卫发表于 2023-11-22 14:31
卢老师，您好。我按照您的建议，采用B3LYP结合aug-cc-pVTZ基组对金刚烷及其相关的衍生结构进行TDDFT计算 ...

一看SCF收敛情况就知道当前波函数已经很离谱了，导致密度在盒子边缘都有分布，所以有那种提示。B3LYP的SIE问题太大，根本不适合算里德堡激发

不会扣篮的后卫 · 发表于 Post on 2023-11-22 14:31:01

sobereva 发表于 2023-11-21 02:22
1-2eV差得太多

温度不直接影响TDDFT计算结果，仔细看此文

卢老师，您好。我按照您的建议，采用B3LYP结合aug-cc-pVTZ基组对金刚烷及其相关的衍生结构进行TDDFT计算。
TDDFT计算的孤立体系结构都使用B3LYP和6-311G**的结构优化(GEO_OPT), 输入文件统一设置为30个激发态。

TDDFT的输出结果：孤立体系的金刚烷和1-一甲基金刚烷结构的计算分别成功收敛，同时计算过程中没有报错。
但是孤立体系的1,3-二甲基金刚烷结构的计算没有成功收敛，并多次出现以下这种报错。

<div> Density non-zero on the ***</div><div> *** edges of the unit cell: wrong results in WAVELET solver ***</div>

复制代码

我尝试在原先的基础上，把输入文件基组改为6-311++G**。再次提交任务后，孤立体系的1,3-二甲基金刚烷结构的计算可以成功收敛，在TDDFT的计算中没有报错。
请问卢老师，上述报错出现的原因是什么呢？附上孤立体系，两组不同基组的1,3-二甲基金刚烷结构输入输出文件供您参考。

不会扣篮的后卫 · 发表于 Post on 2023-11-22 11:32:20

wzkchem5 发表于 2023-11-21 18:22
如果要研究温度对激发态的影响，看https://mattermodeling.stackexchange.com/questions/10922/does-it-mak ...

谢谢老师的解答！

wzkchem5 · 发表于 Post on 2023-11-21 18:22:50

如果要研究温度对激发态的影响，看https://mattermodeling.stackexch ... e-scf-a/10923#10923

不会扣篮的后卫 · 发表于 Post on 2023-11-21 13:42:30

sobereva 发表于 2023-11-21 02:22
1-2eV差得太多

温度不直接影响TDDFT计算结果，仔细看此文

谢谢卢老师！

sobereva · 发表于 Post on 2023-11-21 02:22:20

1-2eV差得太多

温度不直接影响TDDFT计算结果，仔细看此文
谈谈温度、压力、同位素设定对量子化学计算结果产生的影响
http://sobereva.com/423（http://bbs.keinsci.com/thread-10159-1-1.html）

当前会有很多里德堡态，所以必须加弥散函数。但是仅仅用6-311++G**根本不够，鉴于体系不大，至少应该用aug-cc-pVTZ

不会扣篮的后卫 · 发表于 Post on 2023-11-20 21:42:18

KSeGaSn 发表于 2023-11-20 17:26
孤立体系为啥先想着用cp2k来算了尤其是在别人用Gaussian算出来结果还不错的情况下？我个人觉得激发能和实 ...

谢谢老师的建议（ps：主要是通过这个孤立体系的TDDFT 计算，结合卢老师培训班的讲义来熟悉CP2K的激发态输入文件）。

不会扣篮的后卫 · 发表于 Post on 2023-11-20 21:31:12

zjxitcc 发表于 2023-11-20 16:39
不是“弥撒”，是弥散。基组名6-311++G(d,p)中的++就是弥散函数，不用刻意加什么东西。

谢谢老师的纠错，已修改。

KSeGaSn · 发表于 Post on 2023-11-20 17:26:15

孤立体系为啥先想着用cp2k来算了尤其是在别人用Gaussian算出来结果还不错的情况下？我个人觉得激发能和实验如果差了1~1.5eV那确实不好说明什么问题了……

zjxitcc · 发表于 Post on 2023-11-20 16:39:21

不是“弥撒”，是弥散。基组名6-311++G(d,p)中的++就是弥散函数，不用刻意加什么东西。

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[CP2K] 求助，关于TDDFT输入文件的温度和弥散参数设置