24小时热门版块排行榜    

查看: 3220  |  回复: 6

kent1022

金虫 (著名写手)

[求助] uccsd(t)-f12计算不收敛

对一开壳层体系用rhf  +  ccsd(t)-f12扫描单点能,出现收敛问题,
Starting UCCSD calculation

ITER.      SQ.NORM     CORR.ENERGY   TOTAL ENERGY   ENERGY CHANGE        DEN1      VAR(S)    VAR(P)  DIIS     TIME
   1      1.16866374    -0.31973648  -173.84507884    -0.31973648    -0.04368108  0.52D-01  0.40D-02  1  1     1.95
   2      1.21529237    -0.34466229  -173.87000465    -0.02492582    -0.01618288  0.12D-01  0.59D-02  2  2     2.01
  .....
   8      2.36745749    -0.42799740  -173.95333975    -0.00988771    -0.00090030  0.75D-03  0.35D-03  6  3     2.31
   9      2.45888304    -0.42663949  -173.95198185     0.00135791    -0.00037733  0.16D-03  0.19D-03  6  2     2.37

UNREASONABLE NORM. CALCULATION STOPPED

?ERROR: NO CONVERGENCE. TRIPLES WILL NOT BE DONE. THIS CAN BE AVOIDED USING THE NOCHECK OPTION.

最后显示
? Error
? No convergence. This error exit can be avoided using the NOCHECK option
? The problem occurs in ccsd

增加迭代次数,加上maxit=200,
提示 Input error: Attempt to set Variable MAXIT in command block

加上 NOCHECK 后
{uccsd(t)-f12,df_basis=df,df_basis_exch=jk,ri_basis=ri;wf,19,4,2;NOCHECK }
又提示
Unknown directive NOCHECK for command UCCSD(T)-F12
GLOBAL ERROR checkinp on processor   0                                       

请问怎样解决这个收敛问题?

另:开壳层体系的自旋污染怎样判断是大还是可以忽略?
Spin contamination <S**2-Sz**2-Sz> 的值正常情况下大概是多少?
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luqing6879

木虫 (著名写手)

★ ★
zhou2009: 金币+2 2014-04-03 05:49:53
加上nocheck的结果一般是不可靠的,不建议使用。CCSD(T)不收敛,这种情况一般都是HF没有收敛到正确的态上。
cishan.goodweb.cn/慈善点击,点滴帮助。
2楼2014-04-03 03:51:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kent1022

金虫 (著名写手)

引用回帖:
2楼: Originally posted by luqing6879 at 2014-04-03 03:51:40
加上nocheck的结果一般是不可靠的,不建议使用。CCSD(T)不收敛,这种情况一般都是HF没有收敛到正确的态上。

请问该怎样处理让计算能继续下去,对hf的计算是要作怎样的限制?
3楼2014-04-04 16:26:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yongleli

木虫 (正式写手)

f12据我的经验,只在稳定点附近好使,扫描势能面还是推荐不用f12,加大基组。
4楼2014-04-05 02:07:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luqing6879

木虫 (著名写手)

一个是可以试试显示地写上occ和closed卡片,使hf收敛到符合化学直觉上的电子态。
另一个是在hf之前跑一个casscf,以它来做初猜,喂给hf。
如果非常确定hf收敛到了正确的态,而ccsd不收敛,可以试试给ccsd加level shift。不过如之前所说,如果ccsd不收敛,或几十步内没有收敛,一般还是应该在hf阶段找问题。
cishan.goodweb.cn/慈善点击,点滴帮助。
5楼2014-04-05 05:02:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luqing6879

木虫 (著名写手)

不知道你有没有试过计算对应的闭壳层分子。看看它是否能成功收敛。能收敛的话可以以它为初猜,喂给rohf
cishan.goodweb.cn/慈善点击,点滴帮助。
6楼2014-04-05 05:09:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

moulh

禁虫 (小有名气)

本帖内容被屏蔽

7楼2018-06-24 22:59:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 kent1022 的主题更新
信息提示
请填处理意见