24小时热门版块排行榜    

查看: 3310  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 304求调剂 +5 素年祭语 2026-03-15 5/250 2026-03-16 17:00 by 我的船我的海
[考研] 312求调剂 +3 陌宸希 2026-03-16 4/200 2026-03-16 15:06 by peike
[基金申请] 国自科面上基金字体 +5 iwuli 2026-03-12 6/300 2026-03-16 13:13 by Kamiu_MK
[考研] 326求调剂 +4 上岸的小葡 2026-03-15 5/250 2026-03-16 08:39 by Linda Hu
[考研] 求老师收留调剂 +4 jiang姜66 2026-03-14 5/250 2026-03-15 20:11 by Winj1e
[考研] 22408总分284求调剂 +3 InAspic 2026-03-13 3/150 2026-03-15 11:10 by zhq0425
[考研] 【0703化学调剂】-一志愿华中师范大学-六级475 +5 Becho359 2026-03-11 5/250 2026-03-14 11:35 by 哦哦123
[考研] 266求调剂 +4 学员97LZgn 2026-03-13 4/200 2026-03-14 08:37 by zhukairuo
[考研] 一志愿湖师大化学289求调剂 +6 XMCMM3.14159 2026-03-10 6/300 2026-03-14 00:28 by JourneyLucky
[考研] 26考研调剂 +3 ying123. 2026-03-10 3/150 2026-03-14 00:18 by JourneyLucky
[考研] 材料工程,326分,求调剂 +6 KRSLSR 2026-03-10 6/300 2026-03-13 23:47 by JourneyLucky
[考研] 复试调剂 +9 Copy267 2026-03-10 9/450 2026-03-13 23:45 by userper
[考研] 332求调剂 +3 zjy101327 2026-03-11 6/300 2026-03-13 22:48 by JourneyLucky
[考研] 材料与化工求调剂一志愿 985 总分 295 +8 dream…… 2026-03-12 8/400 2026-03-13 22:17 by 星空星月
[考研] 308求调剂 +5 是Lupa啊 2026-03-11 5/250 2026-03-13 22:13 by JourneyLucky
[考研] 材料工程调剂 +9 咪咪空空 2026-03-12 9/450 2026-03-13 22:05 by 星空星月
[考研] 材料与化工085600调剂求老师收留 +9 jiaanl 2026-03-11 9/450 2026-03-13 20:22 by JourneyLucky
[考博] 26读博 +4 Rui135246 2026-03-12 10/500 2026-03-13 07:15 by gaobiao
[考研] 081200-11408-276学硕求调剂 +3 崔wj 2026-03-12 4/200 2026-03-12 19:33 by 求调剂zz
[硕博家园] 木虫好像不热闹了,是不是? +4 偏振片 2026-03-10 4/200 2026-03-10 09:51 by longwave
信息提示
请填处理意见