跑绝热的MD，追踪state 1的计算时，出现了Cholesky decompose failed

不会扣篮的后卫 · 发表于 Post on 2024-3-4 12:54:45

大家好。如题，我在用孤立体系，跑绝热的NVT MD，追踪激发态state 1的计算时，出现了这个报错

请问大家有什么修改建议吗？目前的打算是跑出这个计算的轨迹文件，用于后续的计算（附上输入，输出文件供大家参考）。

sobereva · 发表于 Post on 2024-3-5 02:27:22

B3LYP对里德堡激发描述很烂，换成wB97X或者全局高HF成份的BHandHLYP
尝试把盒子设得更大。里德堡激发会使得电子高度弥散

不会扣篮的后卫 · 发表于 Post on 2024-3-5 12:01:12

sobereva 发表于 2024-3-5 02:27
B3LYP对里德堡激发描述很烂，换成wB97X或者全局高HF成份的BHandHLYP
尝试把盒子设得更大。里德堡激发会使 ...

谢谢卢老师！

不会扣篮的后卫 · 发表于 Post on 2024-4-15 12:54:11

sobereva 发表于 2024-3-5 02:27
B3LYP对里德堡激发描述很烂，换成wB97X或者全局高HF成份的BHandHLYP
尝试把盒子设得更大。里德堡激发会使 ...

卢老师，您好。我按照您的建议，增大了盒子，外加使用了wB97x泛函，但是出现了新的报错。

您有空的话，请问能帮我看看这个报错原因吗？附上输入，输出文件。

sobereva · 发表于 Post on 2024-4-16 01:07:21

不会扣篮的后卫发表于 2024-4-15 12:54
卢老师，您好。我按照您的建议，增大了盒子，外加使用了wB97x泛函，但是出现了新的报错。

分清楚warning和error
尝试减小EPS_SCF、TDDFT的CONVERGENCE、EPS_DEFAULT、增大CUTOFF
如果还这样，尝试诸如CAM-B3LYP，如果还这样，暂且无视，继续观望

不会扣篮的后卫 · 发表于 Post on 2024-4-17 12:43:20

本帖最后由不会扣篮的后卫于 2024-4-17 12:48 编辑

sobereva 发表于 2024-4-16 01:07
分清楚warning和error
尝试减小EPS_SCF、TDDFT的CONVERGENCE、EPS_DEFAULT、增大CUTOFF
如果还这样，尝 ...

卢老师，您好。我按照您的建议使用CAM-B3LYP杂化泛函，然后尝试修改您之前提到的参数来跑新的激发态MD，但还是出现了新的报错信息（附上相关的文件）。

来自文件夹CAM_B3LYP_S1_excited_MD的info.out文件的报错

还有文件out.log的warning信息以及振子强度为0的情况，

因为这个C10H16的孤立体系，优化激发态结构（文件夹OPT_S1）或者跑TDDFT的MD
（文件夹CAM_B3LYP_TDDFT_MD, 不加以下这个&EXCITED_STATES的部分）计算，都可以成功跑出输出文件。

&EXCITED_STATES
STATE -1
&END EXCITED_STATES

但是跑激发态的MD计算就一直报错。我有两个问题想请教卢老师，
1.请问卢老师可以帮忙看看文件夹CAM_B3LYP_S1_excited_MD出错的原因吗？
2.学生想请问一个愚蠢的问题。文件夹CAM_B3LYP_TDDFT_MD的计算，虽然加了&TDDFT的部分，但是没有加&EXCITED_STATES 的部分，
所以就不会指定感兴趣的激发态或者也不会有跟踪激发态的计算。请问卢老师，这样的计算产生的xyz轨迹文件和基态的MD计算产生的xyz轨迹有区别吗？

sobereva · 发表于 Post on 2024-4-18 07:13:05

用CP2K最新版本再试

实在搞不定就换成ORCA跑TDDFT的AIMD

不会扣篮的后卫 · 发表于 Post on 2024-4-19 22:02:15

sobereva 发表于 2024-4-18 07:13
用CP2K最新版本再试

实在搞不定就换成ORCA跑TDDFT的AIMD

感谢卢老师的建议，CP2K 2024版本可以成功跑出激发态动力学。

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

[CP2K] 跑绝热的MD，追踪state 1的计算时，出现了Cholesky decompose failed