高斯计算受力，同一结构加溶剂化后受力增大1945倍

北大-陶豫 · 发表于 Post on 2025-7-2 16:18:30

最近我同学优化一个钒配合物时发现总是难以收敛，我帮他debug时发现，同一结构加溶剂化后，最大受力会增大1945倍（溶剂化前0.005142，溶剂化后9.999381），似乎溶剂化之后得到的受力很不合理。
看了一下输出文件中各个原子的受力，尤其关注受力巨大的原子，但它们的几何结构似乎并没有不合理之处。
由于已经在同样的理论水平下计算了很多其它结构的 opt freq，因此更换方法推倒重来可能成本过大。

我应该如何解决？能否在气相下 opt freq，然后加上溶剂算单点能？

VO-L1-Substrate-NVOOH-OOH-force.zip (40.97 KB, 下载次数 Times of downloads: 3)

wal · 发表于 Post on 2025-7-2 16:23:59

要不要先带溶剂几何优化一下看看它会优化成啥样

北大-陶豫 · 发表于 Post on 2025-7-2 16:34:37

wal 发表于 2025-7-2 16:23
要不要先带溶剂几何优化一下看看它会优化成啥样

优化不收敛，只能走到比该结构能量还高的结构，或者走几步就SCF不收敛。。

wal · 发表于 Post on 2025-7-2 16:38:07

北大-陶豫发表于 2025-7-2 16:34
优化不收敛，只能走到比该结构能量还高的结构，或者走几步就SCF不收敛。。

有点神奇 stable一下试试?

北大-陶豫 · 发表于 Post on 2025-7-2 21:09:16

wal 发表于 2025-7-2 16:38
有点神奇 stable一下试试?

谢谢！刚算好，结果也挺神奇，如下所示：

***********************************************************************
Stability analysis using <AA,BB:AA,BB> singles matrix:
***********************************************************************
Ground to excited state transition densities written to RWF 633
Eigenvectors of the stability matrix:
Eigenvector 1: Triplet-A Eigenvalue= 0.0414753 <S**2>=2.000
154 ->161 -0.13406
156 ->161 0.56910
157 ->160 -0.18168
The wavefunction is stable under the perturbations considered.
The wavefunction is already stable.
**********************************************************************

复制代码

这是说激发能只有 0.04 eV 吗？

laocuile · 发表于 Post on 2025-7-2 21:16:11

北大-陶豫发表于 2025-7-2 21:09
谢谢！刚算好，结果也挺神奇，如下所示：

stable之后继承波函数可以收敛？这样似乎波函数是稳定的

wal · 发表于 Post on 2025-7-2 21:17:52

本帖最后由 wal 于 2025-7-2 21:28 编辑

北大-陶豫发表于 2025-7-2 21:09
谢谢！刚算好，结果也挺神奇，如下所示：

好像是这样看会不会是带溶剂走几步之后这个三重态掉到单重态下面了 XD 然后就波函数不稳定，不收敛了
不对这个单位是eV还是au啊

编辑：是au，算了个水的是这样

那看来不像波函数不稳定

北大-陶豫 · 发表于 Post on 2025-7-2 21:33:51

本帖最后由北大-陶豫于 2025-7-2 21:39 编辑

wal 发表于 2025-7-2 21:17
好像是这样看会不会是带溶剂走几步之后这个三重态掉到单重态下面了 XD 然后就波函数不稳定，不收敛了
...

补充一下：刚用Multiwfn看了一下fchk文件里的波函数，HOMO 是 -0.219020 au，-5.9598 eV，是硫的孤对；LUMO 是 -0.121753 au，-3.3131 eV，是钒的 t2g 轨道。
HOMO-LUMO gap: 0.097267 a.u. 2.646768 eV 255.374330 kJ/mol

VO-L1-Substrate-NVOOH-OOH-force.zip (40.97 KB, 下载次数 Times of downloads: 0)

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

[Gaussian/gview] 高斯计算受力，同一结构加溶剂化后受力增大1945倍

评分 Rate

浏览过的版块