| 查看: 931 | 回复: 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 ] |
回复此楼» 猜你喜欢
26年材料/化学相关专业博士研究生招生 温州大学化学与材料工程学院(26年秋季入学)
已经有0人回复
-
求助电催化剂制备，怎么在碳布表面生长PdCo合金
已经有0人回复
-
物理化学论文润色/翻译怎么收费?
已经有100人回复
-
求助电催化剂制备，怎么在碳布表面生长PdCo合金
已经有1人回复
-
304求调剂
已经有5人回复
-
上海华东理工大学机械与动力工程学院先进材料与制造研究组招聘公告
已经有0人回复
-
304求调剂
已经有6人回复
-
304求调剂
已经有1人回复
-
中药学306，求调剂
已经有0人回复
-
河南一本高校化学，材料，化工专业招收研究生
已经有0人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 使用AdNDP方法以及ELF/LOL、多中心键级研究多中心键
已经有14人回复
4楼2013-03-16 20:45:48
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +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
