24小时热门版块排行榜    

查看: 1365  |  回复: 18
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

xiemin8761

铁杆木虫 (著名写手)

[求助] 2D-scan求助

根据一篇文章,我尝试用2D-scan扫描一个分子,计算的都失败了不知什么原因,请各位高手指教。
我的输入文件是:
# opt(maxcyc=100,addredundant) upbe1pbe/6-31+g* scf=(qc,nosymm) geom=connectivity

Title Card Required

0 1
C              
H                  1            B1
H                  1            B2    2            A1
H                  1            B3    2            A2    3            D1    0
C                  1            B4    2            A3    3            D2    0
C                  5            B5    1            A4    2            D3    0
H                  6            B6    5            A5    1            D4    0
H                  6            B7    5            A6    1            D5    0
H                  6            B8    5            A7    1            D6    0
O                  5            B9    1            A8    6            D7    0
O                 10           B10    5            A9    1            D8    0
H                 11           B11   10           A10    5            D9    0
H                 11           B12   10           A11    5           D10    0

   B1             1.09041428
   B2             1.09834559
   B3             1.09588262
   B4             1.50403959
   B5             1.48882226
   B6             1.08944121
   B7             1.10029571
   B8             1.09377481
   B9             1.22943504
   B10            2.25904576
   B11            0.97274001
   B12            0.97271156
   A1           109.92873208
   A2           111.12632155
   A3           112.20520707
   A4           120.67941448
   A5           111.38556342
   A6           106.61508768
   A7           109.59825862
   A8           115.77988206
   A9           112.26541245
   A10           97.56795174
   A11           99.56438207
   D1          -119.61623552
   D2           118.83803220
   D3           169.57863237
   D4           161.78293646
   D5           -79.13473791
   D6            36.72209100
   D7          -179.31990089
   D8          -177.29025062
   D9           129.69444934
   D10         -123.23665237

1 2 1.0 3 1.0 4 1.0 5 1.0
2
3
4
5 6 1.0 10 2.0
6 7 1.0 8 1.0 9 1.0
7
8
9
10
11 12 1.0 13 1.0
12
13

7 11 2.33673 s 10 -0.153
7 10 2.5872 s 10 -0.178

输出文件显示的错误为:

Iteration100 RMS(Cart)=  0.00002839 RMS(Int)=  0.03650170
New curvilinear step not converged.
RedCar failed.
Error termination via Lnk1e in /usr/local/intel/g03/l103.exe at Thu Nov  8 21:54:38 2012.
Job cpu time:  0 days  3 hours  1 minutes 54.0 seconds.
File lengths (MBytes):  RWF=     29 Int=      0 D2E=      0 Chk=     19 Scr=      1
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiemin8761

铁杆木虫 (著名写手)

★ ★
cenwanglai: 金币+2, 谢谢分享~ 2012-11-18 16:55:22
我昨天收到高斯客服的邮件,
Thank you for your inquiry. The error here indicates that the program was not able to make the modification of the structure to compute the next scan point without changing other internal coordinates. The redundant internal coordinates are very convenient but the drawback of having the redundancy in the internal coordinates is that this type of situations can happen, i.e. a change in one internal coordinate cannot be done without altering other internal coordinates and/or there are multiple solutions and the program cannot decide for the most appropriate one.

What this means is that, unfortunately, in some cases, one needs to do scans in smaller sections or even computing some points in the scan separately. It sometimes helps using smaller step sizes, you may try for example scan steps of ~ 0.1 Angstrom (yours are > 0.15 Angstrom, which can be too long for certain geometries). Of course, that would also mean that you would need to compute more scan points than what you are doing here. By having somewhat smaller changes, there may be a smaller chance of failure when the program needs to generate the geometry of the next scan point.
Sincerely,

但是如果按原先的方法的话,多加步数好像也不行
15楼2012-11-16 08:23:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 19 个回答

xiemin8761

铁杆木虫 (著名写手)

还有时候出错显示
。。。
Iteration 96 RMS(Cart)=  0.00018544 RMS(Int)=  0.47065122
Iteration 97 RMS(Cart)=  0.00016406 RMS(Int)=  0.47055968
Iteration 98 RMS(Cart)=  0.00014510 RMS(Int)=  0.47047863
Iteration 99 RMS(Cart)=  0.00012829 RMS(Int)=  0.47040689
Iteration100 RMS(Cart)=  0.00011339 RMS(Int)=  0.47034342
New curvilinear step not converged.
RedCar failed in NewRed.
Error termination via Lnk1e in /usr/local/intel/g03/l101.exe at Mon Nov 12 12:59:51 2012.
Job cpu time:  0 days  0 hours  0 minutes  3.3 seconds.
File lengths (MBytes):  RWF=      8 Int=      0 D2E=      0 Chk=      1 Scr=      1
2楼2012-11-12 12:08:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lishijunzong

铁杆木虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
你先看一下这个是否为优化好的体系
3楼2012-11-12 19:43:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiemin8761

铁杆木虫 (著名写手)

引用回帖:
3楼: Originally posted by lishijunzong at 2012-11-12 19:43:50
你先看一下这个是否为优化好的体系

已经是优化好了的,而且我用高斯09计算他竟然可以计算完成,但是不scan奇怪了。
4楼2012-11-13 08:13:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见