计算化学公社

标题: 不同计算级别下的Gaussian波函数稳定性分析是否会得到同一个态? [打印本页]

作者
Author:
yueyue1765    时间: 2024-3-30 15:49
标题: 不同计算级别下的Gaussian波函数稳定性分析是否会得到同一个态?
各位老师好!我在用Gaussian不同的泛函计算同一个过渡金属体系的单点能,并做稳定性分析。我的疑惑是,理论上稳定波函数的定义是:在轨道正交归一的限制下,能量对轨道系数的一阶导数等于零,能量对轨道系数的Hessian矩阵本征值全为正。那么,对于不同的泛函,stable=opt的结果是否会得到同一个电子态的能量?如果换用不同的基组,又能否得到同一个电子态?
另外,在Gaussian的log文件中,是否有办法判断【同一体系,不同泛函,同一基组】的两次计算收敛得到的是否是同一个电子态?
谢谢!

作者
Author:
sobereva    时间: 2024-3-30 16:09
拿闭壳层单重态体系来说,HF成份较高的泛函更倾向于对称破缺(自旋极化单重态)。这在论坛里之前有很多讨论,可以用首页google框搜索。比如1,3-丁二烯在HF下做波函数稳定性测试都会提示不稳定,而主流的DFT泛函就没这个问题。
基组影响也有,但没理论方法那么显著
作者
Author:
wzkchem5    时间: 2024-3-30 16:56
甚至同一个泛函、同一个基组的stable=opt结果都未必收敛到同一个电子态,因为给定体系的稳定波函数未必是唯一的。典型例子如两个Fe(III)反铁磁耦合的体系,一般而言,两个Fe都是高自旋的态,和两个Fe都是低自旋的态,都是稳定的,但两者不等价




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