请选择 进入手机版 | 继续访问电脑版

计算化学公社

 找回密码
 现在注册!
查看: 173|回复: 4

[Gaussian/gview] 前辈们,请问在用内坐标输入格式部分优化构型计算时为什么老是报错……

[复制链接]

3

帖子

0

威望

37

eV
积分
40

Level 2 能力者

发表于 2018-10-8 21:15:20 | 显示全部楼层 |阅读模式
大神们:请问一下,我在体系中添加了两个虚原子,想固定前面几个原子的相对位置,进行部分构型优化,这样设置了输入文件,但老是报错(错误提示见最后)……烦请指点迷津!多谢各位!!!
输入文件:%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









127

帖子

1

威望

1340

eV
积分
1487

Level 4 (黑子)

发表于 2018-10-8 21:45:11 | 显示全部楼层
老老实实用XYZ坐标不好么。。内坐标下面那些变量是要有数据的,从这个来看你一点都不懂内坐标,而且最后报错信息写得很清楚了,后面需要一个浮点数,你却啥都没有

3

帖子

0

威望

37

eV
积分
40

Level 2 能力者

 楼主| 发表于 2018-10-9 09:06:19 | 显示全部楼层
niobium 发表于 2018-10-8 21:45
老老实实用XYZ坐标不好么。。内坐标下面那些变量是要有数据的,从这个来看你一点都不懂内坐标,而且最后报 ...

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

1万

帖子

25

威望

1万

eV
积分
34922

管理员

公社社长

发表于 2018-10-9 18:26:23 | 显示全部楼层
FANYI 发表于 2018-10-9 09:06
非常谢谢您的回复~!请问用直角坐标的话,虚原子也可以通过添加-1 来固定其位置吗?然后,这个内坐标的方 ...

固定虚原子毫无意义,本来也不影响结果
仔细看
在Gaussian中做限制性优化的方法
http://sobereva.com/404
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)
计算化学公社论坛:http://bbs.keinsci.com(高水平、高人气、综合性计算化学交流论坛)
思想家公社QQ群1号:18616395,2号:466017436。用于讨论理论、计算化学,两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。研究方向和理论、计算化学无关者勿加,以免浪费宝贵的空位

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

3

帖子

0

威望

37

eV
积分
40

Level 2 能力者

 楼主| 发表于 2018-10-11 11:11:34 | 显示全部楼层
sobereva 发表于 2018-10-9 18:26
固定虚原子毫无意义,本来也不影响结果
仔细看
在Gaussian中做限制性优化的方法

非常感谢社长的回复与指点!!
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

手机版|北京科音自然科学研究中心|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949-1号 )

GMT+8, 2018-10-20 21:58 , Processed in 0.205444 second(s), 24 queries .

快速回复 返回顶部 返回列表