过渡态搜索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文件里的波函数当初猜。