24小时热门版块排行榜    

查看: 3607  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 25/1250 2026-08-16 10:53 by hmhminy
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 时间戳又变了8-15 +12 archvillain 2026-08-15 23/1150 2026-08-15 20:47 by foyo404
[基金申请] filecode=后面第一个是大写字母 +6 wangze12014 2026-08-14 7/350 2026-08-15 20:12 by gltch
[基金申请] filecode +8 cratir 2026-08-14 12/600 2026-08-15 18:08 by zyfgau
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 是这周出结果还是下周出结果? +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +23 Lanmanbaby 2026-08-09 37/1850 2026-08-14 13:40 by Equinoxhua
[硕博家园] 请教兼职经验 +3 是阿文鸭 2026-08-09 3/150 2026-08-14 12:07 by HER12025
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] filecode +15 documentary 2026-08-10 17/850 2026-08-14 10:08 by kissu88
[基金申请] 静等基金结果 +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 by 且听虎啸
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
信息提示
请填处理意见