可以尝试用wB2GP-PLYP或RSX-QIDH看看,本身其杂化部分就是范围分离泛函 或者用w调控的范围分离泛函,更稳 |
清明雨上 发表于 2024-10-22 15:47 1eV已经比较大了,尤其在激发能本身都不到1eV的情况下,更是相当大了。一个校正和被校正的量本身差不多,这个在大部分情况下都说明这个校正很大 |
不能轻易下这样的结论。 可以在相同结构下 用普通的范围分离泛函和常见的多种双杂化都算单点横向比较一下。也许确实DSD-PBEP86对于你的体系非常不合适。 |
wzkchem5 发表于 2024-10-22 15:54 非常感谢老师教的方法,手册里的同义词为(D)correction。在输出文件中搜出来的结果为 Contributions to the (D) correction(s) (in eV) STATE Final Init. Total v(ext) v(int) v(3) u omega omega (D) ----------------------------------------------------------------------------- 0 0.387 1.460 -1.073 1.207 2.100 -0.051 -4.329 1 1.656 3.442 -1.785 1.248 2.081 -0.041 -5.073 2 1.359 3.491 -2.132 1.242 2.177 -0.034 -5.517 发现第一个态的doubles Correction没有太大,私以为可能是范围分离泛函激发态结构优化不合理导致双杂化泛函单点计算不合理。 |
清明雨上 发表于 2024-10-22 08:50 应该在手册里搜。 以后遇到别人叫你搜某关键词A,但是A搜不到的情况,应当先去了解A都有哪些同义词,然后在输出文件里依次搜这些同义词。了解A有哪些同义词的方法之一就是在手册里搜A,然后仔细阅读相应章节 |
感谢老师的回复,由于我对Doubles Correction是什么缺乏了解,在输出文件中搜索Doubles Correction,并未发现直接的数值,只注意到CI dimension:23715 同时看了一下total energy的组成,components和Virial components和范围分离泛函计算出来的几乎一致,而DSD-PBEP86/def2-TZVP算出来的E(XC)=-106与cam-b3lyp/def2-TZVP算出来的E(XC)=-286相差较大,符合预期。 最后检查DSD-PBEP86/def2-TZVP的out文件发现了 CIS/TD-DFT TOTAL ENERGY ----------------------- E(SCF) = -5462.034020634 Eh DE(CIS) = 0.014218971 Eh (Root 1) ----------------------------- --------- E(tot) = -5462.019801663 Eh 不知道和您说的doubles correction有没有关联 |
可以检查一下第一个态的doubles correction是不是太大。双杂化泛函容易有doubles correction发散的问题,虽然大部分情况下比普通杂化泛函、范围分离泛函准,但是少部分情况下误差可以很大 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2024-11-27 06:45 , Processed in 0.240362 second(s), 26 queries , Gzip On.