四点法计算重组能E2小于E4是为什么？

身在盘丝洞 · 发表于 Post on 2025-4-25 16:06:23

本帖最后由身在盘丝洞于 2025-4-25 17:51 编辑

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

wcy鱼饵 · 发表于 Post on 2025-4-25 16:48:08

泛函不一致吧，因为你E2用的M062x E4用的BELYP。你统一一下泛函比如如果都用M062X E4优化完结构用M062X算个能量

身在盘丝洞 · 发表于 Post on 2025-4-25 16:50:40

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

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

wcy鱼饵 · 发表于 Post on 2025-4-25 16:57:18

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

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

wal · 发表于 Post on 2025-4-25 17:08:56

简单来说，B3LYP/6-311g(d,p)势能面上S0/S1极小点结构未必对应M062X/6-311g(d,p)势能面上的S0/S1极小点结构所以目前你算的点不对应图上的E2/E4

身在盘丝洞 · 发表于 Post on 2025-4-25 17:21:38

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

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

陈AG · 发表于 Post on 2025-4-27 15:57:48

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

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

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[综合交流] 四点法计算重组能E2小于E4是为什么？

评分 Rate

浏览过的版块