24小时热门版块排行榜    

查看: 3613  |  回复: 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的回帖

moulh

禁虫 (小有名气)

本帖内容被屏蔽

7楼2018-06-24 22:59:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +32 医学老男孩 2026-08-13 72/3600 2026-08-16 13:15 by 医学老男孩
[基金申请] 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
[基金申请] 关于Filecode分析方法 +10 majunge000 2026-08-10 13/650 2026-08-15 12:26 by hanpeng972
[基金申请] 是这周出结果还是下周出结果? +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +8 布布和一二 2026-08-10 11/550 2026-08-14 14:58 by Equinoxhua
[硕博家园] 请教兼职经验 +3 是阿文鸭 2026-08-09 3/150 2026-08-14 12:07 by HER12025
[基金申请] 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 +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +11 yufeiwaner 2026-08-09 13/650 2026-08-12 08:17 by yufeiwaner
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 什么时候出结果,有咨询渠道??? +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
信息提示
请填处理意见