计算化学公社

标题: 前辈们，请问在用内坐标输入格式部分优化构型计算时为什么老是报错…… [打印本页]

作者
Author: FANYI 时间: 2018-10-8 21:15
标题: 前辈们，请问在用内坐标输入格式部分优化构型计算时为什么老是报错……
大神们：请问一下，我在体系中添加了两个虚原子，想固定前面几个原子的相对位置，进行部分构型优化，这样设置了输入文件，但老是报错（错误提示见最后）……烦请指点迷津！多谢各位！！！
输入文件：%chk=r2-453.chk
%nprocshared=14
%mem=4GB
#p opt b3lyp/genecp

r2-453

0 2
X
X                1 1.00000000
C                2 1.42630500       1 90.0
H                2 1.42630500       1 90.0 3    180.0
O                2 4.53000000       1 a1    3    d1 0
O                5 B1                1 a2    2    d2 0
Cu                1 B2                2 a3    3    d3 0
O                7 B3                1 a4    2    d4 0
O                7 B4                1 a5    2    d5 0
S                9 B5                7 a6    1    d6 0
S                8 B6                7 a7    1    d7 0
O                11 B7                8 a8    7    d8 0
O                11 B8                8 a9    7    d9 0
O                10 B9                9 a10 7    d10 0
O                10 B10                9 a11 7    d11 0
C                11 B11                8 a12 7    d12 0
C                10 B12                9 a13 7    d13 0
F                16 B13             11 a14 8    d14 0
F                16 B14             11 a15 8    d15 0
F                16 B15             11 a16 8    d16 0
F                17 B16             10 a17 9    d17 0
F                17 B17             10 a18 9    d18 0
F                17 B18             10 a19 9    d19 0
C                7 B19                1 a20 2    d20 0
N                24 B20                7 a21 1    d21 0
H                25 B21             24 a22 7    d22 0
N                24 B22                7 a23 1    d23 0
N                24 B23                7 a24 1    d24 0
C                27 B24             24 a25 7    d25 0
H                29 B25             27 a26 24 d26 0
H                29 B26             27 a27 24 d27 0
H                29 B27             27 a28 24 d28 0
C                27 B28             24 a29    7 d29 0
H                33 B29             27 a30 24 d30 0
H                33 B30             27 a31 24 d31 0
H                33 B31             27 a32 24 d32 0
C                28 B32             24 a33    7 d33 0
H                37 B33             28 a34 24 d34 0
H                37 B34             28 a35 24 d35 0
H                37 B35             28 a36 24 d36 0
C                28 B36             24 a37    7 d37 0
H                41 B37             28 a38 24 d38 0
H                41 B38             28 a39 24 d39 0
H                41 B39             28 a40 24 d40 0
C                3 B40                2 a41    1 d41 0
C                45 B41                3 a42 2    d42 0
C                46 B42             45 a43    3 d43 0
C                47 B43             46 a44 45 d44 0
C                48 B44             47 a45 46 d45 0
C                45 B45                3 a46    2 d46 0
H                45 B46                3 a47    2 d47 0
H                46 B47             45 a48    3 d48 0
H                48 B48             47 a49 46 d49 0
H                49 B49             48 a50 47 d50 0
H                50 B50             45 a51    3 d51 0
C                3 B51                2 a52 1    d52 0
H                56 B52                3 a53    2 d53 0
H                56 B53                3 a54    2 d54 0
H                56 B54                3 a55    2 d55 0
C                2 B55                1 a56 3    d56 0
H                60 B56                2 a57    1 d57 0
C                2 B57                1 a58 60 d58 0
H                62 B58                2 a59    1 d59 0
C                62 B59                2 a60    1 d60 0
O                64 B60             62 a61    2 d61 0
O                64 B61             62 a62    2 d62 0
C                66 B62             64 a63 62 d63 0
H                67 B63             66 a64 64 d64 0
H                67 B64             66 a65 64 d65 0
H                67 B65             66 a66 64 d66 0
variables:
B1
B2
B3
B4
B5
B6
B7
B8
B9
B10
B11
B12
B13
B14
B15
B16
B17
B18
B19
B20
B21
B22
B23
B24
B25
B26
B27
B28
B29
B30
B31
B32
B33
B34
B35
B36
B37
B38
B39
B40
B41
B42
B43
B44
B45
B46
B47
B48
B49
B50
B51
B52
B53
B54
B55
B56
B57
B58
B59
B60
B61
B62
B63
B64
B65
a1
a2
a3
a4
a5
a6
a7
a8
a9
a10
a11
a12
a13
a14
a15
a16
a17
a18
a19
a20
a21
a22
a23
a24
a25
a26
a27
a28
a29
a30
a31
a32
a33
a34
a35
a36
a37
a38
a39
a40
a41
a42
a43
a44
a45
a46
a47
a48
a49
a50
a51
a52
a53
a54
a55
a56
a57
a58
a59
a60
a61
a62
a63
a64
a65
a66
d1
d2
d3
d4
d5
d6
d7
d8
d9
d10
d11
d12
d13
d14
d15
d16
d17
d18
d19
d20
d21
d22
d23
d24
d25
d26
d27
d28
d29
d30
d31
d32
d33
d34
d35
d36
d37
d38
d39
d40
d41
d42
d43
d44
d45
d46
d47
d48
d49
d50
d51
d52
d53
d54
d55
d56
d57
d58
d59
d60
d61
d62
d63
d64
d65
d66

Cu 0
sdd
****
C H O N S F 0
6-31g(d)
****

Cu 0
sdd

报错：
   Variables:
  WANTED A FLOATING POINT NUMBER AS INPUT.
  FOUND AN END-OF-LINE FOR INPUT.
B1

   ?
Error termination via Lnk1e in /home/test/gaussian/g09/l101.exe at Tue Oct  9 04:56:14 2018.
Job cpu time:  0 days  0 hours  0 minutes  0.9 seconds.
File lengths (MBytes):  RWF=    5 Int=    0 D2E=    0 Chk=    1 Scr=    1

作者
Author: niobium 时间: 2018-10-8 21:45
老老实实用XYZ坐标不好么。。内坐标下面那些变量是要有数据的，从这个来看你一点都不懂内坐标，而且最后报错信息写得很清楚了，后面需要一个浮点数，你却啥都没有

作者
Author: FANYI 时间: 2018-10-9 09:06

niobium 发表于 2018-10-8 21:45
老老实实用XYZ坐标不好么。。内坐标下面那些变量是要有数据的，从这个来看你一点都不懂内坐标，而且最后报 ...

非常谢谢您的回复~！请问用直角坐标的话，虚原子也可以通过添加-1 来固定其位置吗？然后，这个内坐标的方法是别人告诉我的，或许是我理解错误吧

作者
Author: sobereva 时间: 2018-10-9 18:26

FANYI 发表于 2018-10-9 09:06
非常谢谢您的回复~！请问用直角坐标的话，虚原子也可以通过添加-1 来固定其位置吗？然后，这个内坐标的方 ...

固定虚原子毫无意义，本来也不影响结果
仔细看
在Gaussian中做限制性优化的方法
http://sobereva.com/404

作者
Author: FANYI 时间: 2018-10-11 11:11

sobereva 发表于 2018-10-9 18:26
固定虚原子毫无意义，本来也不影响结果
仔细看
在Gaussian中做限制性优化的方法

非常感谢社长的回复与指点！！

欢迎光临计算化学公社 (http://bbs.keinsci.com/)