计算化学公社

标题: 求助，ccsd(t)算单点，rhf不收敛 [打印本页]

作者
Author: xuejunfei1987 时间: 2022-9-27 21:47
标题: 求助，ccsd(t)算单点，rhf不收敛
以前没用过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)

作者
Author: xuejunfei1987 时间: 2022-9-27 21:48
如果可以，希望有人帮我跑一下，几分钟应该就能算完

作者
Author: zjxitcc 时间: 2022-9-27 22:00
本帖最后由 zjxitcc 于 2022-9-27 22:22 编辑

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

%chk=NH3_OH_AVTZ.chk
%mem=56GB
%nprocshared=28
#p ROHF/aug-cc-pVTZ nosymm int=nobasistransform
title
0 2
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

复制代码

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

fch2com NH3_OH_AVTZ.chk

复制代码

或

fch2com NH3_OH_AVTZ.fch

复制代码

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

作者
Author: dlldwtz 时间: 2022-9-28 08:54
本帖最后由 dlldwtz 于 2022-9-28 08:57 编辑

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

作者
Author: xuejunfei1987 时间: 2022-9-28 19:12

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

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

作者
Author: xuejunfei1987 时间: 2022-9-28 19:13

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

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

作者
Author: zjxitcc 时间: 2022-9-28 19:16

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

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

作者
Author: xuejunfei1987 时间: 2022-9-28 19:30

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

嗯好的明白啦太感谢你了

作者
Author: dlldwtz 时间: 2022-9-29 09:05

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

更新一下版本呗

欢迎光临计算化学公社 (http://bbs.keinsci.com/)