stable=opt稳定后用stable关键字检测显示不稳定

Stepheee

最近在复现文献中的一个数据时遇到了奇怪的错误，计算的是一个五重态的FeO分子，在进行CCSD(T)计算之前先测试了一下HF波函数的稳定性，然后发现了如下的规律：

(1) stable=opt后再次用stable=opt会提示"The wavefunction is already stable."
(2) stable=opt后用stable会提示"The wavefunction has an internal instability."
(3) 第(2)步后再用stable=opt同样会提示instability然后优化到稳定。
(4) (1),(2),(3)步可以无限循环

因此我的问题如下：
1. stable=opt和stable检测稳定性采用的算法不一样吗，为什么stable=opt优化到稳定的波函数在两种关键字下再次检测出现了不同的结果？
2. guess=read和stable关键字会改变波函数吗？为什么经过(1)(2)(3)步骤后，原本stable=opt提示稳定的波函数在stable=opt关键字下也不稳定了？（在guess=read之后，各个步骤的scf过程都只有一步）
3. 对于当前这个体系来说，如何保证波函数的稳定性？


附输入文件：

(输入文件命名为FeO_5.gjf)

1. %nprocshared=32
   
2. %mem=40GB
   
3. %chk=FeO_5.chk
   
4. #p hf/6-311++g(3df,3pd) scf=fermi stable
   
5. 
   
6. FeO_5
   
7. 
   
8. 0 5
   
9.   Fe
   
10.   O 1 B1
    
11. 
    
12. B1=1.618
    
13. 
    
14. 
    

复制代码

提示"internal instability"，因此将上述输入文件改为stable=opt：

1. %nprocshared=32
   
2. %mem=40GB
   
3. %chk=FeO_5.chk
   
4. #p hf/6-311++g(3df,3pd) scf=fermi stable=opt guess=read
   
5. 
   
6. FeO_5
   
7. 
   
8. 0 5
   
9. Fe
   
10. O 1 B1
    
11. 
    
12. B1=1.618
    
13. 
    

复制代码

继续计算后提示"stable under perturbations"，此时如果将stabl=opt改为stable后再次检测稳定性：

1. %nprocshared=32
   
2. %mem=40GB
   
3. %chk=FeO_5.chk
   
4. #p hf/6-311++g(3df,3pd) scf=fermi stable guess=read
   
5. 
   
6. FeO_5
   
7. 
   
8. 0 5
   
9. Fe
   
10. O 1 B1
    
11. 
    
12. B1=1.618
    
13. 
    

复制代码

则会再次提示"internal instability"。

zjxitcc

CCSD(T)计算最好采用相关一致性基组，如cc-pVTZ基组。如果机器算力较强，可以用cc-pVQZ，高角动量极化函数远比6-311++g(3df,3pd)加的那些东西重要很多。stable或stable=opt最好搭配nosymm一起使用。读取已收敛的波函数时一般不需要加scf=fermi，因为如果读进来的波函数是稳定的，检验稳定性后计算就正常结束了，用不到scf=fermi；如果读进来的波函数不稳定，stable=opt会调用scf=qc做轨道优化，这个收敛比较强劲，压根不需要scf=fermi凑热闹。

hebrewsnabla

我本来想复现一下这个例子，但是我第一步用了scf(vshift=200,xqc)以后直接就stable了……

又用scf(fermi)算了一遍，还是stable

Stepheee

hebrewsnabla 发表于 2022-9-19 14:09
我本来想复现一下这个例子，但是我第一步用了scf(vshift=200,xqc)以后直接就stable了……

又用scf(fermi ...

老师好，这样就奇怪了，这个输入文件我在集群上复现了无数次，集群上安装的是gaussian16 C.01版本

zjxitcc

Stepheee 发表于 2022-9-19 16:00
老师好，这样就奇怪了，这个输入文件我在集群上复现了无数次，集群上安装的是gaussian16 C.01版本

开壳层过渡金属用不同版本的量化程序、不同机器甚至同一机器提交多次确实可能出现无法复现的现象，但一般是初猜和计算过程无法复现，稳定的波函数计算结果是可以复现的。不必纠结于是否能复现，结果合理就行。

Stepheee

zjxitcc 发表于 2022-9-19 16:37
开壳层过渡金属用不同版本的量化程序、不同机器甚至同一机器提交多次确实可能出现无法复现的现象，但一般 ...

好的，多谢老师指点。