24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 3355  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂 +8 xzghyuj 2026-04-04 8/400 2026-04-06 06:47 by houyaoxu
[考研] 308求调剂 +3 终不似从前 2026-04-05 3/150 2026-04-05 22:23 by hemengdong
[考研] 求助071001调剂!!! +3 黄守松 2026-04-05 4/200 2026-04-05 21:01 by barlinike
[考研] 材料调剂 +13 一样YWY 2026-04-03 14/700 2026-04-05 18:20 by 蓝云思雨
[考研] 材料0856 英一数二 323 求调剂 +14 袁sy 2026-04-01 14/700 2026-04-05 18:18 by cql1109
[考研] 0703化学 +12 妮妮ninicgb 2026-04-04 13/650 2026-04-05 10:46 by 啊俊!
[考研] 一志愿华中农业大学0710(A)初试329分 求调剂 +4 一名26考研生 2026-04-04 4/200 2026-04-05 10:01 by barlinike
[考研] 材料调剂 +12 一样YWY 2026-04-04 12/600 2026-04-05 08:24 by 544594351
[考研] 309求调剂 +4 快乐的小白鸽 2026-04-04 5/250 2026-04-04 15:55 by cql1109
[考研] 材料295 +13 小英11 2026-04-03 14/700 2026-04-04 09:02 by 来看流星雨10
[考研] 考研调剂 +5 小sun要好运 2026-04-03 5/250 2026-04-03 21:43 by 啵啵啵0119
[考研] 11408,284分,二战真诚求调剂 +4 12.27 2026-04-02 4/200 2026-04-03 14:14 by dxiaoxin
[考研] 0705理学294求调剂 +3 成果成果cg5 2026-04-03 3/150 2026-04-03 14:04 by simons1972
[考研] 285求调剂 +6 FZAC123 2026-03-30 6/300 2026-04-03 12:22 by xingguangj
[考研] 0703化学 +7 goldtt 2026-04-02 9/450 2026-04-03 09:36 by 蓝云思雨
[考研] 一志愿陕西师范大学生物学317分 +5 1563日。 2026-04-02 5/250 2026-04-03 06:58 by ilovexiaobin
[考研] 282求调剂 +13 呼吸都是减肥 2026-04-01 13/650 2026-04-02 14:10 by baoball
[考研] 化学工程专硕324分,一志愿中国矿业大学求调剂 +7 耿耿1314 2026-04-01 7/350 2026-04-02 07:40 by 尚水阁主
[考研] 材料调剂 +12 一样YWY 2026-04-01 12/600 2026-04-02 00:21 by 百秒光年
[考研] 085602化学工程268分蹲调剂 +8 月照花林。 2026-04-01 8/400 2026-04-01 22:08 by 无际的草原
信息提示
请填处理意见