计算化学公社

标题: 四点法计算重组能E2小于E4是为什么？ [打印本页]

作者
Author: 身在盘丝洞 时间: 2025-4-25 16:06
标题: 四点法计算重组能E2小于E4是为什么？
本帖最后由身在盘丝洞于 2025-4-25 17:51 编辑

我先使用dushin分解了重组能，想在试一下四点法计算，我直接在分解重组能的结构上计算了单点能，以total energy读数发现E2小于E4，HF为相应基态的能量。
optfreqs0和optfreqs1为分解重组能用到的s0和s1的结构和频率优化，在此基础上做了相应的单点计算es0和es1。

作者
Author: wcy鱼饵 时间: 2025-4-25 16:48
泛函不一致吧，因为你E2用的M062x E4用的BELYP。你统一一下泛函比如如果都用M062X E4优化完结构用M062X算个能量

作者
Author: 身在盘丝洞 时间: 2025-4-25 16:50

wcy鱼饵发表于 2025-4-25 16:48
泛函不一致吧，因为你E2用的M062x E4用的BELYP。你统一一下泛函比如如果都用M062X E4优化完结构用M062X算 ...

单点都是M062X，哪里用了b3lyp?

作者
Author: wcy鱼饵 时间: 2025-4-25 16:57

身在盘丝洞发表于 2025-4-25 16:50
单点都是M062X，哪里用了b3lyp?

您基态结构和激发态结构优化的时候用的是B3LYP，我如果没看错的话应该使用的log结构去进行的M062X的单点计算，你如果用的是优化结构的HF值为E4 那么计算的单点能量应该就是B3LYP泛函下的能量

作者
Author: wal 时间: 2025-4-25 17:08
简单来说，B3LYP/6-311g(d,p)势能面上S0/S1极小点结构未必对应M062X/6-311g(d,p)势能面上的S0/S1极小点结构所以目前你算的点不对应图上的E2/E4

作者
Author: 身在盘丝洞 时间: 2025-4-25 17:21

wcy鱼饵发表于 2025-4-25 16:57
您基态结构和激发态结构优化的时候用的是B3LYP，我如果没看错的话应该使用的log结构去进行的M062X的单点 ...

我没说用优化结构的HF值为E4

作者
Author: 陈AG 时间: 2025-4-27 15:57

身在盘丝洞发表于 2025-4-25 17:21
我没说用优化结构的HF值为E4

你可以列一下，E1-E4分别来自哪个文件

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