24小时热门版块排行榜    

Znn3bq.jpeg
查看: 964  |  回复: 3

kent1022

金虫 (著名写手)

[求助] 闭壳层体系单点能扫描出现的问题

对一个简单的负一价双原子分子体系(闭壳层)用ccsd(t)/avdz进行单点能扫描,发现在某些核间距时迭代次数明显逐渐增大,而且出现若干不同的Spin Case,如下。此时计算还在进行,把输出文件在gview中打开发现势能曲线变得不光滑,单点能的能量有上下起伏的情况。问:  出现这种情况是什么原因,怎样避免?
………
Iteration Nr. 228
**********************
DD1Dir will call FoFMem   1 times, MxPair=         6
NAB=     3 NAA=     0 NBB=     0.
DE(Corr)= -0.23891926     E(CORR)=     -21.857781037     Delta= 5.22D-08
NORM(A)=   0.34013606D+01
Dominant configurations:
***********************
Spin Case        I    J    A    B          Value
    AA            3         7           -0.106518D+00
    AA            3        15           -0.146622D+00
    AA            4         8           -0.106518D+00
    AA            4        16           -0.146622D+00
   ABAB           3    3    5    5      -0.189012D+01
   ABAB           3    3    5    6       0.339414D+00
   ABAB           3    3    5    9      -0.302951D+00
   ABAB           3    3    5   13      -0.199322D+00
   ABAB           3    3    5   14      -0.318367D+00
   ABAB           3    3    5   17       0.163331D+00
   ABAB           3    3    5   28       0.224230D+00
   ABAB           3    3    6    5       0.339414D+00
   ABAB           3    3    9    5      -0.302951D+00
   ABAB           3    3   13    5      -0.199322D+00
   ABAB           3    3   13   13       0.131320D+00
   ABAB           3    3   13   14      -0.133078D+00
   ABAB           3    3   13   17      -0.155731D+00
   ABAB           3    3   14    5      -0.318367D+00
   ABAB           3    3   14   13      -0.133078D+00
   ABAB           3    3   14   14       0.173566D+00
   ABAB           3    3   14   17       0.119876D+00
   ABAB           3    3   14   18      -0.165972D+00
   ABAB           3    3   14   28      -0.196664D+00
   ABAB           3    3   17    5       0.163331D+00
   ABAB           3    3   17   13      -0.155731D+00
   ABAB           3    3   17   14       0.119876D+00
   ABAB           3    3   17   17       0.176752D+00
   ABAB           3    3   18   14      -0.165972D+00
   ABAB           3    3   18   28       0.139280D+00
   ABAB           3    3   28    5       0.224230D+00
   ABAB           3    3   28   14      -0.196664D+00
   ABAB           3    3   28   18       0.139280D+00
   ABAB           3    3   28   28       0.195483D+00
   ABAB           4    4    5    5      -0.189013D+01
   ABAB           4    4    5    6       0.339415D+00
   ABAB           4    4    5    9      -0.302952D+00
   ABAB           4    4    5   13      -0.199323D+00
   ABAB           4    4    5   14      -0.318367D+00
   ABAB           4    4    5   17       0.163331D+00
   ABAB           4    4    5   28       0.224230D+00
   ABAB           4    4    6    5       0.339415D+00
   ABAB           4    4    9    5      -0.302952D+00
   ABAB           4    4   13    5      -0.199323D+00
   ABAB           4    4   13   13       0.131320D+00
   ABAB           4    4   13   14      -0.133078D+00
   ABAB           4    4   13   17      -0.155731D+00
   ABAB           4    4   14    5      -0.318367D+00
   ABAB           4    4   14   13      -0.133078D+00
   ABAB           4    4   14   14       0.173566D+00
   ABAB           4    4   14   17       0.119876D+00
   ABAB           4    4   14   18      -0.165972D+00
   ABAB           4    4   14   28      -0.196664D+00
   ABAB           4    4   17    5       0.163331D+00
   ABAB           4    4   17   13      -0.155731D+00
   ABAB           4    4   17   14       0.119876D+00
   ABAB           4    4   17   17       0.176752D+00
   ABAB           4    4   18   14      -0.165972D+00
   ABAB           4    4   18   28       0.139280D+00
   ABAB           4    4   28    5       0.224230D+00
   ABAB           4    4   28   14      -0.196664D+00
   ABAB           4    4   28   18       0.139280D+00
   ABAB           4    4   28   28       0.195484D+00
Largest amplitude= 1.89D+00
Time for triples=        8.56 seconds.
T4(CCSD)= -0.33481841D+00
T5(CCSD)=  0.74039746D-02
CCSD(T)= -0.22185195476D+02

