24小时热门版块排行榜    

查看: 721  |  回复: 6

weilikang

新虫 (初入文坛)

[求助] 势能面扫描

我之前用这种格式扫描二面角来算旋转势垒,算对了一个,结果其他的都是103错误,客气帮我看看是什么原因吗?怎么改
%chk=Scan2.chk
%mem=200mb
# opt=modredundant b3lyp/6-31g(d) nosymm geom=connectivity

Title Card Required

-1 4
B
C                  1              B1
B                  2              B2    1              A1
B                  3              B3    2              A2    1              D1
C                  1              B4    2              A3    3              D2
B                  1              B5    5              A4    4              D3
C                  6              B6    1              A5    5              D4
B                  7              B7    6              A6    1              D5
B                  8              B8    7              A7    6              D6
C                  6              B9    1              A8    5              D7
Cr                 4             B10    3              A9    2              D8
B                  4             B11    3             A10    2              D9
B                  9             B12    8             A11    7             D10
H                  7             B13    6             A12    1             D11
H                  8             B14    7             A13    6             D12
H                  9             B15    8             A14    7             D13
H                 10             B16    6             A15    1             D14
H                  6             B17    1             A16    5             D15
H                 13             B18    9             A17    8             D16
H                  2             B19    1             A18    5             D17
H                  1             B20    5             A19    4             D18
H                  5             B21    1             A20    2             D19
H                 12             B22    4             A21    3             D20
H                  4             B23    3             A22    2             D21
H                  3             B24    2             A23    1             D22

   B1             1.52185804
   B2             1.52186371
   B3             1.51613597
   B4             1.51862702
   B5             3.04803372
   B6             1.52185804
   B7             1.52186371
   B8             1.51613597
   B9             1.51862702
   B10            1.96437355
   B11            1.57755480
   B12            1.57755480
   B13            1.06912176
   B14            1.17993581
   B15            1.17290604
   B16            1.06636655
   B17            1.17983204
   B18            1.18000000
   B19            1.06912176
   B20            1.17983204
   B21            1.06636655
   B22            1.18000000
   B23            1.17290604
   B24            1.17993581
   A1           107.94017645
   A2           107.98520653
   A3           107.99711290
   A4            90.00000000
   A5            90.00000000
   A6           107.94017645
   A7           107.98520653
   A8            90.00000000
   A9            67.39249694
   A10           61.40105857
   A11           61.40105857
   A12          124.85748594
   A13          125.27852871
   A14          126.55079742
   A15          124.73005829
   A16           90.00000000
   A17          127.42921037
   A18          124.85748594
   A19          125.98865070
   A20          124.73005829
   A21          127.42921037
   A22          126.55079742
   A23          125.27852871
   D1             0.00000000
   D2             0.00000000
   D3            90.00000000
   D4          -107.99711290
   D5            90.00000000
   D6             0.00000000
   D7             0.00000000
   D8            57.18294038
   D9           -44.80048044
   D10           44.80048044
   D11          -90.00000000
   D12          180.00000000
   D13          180.00000000
   D14           90.00000000
   D15          125.98865070
   D16          112.94380877
   D17          180.00000000
   D18          180.00000000
   D19          180.00000000
   D20         -112.94380877
   D21          180.00000000
   D22          180.00000000

1 2 1.0 5 1.0 11 1.0 12 1.0 21 1.0
2 3 1.0 11 1.0 12 1.0 20 1.0
3 4 1.0 11 1.0 12 1.0 25 1.0
4 5 1.0 11 1.0 12 1.0 24 1.0
5 11 1.0 12 1.0 22 1.0
6 7 1.0 10 1.0 11 1.0 13 1.0 18 1.0
7 8 1.0 11 1.0 13 1.0 14 1.0
8 9 1.0 11 1.0 13 1.0 15 1.0
9 10 1.0 11 1.0 13 1.0 16 1.0
10 11 1.0 13 1.0 17 1.0
11
12 23 1.0
13 19 1.0
14
15
16
17
18
19
20
21
22
23
24
25

D 1 12 13 6 0.000000 S 180 10.000000
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gauss98

禁虫 (正式写手)

★ ★
感谢参与,应助指数 +1
cenwanglai: 金币+2, 谢谢回复~ 2012-11-28 10:12:12
本帖内容被屏蔽

2楼2012-11-27 15:13:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

weilikang

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by gauss98 at 2012-11-27 15:13:58
最后一行  的 180 , 你是要扫到180度 还是要扫180步?

扫到 180度的话,要用 180.0 , 要扫 180 步,每部 10。0度 的话,您打算扫 5圈 360 度?


谨慎估计 您是要扫 18 步,笔误成 180

嘿嘿 就是啊 谢谢
3楼2012-11-28 09:25:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

weilikang

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by gauss98 at 2012-11-27 15:13:58
最后一行  的 180 , 你是要扫到180度 还是要扫180步?

扫到 180度的话,要用 180.0 , 要扫 180 步,每部 10。0度 的话,您打算扫 5圈 360 度?


