计算化学公社

标题: Gaussian计算过渡金属配合物不同自旋态能量差的问题 [打印本页]

作者
Author: yueyue1765 时间: 2023-11-29 17:50
标题: Gaussian计算过渡金属配合物不同自旋态能量差的问题
各位老师好！我在用Gaussian计算过渡金属配合物体系的两种自旋态（高自旋和低自旋）的能量差时，碰到了一些问题：
文献中用的方法是，对于单重态体系，使用RDFT；对于多重态体系，使用UDFT。然而，用这种方法做波函数稳定性分析后，单重态体系大量出现RHF->UHF instability的提示。如果用stable=opt，最终得到的UDFT能量和原先的RDFT能量相差很多。
我感觉这种处理方法有点问题，但不确定怎么做更合适。请问我是否应该统一用UDFT计算，再加上stable=opt？

作者
Author: zjxitcc 时间: 2023-11-29 20:33
本帖最后由 zjxitcc 于 2023-11-29 20:42 编辑

你是应该加上stable=opt，但单重态UDFT nosymm guess=mix stable=opt几个应该一起写。如果你少写了guess=mix，会先用UDFT代码算出与RDFT一样的结果，浪费一次机时；随后stable=opt从闭壳层优化至开壳层，这个过程往往比较坎坷，又浪费一次机时。
另外，在发表文章时应在SI里给出每个结构的UDFT能量和<S**2>的值。

作者
Author: yueyue1765 时间: 2023-12-1 00:35
谢谢@zjxitcc 老师！已经在用这种方法计算了。
还想请教一个问题，我在算某些体系的时候，用scf=(novaracc,noincfock,vshift=500)后仍然没有收敛，能量出现振荡的情况。我分别用两种方法处理：
1. 改用二次收敛方法scf=qc
2. 放宽收敛限scf=conver=6（配合IOp(8/117=-99)）
这两者最终都能收敛，前者稳定性分析时直接得到稳定的波函数，后者程序判断internal stability，通过stable=opt也得到稳定波函数。但是，前者的能量比后者高了0.0004Hartree左右。我不确定哪种方法是更合适的，或者有没有更好的方法。

作者
Author: sobereva 时间: 2023-12-1 06:22

yueyue1765 发表于 2023-12-1 00:35
谢谢@zjxitcc 老师！已经在用这种方法计算了。
还想请教一个问题，我在算某些体系的时候，用scf=(novaracc ...

后者

作者
Author: yueyue1765 时间: 2023-12-1 10:05

sobereva 发表于 2023-12-1 06:22
后者

谢谢sob老师！请问选择后者的原因是什么呢？是因为qc方法应该迫不得已才使用，还是说单纯因为这个体系中后者能得到能量更低的结果？

作者
Author: sobereva 时间: 2023-12-2 11:00

yueyue1765 发表于 2023-12-1 10:05
谢谢sob老师！请问选择后者的原因是什么呢？是因为qc方法应该迫不得已才使用，还是说单纯因为这个体系中 ...

能量更低

作者
Author: zjxitcc 时间: 2023-12-3 11:05

yueyue1765 发表于 2023-12-1 00:35
谢谢@zjxitcc 老师！已经在用这种方法计算了。
还想请教一个问题，我在算某些体系的时候，用scf=(novaracc ...

qc.com文件里不应该加guess=mix。这个文件目前得到的是稳定波函数，但没有收敛到Fe的3d^6组态上，具体是什么组态，需要看UNO轨道判断，不过这倒也不是很重要，因为能量比另一个任务的能量高。

conver.com这个任务最好读取当前轨道继续算至正常收敛阈值，结果应是-1602.889776 a.u.，可通过观看UNO轨道图形和轨道占据数判断Fe是3d^6组态。

作者
Author: yueyue1765 时间: 2023-12-4 14:04

zjxitcc 发表于 2023-12-3 11:05
qc.com文件里不应该加guess=mix。这个文件目前得到的是稳定波函数，但没有收敛到Fe的3d^6组态上，具体是 ...

谢谢老师！还想问一下不加guess=mix的原因是什么？是所有用scf=qc的都不应该加guess=mix，还是说多重态的体系不用加guess=mix？
另外，可以直接用gaussview读取fchk文件来判断有没有收敛到我想要的态吗？

作者
Author: zjxitcc 时间: 2023-12-4 14:07

yueyue1765 发表于 2023-12-4 14:04
谢谢老师！还想问一下不加guess=mix的原因是什么？是所有用scf=qc的都不应该加guess=mix，还是说多重态的 ...

guess=mix仅给对称破缺单重态计算使用，其他情况均不使用。
我就是用GaussView打开UNO的fch文件看轨道和轨道占据数判断的（前提是要产生UNO的fch文件）。

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