24小时热门版块排行榜    

Znn3bq.jpeg
汕头大学海洋科学接受调剂
查看: 970  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿华南理工大学331分材料求调剂 +10 天下ww 2026-04-09 11/550 2026-04-13 23:25 by pies112
[考研] 一志愿厦大生物学332求调剂 +11 池池池池池池 2026-04-08 11/550 2026-04-13 14:10 by 科研论
[考研] 考研二轮调剂 +10 故人?? 2026-04-09 10/500 2026-04-13 09:55 by szhize
[考研] 322求调剂,08工科 +4 今天是个小号 2026-04-08 4/200 2026-04-13 00:20 by baobaoye
[考研] 344 材料专业 求调剂211 无地域要求 +8 hualkop 2026-04-11 8/400 2026-04-12 22:24 by fqwang
[考研] 291求调剂 +11 关忆北. 2026-04-09 12/600 2026-04-12 10:32 by 逆水乘风
[考研] 一志愿西北工业大学289 085602 +33 yang婷 2026-04-10 34/1700 2026-04-12 08:11 by Art1977
[考研] 277 数一104,学硕,求调剂 +21 瓶子PZ 2026-04-09 23/1150 2026-04-11 23:12 by labixiaoqiao
[考研] 284求调剂 +11 archer.. 2026-04-09 12/600 2026-04-11 20:23 by 蓝云思雨
[考研] 一志愿211,化学310分,本科重点双非,求调剂 +23 努力奋斗112 2026-04-08 23/1150 2026-04-10 23:29 by 314126402
[考研] 309求调剂 +14 wdhw 2026-04-10 15/750 2026-04-10 21:06 by zhouxiaoyu
[考研] 263能源动力专硕求调剂 +4 加大号饭盒袋 2026-04-10 4/200 2026-04-10 20:52 by gong120082
[考研] 吉大计算机技术331分,英语六级,求调剂 +3 峰峰021116 2026-04-09 3/150 2026-04-10 20:01 by chemisry
[考研] 一志愿京区985,085401电子信息,本科电子信息 +3 阳光开朗的男孩 2026-04-10 3/150 2026-04-10 16:29 by sophia_93
[考研] 071000生物学调剂求助 +17 zzzzwww 2026-04-09 20/1000 2026-04-10 15:55 by 求调剂zz
[考研] 江苏大学 工科调剂 捡漏 +3 Evan_Liu 2026-04-09 5/250 2026-04-10 10:22 by Evan_Liu
[考研] 331求调剂 +5 luoxin0706. 2026-04-08 5/250 2026-04-08 22:15 by zhouyuwinner
[考研] 275 求调剂 +8 Lei812514 2026-04-07 8/400 2026-04-08 12:46 by chemisry
[考研] 电子信息346 +4 zuoshaodian 2026-04-08 4/200 2026-04-08 11:54 by zzucheup
[考研] 307求调剂 +3 Youth@@ 2026-04-07 3/150 2026-04-07 22:00 by hemengdong
信息提示
请填处理意见