新手关于四点法计算重组能求助

科研小白0126

各位老师好，我化学和理论计算基础较为薄弱，对重组能这里的一些步骤不太理解，希望各位老师能给予指导。
1.这里面E1是中性态分子能量是out文件里的SCF done吗
2.关于电荷和自旋多重度更改，若我将原来的0  1改为阳离子是不是+1   2，阴离子是不是-1   2，这个更改不是很理解
3.E2是指我直接改变优化后中性态分子的电荷和自旋多重度，然后直接计算单点能吗，然后在优化改过的电荷的分子结构，再计算对应单点能E3、E4吗
可能我的问题有的比较多余，还是希望有老师能予以解答，感谢

sgwzq

如果你求的是电子能量，那么读SCF Done，
如果你想算热力学数据（焓变/吉布斯自由能变）这个需要做振动频率分析。（当然他应该不是这个意思）

原来的体系的自旋多重度为1, 加入或减少一个电子后，没有与这个电子自旋相反的匹配
（相当于这个轨道只有一个电子）所以自旋多重度为2.

这几个能量，解释起来比较费劲，或许你看看sober老师博客里面的电子跃迁过程的图片就会理解了
http://sobereva.com/314

Gaussian中用TDDFT计算激发态和吸收、荧光、磷光光谱的方法。

多说一嘴热力学数据（因为我自己在这里也想了好久）
几何优化是指将结构通过一些算法（比如高数讲过的牛顿切线法），将分子中各原子放在附近受力为0、势能极小的位置。
这个位置是势能面的极小点，一阶导数（负梯度即受力）为零，二阶导数为正数，形状类似于开口向正方向的抛物线。
在0K下，原子核不会处于静态，否则位置坐标同时确定，违反不确定原理，所以它一定有一些微小振动。
这些振动可以用一些方法求
（比如把势能，想象成有一个拉着原子核的弹簧，弹簧的劲度系数等于势能对原子坐标的二阶偏导数。
弹簧拉着原子核往复运动有一个固定的频率，这个叫作谐振频率，虽然有些时候不能简单地这么考虑）
求出这些振动的频率后，可以确定分子体系在0K时具有的能量，即零点能Zero Point Energy，
如果再考虑298.15K分子振动构象所具有的能量，就能确定分子的内能，
如果再加上pV即再加上nRT(此处用到理想气体近似)，就能确定分子焓，
如果再减去TS，就能得到吉布斯自由能。（S的算法我也不太清楚，可能是先求配分函数）

这几个数据都会在提交振动频率分析任务的时候，同步输出。找一下就能找得到。

sobereva

1 认真把此文看了
谈谈该从Gaussian输出文件中的什么地方读电子能量
http://sobereva.com/488（http://bbs.keinsci.com/thread-13450-1-1.html）

2 是。这都是最基本的量子化学常识，这都不理解的话强烈建议从头系统性学学
谈谈学量子化学如何正确地入门
http://sobereva.com/355（http://bbs.keinsci.com/thread-4447-1-1.html）

3 是

科研小白0126

sgwzq 发表于 2026-1-9 00:04
如果你求的是电子能量，那么读SCF Done，
如果你想算热力学数据（焓变/吉布斯自由能变）这个需要做振动频 ...

太感谢您了

科研小白0126

sobereva 发表于 2026-1-9 00:19
1 认真把此文看了
谈谈该从Gaussian输出文件中的什么地方读电子能量
http://sobereva.com/488（http://bb ...

感谢老师指导