计算化学公社

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

我应该如何解决？能否在气相下 opt freq，然后加上溶剂算单点能？
(, 下载次数 Times of downloads: 3)



作者
Author: wal    时间: 2025-7-2 16:23
要不要先带溶剂几何优化一下看看它会优化成啥样
作者
Author: 北大-陶豫    时间: 2025-7-2 16:34
wal 发表于 2025-7-2 16:23
要不要先带溶剂几何优化一下看看它会优化成啥样

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

有点神奇 stable一下试试?
作者
Author: 北大-陶豫    时间: 2025-7-2 21:09
wal 发表于 2025-7-2 16:38
有点神奇 stable一下试试?

谢谢！刚算好，结果也挺神奇，如下所示：
  1. ***********************************************************************
  2. Stability analysis using <AA,BB:AA,BB> singles matrix:
  3. ***********************************************************************
  4. Ground to excited state transition densities written to RWF  633

  6. Eigenvectors of the stability matrix:

  8. Eigenvector   1:      Triplet-A    Eigenvalue= 0.0414753  <S**2>=2.000
  9.      154 ->161        -0.13406
  10.      156 ->161         0.56910
  11.      157 ->160        -0.18168
  12. The wavefunction is stable under the perturbations considered.
  13. The wavefunction is already stable.

  15. **********************************************************************
复制代码


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

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

stable之后继承波函数可以收敛？这样似乎波函数是稳定的
作者
Author: wal    时间: 2025-7-2 21:17
本帖最后由 wal 于 2025-7-2 21:28 编辑
北大-陶豫 发表于 2025-7-2 21:09
谢谢！刚算好，结果也挺神奇，如下所示：

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


编辑：是au，算了个水的是这样
(, 下载次数 Times of downloads: 60)

那看来不像波函数不稳定
作者
Author: 北大-陶豫    时间: 2025-7-2 21:33
本帖最后由 北大-陶豫 于 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
(, 下载次数 Times of downloads: 0)





欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3