谨慎估计 您是要扫 18 步,笔误成 180

不好意思啊 我改了下还是103错误么 我发给你你帮我看一下啊 谢谢
输入文件%chk=Scan3.chk
%mem=200mb
# opt=modredundant b3lyp/6-31g(d) nosymm geom=connectivity

Title Card Required

0 3
C
B                  1              B1
B                  2              B2    1              A1
B                  3              B3    2              A2    1              D1
C                  1              B4    2              A3    3              D2
C                  1              B5    5              A4    4              D3
B                  6              B6    1              A5    5              D4
B                  7              B7    6              A6    1              D5
B                  8              B8    7              A7    6              D6
C                  6              B9    1              A8    5              D7
Fe                 4             B10    3              A9    2              D8
B                  4             B11    3             A10    2              D9
B                  9             B12    8             A11    7             D10
H                  7             B13    6             A12    1             D11
H                  8             B14    7             A13    6             D12
H                  9             B15    8             A14    7             D13
H                 10             B16    6             A15    1             D14
H                  6             B17    1             A16    5             D15
H                 13             B18    9             A17    8             D16
H                  2             B19    1             A18    5             D17
H                  1             B20    5             A19    4             D18
H                  5             B21    1             A20    2             D19
H                 12             B22    4             A21    3             D20
H                  4             B23    3             A22    2             D21
H                  3             B24    2             A23    1             D22

   B1             1.52185804
   B2             1.52186371
   B3             1.51613597
   B4             1.51862702
   B5             3.04803372
   B6             1.52185804
   B7             1.52186371
   B8             1.51613597
   B9             1.51862702
   B10            1.96437355
   B11            1.57755480
   B12            1.57755480
   B13            1.06912176
   B14            1.17993581
   B15            1.17290604
   B16            1.06636655
   B17            1.17983204
   B18            1.18000000
   B19            1.06912176
   B20            1.17983204
   B21            1.06636655
   B22            1.18000000
   B23            1.17290604
   B24            1.17993581
   A1           107.94017645
   A2           107.98520653
   A3           107.99711290
   A4            90.00000000
   A5            90.00000000
   A6           107.94017645
   A7           107.98520653
   A8            90.00000000
   A9            67.39249694
   A10           61.40105857
   A11           61.40105857
   A12          124.85748594
   A13          125.27852871
   A14          126.55079742
   A15          124.73005829
   A16           90.00000000
   A17          127.42921037
   A18          124.85748594
   A19          125.98865070
   A20          124.73005829
   A21          127.42921037
   A22          126.55079742
   A23          125.27852871
   D1             0.00000000
   D2             0.00000000
   D3            90.00000000
   D4          -107.99711290
   D5            90.00000000
   D6             0.00000000
   D7             0.00000000
   D8            57.18294038
   D9           -44.80048044
   D10           44.80048044
   D11          -90.00000000
   D12          180.00000000
   D13          180.00000000
   D14           90.00000000
   D15          125.98865070
   D16          112.94380877
   D17          180.00000000
   D18          180.00000000
   D19          180.00000000
   D20         -112.94380877
   D21          180.00000000
   D22          180.00000000

1 2 1.0 5 1.0 11 1.0 12 1.0 21 1.0
2 3 1.0 11 1.0 12 1.0 20 1.0
3 4 1.0 11 1.0 12 1.0 25 1.0
4 5 1.0 11 1.0 12 1.0 24 1.0
5 11 1.0 12 1.0 22 1.0
6 7 1.0 10 1.0 11 1.0 13 1.0 18 1.0
7 8 1.0 11 1.0 13 1.0 14 1.0
8 9 1.0 11 1.0 13 1.0 15 1.0
9 10 1.0 11 1.0 13 1.0 16 1.0
10 11 1.0 13 1.0 17 1.0
11
12 23 1.0
13 19 1.0
14
15
16
17
18
19
20
21
22
23
24
25

D 1 12 13 6 0.000000 S 10 18.000000
输出文件New curvilinear step not converged.
Error imposing constraints
Error termination via Lnk1e in /scratch/g03/l103.exe at Thu Nov 29 13:14:21 2012.
Job cpu time:  0 days 12 hours 16 minutes  3.3 seconds.
File lengths (MBytes):  RWF=     85 Int=      0 D2E=      0 Chk=     19 Scr=      1
4楼2012-11-29 10:16:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gauss98

禁虫 (正式写手)

本帖内容被屏蔽

5楼2012-11-29 13:07:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

weilikang

新虫 (初入文坛)

引用回帖:
5楼: Originally posted by gauss98 at 2012-11-29 13:07:53
你都算了12个小时了应该差不多了吧?

差多了,根本就没有扫描出来任何结果啊。你以前是怎么算的啊,有没有出现过这种问题啊
6楼2012-12-10 11:19:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luosizuo

禁虫 (小有名气)

本帖内容被屏蔽

7楼2012-12-10 14:03:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 weilikang 的主题更新
信息提示
请填处理意见