| 查看: 973 | 回复: 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 ] |
回复此楼» 猜你喜欢
调剂招生
已经有0人回复
-
调剂招生
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有187人回复
-
调剂招生
已经有0人回复
-
调剂招生
已经有1人回复
-
实际气体内压强的微观机制或将改写--------从一个基本的力学佯谬说起
已经有1人回复
-
调剂招生
已经有0人回复
-
调剂招生
已经有0人回复
-
调剂招生
已经有5人回复
-
招收0703开头的学硕调剂生 线上复试
已经有2人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 使用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
