求助，ccsd(t)算单点，rhf不收敛

xuejunfei1987

以前没用过molpro，最近做ccsd(t)算单点，我用B3lyp/6-311++g(3df,3pd)优化的结构，没几秒钟就出错，总是说?Error: RHF not converged. This error exit can be avoided using the IGNORE_ERROR option on the ORBITAL directive。能级移动，改变迭代也试了，离收敛还差远呢，scf收敛趋势近乎一条直线，希望帮我解决一下，不明白为什么单点都会出错。下面是我的输入文件
***,NH3-HOmemory,4000,m

geometry={
  N        -1.260766000000     -0.000011000000     -0.000037000000
  H        -1.640632000000     -0.848156000000     -0.402540000000
  H        -1.640413000000      0.774001000000     -0.531506000000
  H        -1.639505000000      0.074739000000      0.936334000000
  O         1.646329000000     -0.000051000000     -0.000184000000
  H         0.661544000000     -0.000087000000     -0.000288000000
}

basis=avtz

hf

uccsd(t)

xuejunfei1987

如果可以，希望有人帮我跑一下，几分钟应该就能算完

zjxitcc

（1）应该自己在原贴中编辑，而非自己回复自己多加一层楼。
（2）这个实际上是ROHF不收敛，不是RHF不收敛。Molpro把RHF和ROHF方法统称使用rhf关键词，我认为这在其他量化程序里是比较少见的。
（3）用高斯做一个一样的ROHF计算，然后传轨道给Molpro读入即可解决。高斯输入文件如下

1. %chk=NH3_OH_AVTZ.chk
   
2. %mem=56GB
   
3. %nprocshared=28
   
4. #p ROHF/aug-cc-pVTZ nosymm int=nobasistransform
   
5. 
   
6. title
   
7. 
   
8. 0 2
   
9. N        -1.260766000000     -0.000011000000     -0.000037000000
   
10. H        -1.640632000000     -0.848156000000     -0.402540000000
    
11. H        -1.640413000000      0.774001000000     -0.531506000000
    
12. H        -1.639505000000      0.074739000000      0.936334000000
    
13. O         1.646329000000     -0.000051000000     -0.000184000000
    
14. H         0.661544000000     -0.000087000000     -0.000288000000
    
15. 
    

复制代码

用MOKIT中的fch2com小程序传轨道（Release里也有Win版预编译的可以直接使用）

1. fch2com NH3_OH_AVTZ.chk

复制代码

或

1. fch2com NH3_OH_AVTZ.fch

复制代码

生成molpro输入文件NH3_OH_AVTZ.com和轨道文件nh3_oh_avtz.a，ROHF会在1-2圈收敛，马上进入CCSD(T)计算。你可以自己在molpro输入文件末尾加上你要的CCSD(T)之类的。

dlldwtz

molpro2022版本提供了hf的一个新选项，利用hf,so-sci可以明显帮助收敛，你这个结构40次迭代即可收敛

xuejunfei1987

zjxitcc 发表于 2022-9-27 22:00
（1）应该自己在原贴中编辑，而非自己回复自己多加一层楼。
（2）这个实际上是ROHF不收敛，不是RHF不收敛 ...

（1）好的，下次我注意
（2）非常感谢你的帮助，我学习了很多，按照你的步骤我试了一下真的可以。还想请教一下，之前没接触过molpro，那个生成的com文件中的basis后面一堆信息，是什么意思啊，我之前提交作业中的basis=avtz，和这个有什么关系吗

xuejunfei1987

dlldwtz 发表于 2022-9-28 08:54
molpro2022版本提供了hf的一个新选项，利用hf,so-sci可以明显帮助收敛，你这个结构40次迭代即可收敛

谢谢，可是我这不是2022版本的呀

zjxitcc

xuejunfei1987 发表于 2022-9-28 19:12
（1）好的，下次我注意
（2）非常感谢你的帮助，我学习了很多，按照你的步骤我试了一下真的可以。还想请 ...

fch2com直接生成basis详细数据，用户不需要再写basis = avtz。你平时写basis = avtz，molpro程序内部去它自己基组目录下找基组数据，两种做法是等价的。

xuejunfei1987

zjxitcc 发表于 2022-9-28 19:16
fch2com直接生成basis详细数据，用户不需要再写basis = avtz。你平时写basis = avtz，molpro程序内部去它 ...

嗯 好的 明白啦 太感谢你了  

dlldwtz

xuejunfei1987 发表于 2022-9-28 19:13
谢谢，可是我这不是2022版本的呀

更新一下版本呗