Molpro算CCSD(T)出现HF不收敛，尝试降低精度和增加迭代次数仍无法解决

Scx

Molpro算CCSD(T)出现HF不收敛，用了几种方法修改过后还是不行，想问下，输出文件显示信息是否说明基组用的不合适或者是还需要定义开壳层或者闭壳层轨道之类的信息。如下：
Convergence thresholds:  THRVAR = 1.00D-08  THRDEN = 1.00D-06

CCSD(T)     terms to be evaluated (factor= 1.000)


Number of core orbitals:          29 (  29 )
Number of closed-shell orbitals:  14 (  14 )
Number of active  orbitals:        2 (   2 )
Number of external orbitals:     278 ( 278 )

Memory could be reduced to 587.50 Mwords without degradation in triples

Number of N-1 electron functions:              30
Number of N-2 electron functions:             435
Number of singly external CSFs:              8400
Number of doubly external CSFs:          25803260
Total number of CSFs:                    25811660

Molecular orbitals read from record     2100.2  Type=RHF/CANONICAL (state 1.1)

?Error: RHF not converged. This error exit can be avoided using the IGNORE_ERROR option on the ORBITAL directive



输入文件：
***,
memory,2000,m
print,orbitals=5

angstrom,cartesian
symmetry,nosym
orient,noorient
geometry={
坐标
}

basis,Cr=aug-cc-pVTZ,O=6-311+G*,N=6-311+G*
{hf;maxit,200
wf,88,1,2}
{rccsd(t)}


DFT下state=3-A，对称性Cs，多重度为3，所以找的不可约序号1；
写总电子数88时，出现HF不收敛，降低收敛精度和增加迭代次数都无法解决；
只写最外层电子总数时，没有HF不收敛，正常结束。
我想问下，wf后面elec必须是全部电子总数吗？什么情况下是不写电子总数的呢？

sobereva

读取不带弥散函数的收敛的结果当初猜再试，或者其它程序自带的帮助收敛的算法。
你这基组实在不怎么搭。6-311+G*给CCSD(T)用可谓烂的一塌糊涂。如果体系不带负电荷，就全都用cc-pVTZ就完了，没弥散函数时候容易收敛得多
肯定得写总电子数，要不然只写价电子数就成了算阳离子了，阳离子状态肯定SCF容易收敛得多，但结果并不是你要的

Mikasa

看手册16.12 Handling difficult cases: When SCF does not converge
这是网页版
https://www.molpro.net/info/release/doc/manual/node180.html

Scx

sobereva 发表于 2019-4-24 21:48
读取不带弥散函数的收敛的结果当初猜再试，或者其它程序自带的帮助收敛的算法。
你这基组实在不怎么搭。6- ...

谢谢老师！

Scx

Mikasa 发表于 2019-4-24 21:49
看手册16.12 Handling difficult cases: When SCF does not converge
这是网页版
https://www.molpro.net ...

谢谢！