计算化学公社

标题: 高斯计算三重态的自旋不对，如何解决 [打印本页]

作者
Author: LYFE 时间: 2021-11-5 09:54
标题: 高斯计算三重态的自旋不对，如何解决
如图，高斯计算后的log文件，查找spin后，Cu的自旋为一正一负，我计算的是三重态，按道理来说应该是两个接近0.5的正数，请问各位大佬们，有什么解决办法
(, 下载次数 Times of downloads: 23) (, 下载次数 Times of downloads: 20)

作者
Author: zjxitcc 时间: 2021-11-5 10:00
本帖最后由 zjxitcc 于 2021-11-5 10:12 编辑

（1）scfcon=5算出来结果没有实际意义，不能用于发表，作为粗结果、后续再恢复正常精度计算是可以的，但直接拿来发表的话是严重的垃圾数据。不是说程序不报错、结果就能用，这是错误的观点。
（2）复杂结构应该先算单点，检验波函数稳定性（stable=opt），检查自旋布居与期望的是否有偏差。然后再读取波函数做结构优化。如果波函数不稳定，布居分析是没有意义的。

由于你已经做了一次结构优化，不妨就在此结构上（记得恢复正常收敛精度）算单点、检验波函数稳定性，看看结果如何。

作者
Author: LYFE 时间: 2021-11-5 10:11

zjxitcc 发表于 2021-11-5 10:00
（1）scfcon=5算出来结果没有实际意义，不能用于发表，作为粗结果、后续再恢复正常精度计算是可以的，但直 ...

老师，请问若自旋与期望值有偏差，以什么办法解决呢

作者
Author: zjxitcc 时间: 2021-11-5 10:13

LYFE 发表于 2021-11-5 10:11
老师，请问若自旋与期望值有偏差，以什么办法解决呢

可能性：
（1）自己的化学直觉有误，猜想有误；
（2）自己的化学直觉无误，SCF存在多个解，甚至多个稳定的解，需要找出其他的解。除了确保波函数稳定外，还可以借助《谈谈片段组合波函数与自旋极化单重态》http://sobereva.com/82 找出其他的SCF解

作者
Author: LYFE 时间: 2021-11-5 10:16

zjxitcc 发表于 2021-11-5 10:13
可能性：
（1）自己的化学直觉有误，猜想有误；
（2）自己的化学直觉无误，SCF存在多个解，甚至多个稳 ...

老师，还有一个问题，高斯计算中scfcon是代表精度嘛，scfcon等于的数值越高越好嘛

作者
Author: zjxitcc 时间: 2021-11-5 10:21

LYFE 发表于 2021-11-5 10:16
老师，还有一个问题，高斯计算中scfcon是代表精度嘛，scfcon等于的数值越高越好嘛

高斯官网都有http://gaussian.com/scf，相当于scf(conver=N)，N越大越精确，太大没必要，徒增耗时，但不要比默认值粗糙。

作者
Author: LYFE 时间: 2021-11-5 10:22

zjxitcc 发表于 2021-11-5 10:21
高斯官网都有http://gaussian.com/scf，相当于scf(conver=N)，N越大越精确，太大没必要，徒增耗时，但不 ...

好的，感谢老师！！！

作者
Author: sobereva 时间: 2021-11-5 16:54
你都不给出体系具体结构，这个问题根本没法回答，别人都不知道两个Cu的化学环境是否相同、位置是否是对称的
但凡用gen的时候，具体用的什么基组都必须说清楚，下文明确强调了。仔细看下文再提问
在网上求助计算化学问题的时候必须把问题描述得详细、具体、准确、清楚
http://sobereva.com/620（http://bbs.keinsci.com/thread-25787-1-1.html）

不懂的关键词一律不要随便加，scf=conver=x（等同于scfcon=x）怎么用好好看下文
解决SCF不收敛问题的方法
http://sobereva.com/61

写iop(1/8=6)完全莫名其妙。有maxstep现成的关键词不用用这么抽象的IOp做什么，而且不加notrust的话根本起不到实际作用。maxstep怎么用、什么时候才该用仔细看
量子化学计算中帮助几何优化收敛的常用方法
http://sobereva.com/164

绝对不要从别人那里瞎抄关键词。我担心你在入门阶段看了不少乱七八糟的东西已被误导，建议看
谈谈学量子化学如何正确地入门
http://sobereva.com/355（http://bbs.keinsci.com/thread-4447-1-1.html）

作者
Author: LYFE 时间: 2021-11-7 23:18

sobereva 发表于 2021-11-5 16:54
你都不给出体系具体结构，这个问题根本没法回答，别人都不知道两个Cu的化学环境是否相同、位置是否是对称的 ...

好的，我是个新手，才学习，谢谢老师的教导！！！

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