计算化学公社

标题: 求助，计算时使用MD指定的轨迹作为TDDFT的输入文件的坐标，出现了报错 [打印本页]

作者
Author: 不会扣篮的后卫 时间: 2023-12-8 11:50
标题: 求助，计算时使用MD指定的轨迹作为TDDFT的输入文件的坐标，出现了报错
大家好，我最近在研究使用非绝热动力学算孤立体系的金刚烷（Adamantane, C10H16）的量子产率。
大致流程是通过CP2K和Libra（一个非绝热动力学的开源软件）结合来算。
在走非绝热动力学流程前，我先按照卢老师CP2K培训班讲义的例子，跑了些孤立体系的C10H16，单电子激发的例子，来决定选择合适的泛函和基组。
(, 下载次数 Times of downloads: 0)

黄色背景的例子是先用B3LYP/6-311G**优化，然后再用不同泛函，基组做TDDFT计算。
蓝色背景的例子是先用PBE0/6-311G**优化，然后再用不同泛函，基组做TDDFT计算。
比较下来，我感觉泛函B3LYP，基组6-311G**算出来激发能的结果和实验值差不多，同时，计算时间相对较快。就打算把这个组合应用到后面的计算。

第一步：跑MD计算任务(泛函B3LYP，基组6-311++G**)。
第二步，把MD计算任务的输出轨迹，作为TDDFT计算（泛函B3LYP，基组6-311++G**）的任务输入文件的坐标，来进行激发态的计算，然后跑出time overlap。

但是第二步的任务出现了很多这样的报错， *** WARNING in pw/ps_wavelet_methods.F:238 :: Density non-zero on the ***
*** edges of the unit cell: wrong results in WAVELET solver ***

我后续尝试第一步：跑MD计算任务(泛函PBE0，基组6-311++G**)。
然后第二步的TDDFT计算分别用泛函PBE0，基组6-311++G**和 CAM-B3LYP, 基组6-311++G**,还是会出现上述的报错。
我觉得出现这个报错的原因的是MD的轨迹坐标相对于优化后结构，没有那么稳定。

请问大家这个报错有什么建议吗（附上第一步的轨迹文件和第二步的输出文化供大家参考）？

作者
Author: sobereva 时间: 2023-12-8 12:05
用wavelet时，电子密度不能再在盒子边界有分布，否则结果不靠谱，并且会有这种提示。出现这种提示，要么盒子尺寸太小而电子分布过于弥散（增大盒子解决），要么当前波函数诡异，出现了虚假且弥散的特征，需要确保波函数收敛情况合理、计算方式没有硬伤，并且如果只是SCF初期出现而后期自动消除则可以不用管

作者
Author: 不会扣篮的后卫 时间: 2023-12-10 03:15

sobereva 发表于 2023-12-8 12:05
用wavelet时，电子密度不能再在盒子边界有分布，否则结果不靠谱，并且会有这种提示。出现这种提示，要么盒 ...

谢谢卢老师！

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