过渡态搜索stable检查波函数稳定性问题

西灵壹

我使用g16进行过渡态搜索（#p opt=(calcfc,noeigen,TS) freq m062x/6-311g**），一直报错l502，由于是SCF不收敛，参照http://sobereva.com/61里的内容进行修改，最终只有添加scf=qc才解决问题。但按教程说，需要使用stable进行波函数稳定性检查。那我之后的操作应该是基于加了scf=qc跑完的fch格式和坐标，采用m062x/6-311g** stable guess=read，进行计算（这里的输入文件的坐标用收敛后的最终坐标）。如果波函数稳定性没问题以后，我后续用更高精度CCSD(T)/cc-pVTZ计算单点能，需要读取stable检验后的chk文件么，还是只要使用过渡态最终结构就行？

zjxitcc

问题描述不合适，无法直接回复“是/否/可以”。你先说说输入文件中自旋多重度是什么，你好不容易收敛的波函数对应的<S**2>值是多少

sobereva

CCSD(T)是用HF当参考态的，哪怕你读取了M06-2X的收敛而且判断为稳定的波函数当初猜，也不代表收敛到的HF的波函数就是稳定的。因此你在做CCSD(T)之前最好再在当前结构下做HF的stable任务并且guess=read M06-2X的波函数。如果此时提示波函数也稳定，之后做CCSD(T)的时候guess=read读取这个HF的stable任务的chk文件里的波函数当初猜。

Xx_DerEnTh

sobereva 发表于 2024-6-14 17:02
CCSD(T)是用HF当参考态的，哪怕你读取了M06-2X的收敛而且判断为稳定的波函数当初猜，也不代表收敛到的HF的 ...

sob老师您好！我准备用CCSD(T)算一个结构的单点，我先在这个结构下做了HF计算，并且用stable=opt把HF的波函数弄稳定了，随后我读取这个稳定的HF波函数做CCSD(T)计算，计算级别、结构均没有改变，但是显示SCF用了8步才收敛，这种情况正常吗？这样会不会明显改变之前已经收敛且稳定的HF波函数，有没有可能波函数变得不稳定了？

wzkchem5

Xx_DerEnTh 发表于 2025-6-18 13:16
sob老师您好！我准备用CCSD(T)算一个结构的单点，我先在这个结构下做了HF计算，并且用stable=opt把HF的波 ...

第一步SCF的能量误差是多少？

Xx_DerEnTh

wzkchem5 发表于 2025-6-18 13:17
第一步SCF的能量误差是多少？

0.000000008

wzkchem5

Xx_DerEnTh 发表于 2025-6-18 13:23
0.000000008

最后一步收敛的能量和第一步是否非常接近？如果是，那么问题不大。
判断SCF收敛往往是靠SCF迭代一步的能量、密度矩阵变化多少。但是当打开DIIS的时候，迭代一步所导致的能量、密度矩阵的变化，是和迭代历史有关的。迭代15次能收敛，不代表从最后一步开始重新迭代能1步收敛。但如果是基于轨道梯度判断收敛性的程序，就没有这个问题。

Xx_DerEnTh

wzkchem5 发表于 2025-6-18 13:55
最后一步收敛的能量和第一步是否非常接近？如果是，那么问题不大。
判断SCF收敛往往是靠SCF迭代一步的能 ...

明白了，感谢老师！