计算化学公社

标题: 求助CP2K计算的gap问题 [打印本页]

作者
Author:
pet    时间: 2024-3-8 14:25
标题: 求助CP2K计算的gap问题
本帖最后由 pet 于 2024-3-8 14:26 编辑

因为需要用一个软件计算di/dv map,他不能读高斯的输出,可以读cp2k,所以先用高斯b3lyp def2tzvp优化结构,然后用cp2k计算单点以产生所需波函数。cp2k分别用了PBE/DZVP-MOLOPT-SR-GTH, HSE06/DZVP-MOLOPT-SR-GTH, B3LYP/Ahlrichs-def2-TZVP,输入文件用Multiwfn产生,得到的gap分别是0.29,0.03,8.28 eV,PBE算出来的gap跟其他软件相差较小,但是HSE06和B3LYP算的gap与其他软件相差甚大,HSE06算出来的为啥比PBE还小呢,B3LYP就大的离谱,请各位老师帮我看看是输入文件有问题嘛? (, 下载次数 Times of downloads: 1) (, 下载次数 Times of downloads: 0) (, 下载次数 Times of downloads: 0)
作者
Author:
sobereva    时间: 2024-3-9 06:51
不要用糟糕的PSOLVER ANALYTIC,也绝对不要用过时版本的Multiwfn创建输入文件。Multiwfn对非周期性体系默认用的PSOLVER早已不是ANALYTIC

而且你用了SCREEN_ON_INITIAL_P T却没提供含有收敛的波函数当初猜,结果毫无意义,这是致命硬伤。仔细看下文了解CP2K做杂化泛函计算的基本常识
CP2K做杂化泛函计算的关键要点和简单例子
http://sobereva.com/690http://bbs.keinsci.com/thread-41087-1-1.html

CP2K切勿当黑箱用,完整系统学习一遍对于正确使用CP2K极其重要





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