计算化学公社

标题: 为什么高斯结构优化后，再用相同方法基组算单点，能量差距这么大呢？ [打印本页]

作者
Author: abu123 时间: 2024-12-26 10:15
标题: 为什么高斯结构优化后，再用相同方法基组算单点，能量差距这么大呢？
各位老师好，我将一个结构进行结构优化后，将它的log转成gjf后，用相同方法和基组算了一下单点，发现能量差距特别大，我想问一下，这是为什么呢？第一个文件是结构优化的log文件，第二个结构是计算单点的log文件。

作者
Author: zjxitcc 时间: 2024-12-26 13:49
应当在帖子内容里主动写出哪个能量与哪个能量差距特别大，计算级别是什么，你的期望是什么。不能等回答者下载了文件再去找（他找的数值未必是你看的数值）

作者
Author: Freeman 时间: 2024-12-26 15:21
两次任务的自旋布居不一样，说明没有收敛到同一个电子态。单点的能量比优化的高，说明单点的波函数肯定是错的，但不保证优化的波函数是对的。我建议用stable=opt找到稳定的波函数，再用opt=calcfc重新优化。优化完成后，如果你想再用同样的理论水平算个单点，别忘了读取优化任务的波函数文件。

作者
Author: Loading0760 时间: 2024-12-26 15:31
自旋布局差的太大了。导致体系的能量差了很多，说明两次计算收敛到的波函数不是同一个，而且应该是单点计算的波函数有问题。处理方法楼上说的很详细了，我来解释一下。
首先对要优化的结构使用

pbe1pbe def2tzvpp stable=opt

复制代码

得到一个稳定的波函数
然后

guess=read opt=loose def2tzvpp pbe1pbe symm=loose

复制代码

结构优化
得到的结构这个时候再计算单点，波函数还是有可能会收敛到不稳定，此时可以读取结构优化的波函数计算。

作者
Author: lin-li 时间: 2026-1-14 22:54

Loading0760 发表于 2024-12-26 15:31
自旋布局差的太大了。导致体系的能量差了很多，说明两次计算收敛到的波函数不是同一个，而且应该是单点计算 ...

请问，如果优化使用B3LYP-D3/6-31G* 单点在（1）B3LYP-D3/6-31+G** 及（2）B3LYP-D3/def2TZVP 级别下，（1）计算的单点能量正常，但（2）计算的能量相差很多，大约十几至几十kcal/mol，这种情况如何处理呢？重新用大些的基组做优化吗，或者opt时加tight 可以吗？

作者
Author: sgwzq 时间: 7 day ago

lin-li 发表于 2026-1-14 22:54
请问，如果优化使用B3LYP-D3/6-31G* 单点在（1）B3LYP-D3/6-31+G** 及（2）B3LYP-D3/def2TZVP 级别下 ...

这个不是优化的问题。优化调的是几何构型，你这个明显是波函数自身的问题。有可能是6-31G*没法很好地描述电子分布（给的原子基函数太少了）也有可能是SCF没收敛到最低能量解。

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