Gaussian计算过渡金属配合物不同自旋态能量差的问题

yueyue1765

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

zjxitcc

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

yueyue1765

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

sobereva

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

后者

yueyue1765

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

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

sobereva

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

能量更低

zjxitcc

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组态。

yueyue1765

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

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

zjxitcc

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

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