计算化学公社

标题: stable=opt稳定后用stable关键字检测显示不稳定 [打印本页]
作者
Author: Stepheee    时间: 2022-9-19 11:11
标题: stable=opt稳定后用stable关键字检测显示不稳定
最近在复现文献中的一个数据时遇到了奇怪的错误，计算的是一个五重态的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

  6. FeO_5

  8. 0 5
  9.   Fe
  10.   O 1 B1

  12. B1=1.618


复制代码
提示"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

  6. FeO_5

  8. 0 5
  9. Fe
  10. O 1 B1

  12. B1=1.618

复制代码
继续计算后提示"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

  6. FeO_5

  8. 0 5
  9. Fe
  10. O 1 B1

  12. B1=1.618

复制代码
则会再次提示"internal instability"。



作者
Author: zjxitcc    时间: 2022-9-19 13:12
本帖最后由 zjxitcc 于 2022-9-19 13:16 编辑

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凑热闹。
作者
Author: hebrewsnabla    时间: 2022-9-19 14:09
本帖最后由 hebrewsnabla 于 2022-9-19 14:13 编辑

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

又用scf(fermi)算了一遍，还是stable
作者
Author: Stepheee    时间: 2022-9-19 16:00
hebrewsnabla 发表于 2022-9-19 14:09
我本来想复现一下这个例子，但是我第一步用了scf(vshift=200,xqc)以后直接就stable了……

又用scf(fermi ...

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

开壳层过渡金属用不同版本的量化程序、不同机器甚至同一机器提交多次确实可能出现无法复现的现象，但一般是初猜和计算过程无法复现，稳定的波函数计算结果是可以复现的。不必纠结于是否能复现，结果合理就行。
作者
Author: Stepheee    时间: 2022-9-19 17:28
zjxitcc 发表于 2022-9-19 16:37
开壳层过渡金属用不同版本的量化程序、不同机器甚至同一机器提交多次确实可能出现无法复现的现象，但一般 ...

好的，多谢老师指点。



欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3