计算化学公社

标题: 跑绝热的MD,追踪state 1的计算时,出现了Cholesky decompose failed [打印本页]

作者
Author:
不会扣篮的后卫    时间: 2024-3-4 12:54
标题: 跑绝热的MD,追踪state 1的计算时,出现了Cholesky decompose failed
大家好。如题,我在用孤立体系,跑绝热的NVT MD,追踪激发态state 1的计算时,出现了这个报错
(, 下载次数 Times of downloads: 18)
请问大家有什么修改建议吗?目前的打算是跑出这个计算的轨迹文件,用于后续的计算(附上输入,输出文件供大家参考)。


作者
Author:
sobereva    时间: 2024-3-5 02:27
B3LYP对里德堡激发描述很烂,换成wB97X或者全局高HF成份的BHandHLYP
尝试把盒子设得更大。里德堡激发会使得电子高度弥散
作者
Author:
不会扣篮的后卫    时间: 2024-3-5 12:01
sobereva 发表于 2024-3-5 02:27
B3LYP对里德堡激发描述很烂,换成wB97X或者全局高HF成份的BHandHLYP
尝试把盒子设得更大。里德堡激发会使 ...

谢谢卢老师!
作者
Author:
不会扣篮的后卫    时间: 2024-4-15 12:54
sobereva 发表于 2024-3-5 02:27
B3LYP对里德堡激发描述很烂,换成wB97X或者全局高HF成份的BHandHLYP
尝试把盒子设得更大。里德堡激发会使 ...

卢老师,您好。我按照您的建议,增大了盒子,外加使用了wB97x泛函,但是出现了新的报错。
(, 下载次数 Times of downloads: 17)

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


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

分清楚warning和error
尝试减小EPS_SCF、TDDFT的CONVERGENCE、EPS_DEFAULT、增大CUTOFF
如果还这样,尝试诸如CAM-B3LYP,如果还这样,暂且无视,继续观望
作者
Author:
不会扣篮的后卫    时间: 2024-4-17 12:43
本帖最后由 不会扣篮的后卫 于 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文件的报错

(, 下载次数 Times of downloads: 17)


还有文件out.log的warning信息以及振子强度为0的情况,
(, 下载次数 Times of downloads: 20)
因为这个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轨迹有区别吗?







作者
Author:
sobereva    时间: 2024-4-18 07:13
用CP2K最新版本再试

实在搞不定就换成ORCA跑TDDFT的AIMD
作者
Author:
不会扣篮的后卫    时间: 2024-4-19 22:02
sobereva 发表于 2024-4-18 07:13
用CP2K最新版本再试

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

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




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3