| 查看: 930 | 回复: 3 | |||
[求助]
闭壳层体系单点能扫描出现的问题
|
|
对一个简单的负一价双原子分子体系（闭壳层）用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 ] |
回复此楼» 猜你喜欢
关于引力延迟下双星系统角动量不守恒的疑惑---求高手解答！
已经有10人回复
-
CSC访学博后项目获批，外方学校暂停合作该怎么办？
已经有51人回复
-
物理化学论文润色/翻译怎么收费?
已经有161人回复
-
求标准卡片Cu0.4In0.4Zn1.2S2 PDF#32-0340
已经有2人回复
-
量子点电致发光器件（QLED）：蓝光量子点（QDs）和纳米氧化锌电子传输材料
已经有0人回复
-
需要一个圆偏振光诱导化学反应，请问使用圆偏振光作为光源，应该需要哪些设备搭建
已经有0人回复
-
分子动力学模拟合作需求
已经有2人回复
-
求ENDNOTE破解软件安装包
已经有4人回复
-
河海大学 （211.双一流高校）钙钛矿光电实验室丁勇团队欢迎有志青年！
已经有3人回复
-
黑龙江省原子能研究院化学课题组招收2026级化学专业或化学相关专业硕士生
已经有0人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 使用AdNDP方法以及ELF/LOL、多中心键级研究多中心键
已经有14人回复
 |
2楼2013-03-13 12:34:30
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
kent1022: 金币+20 2013-03-16 20:44:12
kent1022: 金币+30 2013-03-18 10:30:09
感谢参与，应助指数 +1
kent1022: 金币+20 2013-03-16 20:44:12
kent1022: 金币+30 2013-03-18 10:30:09
|
本帖内容被屏蔽 |
3楼2013-03-14 16:48:54
4楼2013-03-16 20:45:48
10