另:同时对该体系用ccsd(t)/avtz计算,发现在核间距较小时就出现超过MAX. CYCLES的情况(输入文件中设置maxcyc=1000),并且列出了更多的Spin Case,同时计算终止。这又是什么原因?

…………
   ABAB           4    4   29   14       0.105618D+00
   ABAB           4    4   34   13      -0.157926D+00
   ABAB           4    4   34   14       0.135242D+00
   ABAB           4    4   34   17       0.129708D+00
   ABAB           4    4   37   14      -0.117016D+00
Largest amplitude= 2.12D+00
Error termination via Lnk1e in /home/g09/l913.exe at Sun Mar 10 18:25:27 2013.

[ Last edited by kent1022 on 2013-3-11 at 19:51 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kent1022

金虫 (著名写手)

2楼2013-03-13 12:34:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangf44

禁虫 (小有名气)

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
kent1022: 金币+20 2013-03-16 20:44:12
kent1022: 金币+30 2013-03-18 10:30:09
本帖内容被屏蔽

3楼2013-03-14 16:48:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kent1022

金虫 (著名写手)

引用回帖:
3楼: Originally posted by wangf44 at 2013-03-14 16:48:54
表明CC方法的结果不可靠。有可能是这些构型下CC不可靠,还有可能是HF结果不合理。

不过你可以检查一下HF能量,看是不是也是不光滑的。

查了一下,用HF计算的结果没有这种情况。
这是不是说明当前体系不适合用cc方法计算?
4楼2013-03-16 20:45:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 kent1022 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 277求调剂 +23 倪建设 2026-04-06 23/1150 2026-04-13 00:40 by 王童子龙
[教师之家] 山东双非院校考核超级无底线,领导幸灾乐祸,教师遭殃恐 +3 qut2026 2026-04-11 7/350 2026-04-12 20:24 by qut2026
[考研] 0831生医工第一轮调剂失败求助 +12 小熊睿睿_s 2026-04-11 16/800 2026-04-12 16:28 by 钰璞
[考研] 电气工程专硕320求调剂 +5 小麻子111 2026-04-10 5/250 2026-04-12 10:47 by zhouyuwinner
[找工作] 山东高校教师考核超级无底线,员工过不下去啦 +4 qut2026 2026-04-09 9/450 2026-04-12 00:54 by qut2026
[考研] 331求调剂 +5 王国帅 2026-04-11 5/250 2026-04-11 22:56 by 溪涧流水
[考研] 086003调剂求助 +21 苏弋万 2026-04-09 22/1100 2026-04-11 20:25 by dongdian1
[考研] 调剂求助 +6 果然有我 2026-04-11 7/350 2026-04-11 16:22 by 明月此时有
[考研] 化学308分求调剂 +22 你好明天你好 2026-04-07 24/1200 2026-04-11 11:14 by ChemPharm
[考研] 一志愿985机械学硕380求调剂 +5 关关雎鸠10 2026-04-11 5/250 2026-04-11 10:10 by 知念。A
[考研] 22408 352分求调剂0854类 +4 努力的夏末 2026-04-09 4/200 2026-04-11 09:57 by zhq0425
[考研] 一志愿211 0703化学 346分求调剂 +22 土豆er? 2026-04-09 23/1150 2026-04-10 10:58 by 高维春
[考研] 调剂申请086000一志愿西北农林科技大学生物与医药320分-本科齐鲁工业大学 +3 美美女士 2026-04-09 3/150 2026-04-10 10:31 by liuhuiying09
[考研] 278求调剂 +27 范婷娜 2026-04-07 31/1550 2026-04-09 20:49 by zhouxiaoyu
[考研] 336求调剂,一志愿中科大 +9 墨彧 yuyu 2026-04-06 9/450 2026-04-08 11:24 by 想读书的菌菌
[考研] 338求调剂 +8 wxygxsaaaaa 2026-04-06 8/400 2026-04-08 06:58 by 无际的草原
[考研] 求调剂 +15 熊二想上岸 2026-04-06 15/750 2026-04-08 04:53 by 无际的草原
[考研] 388求调剂 +6 四川王涛 2026-04-07 8/400 2026-04-08 00:17 by JourneyLucky
[考研] 考研调剂 +3 Wwwwwww哇 2026-04-06 3/150 2026-04-06 20:55 by lbsjt
[考研] 362求调剂一志愿中国石油大学 +4 我要考大 2026-04-06 6/300 2026-04-06 14:11 by 无际的草原
信息提示
请填处理意见