计算化学公社

标题: 氯苯在CuCl2表面的催化反应计算出错 [打印本页]
作者
Author: 升腾    时间: 2016-5-25 21:46
标题: 氯苯在CuCl2表面的催化反应计算出错
请问大神：我想对一个氯苯与CuCl2进行直接优化到过渡态，在523k,1atm,下精心计算。我想问一下，输入文件怎么设置，我的计算的输入文件如下，老是出现错误。
%chk=C:\Users\Administrator\Desktop\gaosi\opt+ts=lvben+CuCl2.chk
#p opt=(calcfc,noeigentest,ts) b3lyp/genecp  Freq=Readisotopes
opt+ts=lvben+CuCl2
0 2
C                 -1.83118742    0.40772532    0.00000000
C                 -1.13209642    1.61890632    0.00000000
C                  0.26786658    1.62822232   -0.00000100
C                  0.93198558    0.40772532   -0.00000200
C                  0.26786658   -0.81277068    0.00000000
C                 -1.13209542   -0.80345568    0.00000000
H                 -2.91197542    0.40772432    0.00000100
H                 -1.66712842    2.55812732    0.00000100
H                  0.81628458    2.55783432    0.00000100
H                  0.81628558   -1.74238268   -0.00000100
H                 -1.66712742   -1.74267768    0.00000100
Cl                 2.77678258    0.40772532    0.00000000
Cu                -0.29897736   -3.33132729    0.29068087
Cl                -2.33371097   -3.65168894    0.33075915
Cl                 1.73575625   -3.01096565    0.25060259
C H Cl 0
6-31g
****
Cu 0
SDD
****
Cu 0
SDD

523 1.0
12
12
12
12
12
12
1
1
1
1
1
35
64
35
35

作者
Author: liyuanhe211    时间: 2016-5-26 00:34
本帖最后由 liyuanhe211 于 2016-5-26 03:17 编辑

错误要说明是什么错误。
试试这个，没测试但差不多了，精度仍然不高但显然比你原来的好。改它的动力是“在523k,1atm,下【精心计算】”，23333，对计算姬充满爱就不报错了。
  1. #p opt=(calcfc,noeigentest,ts) b3lyp/genecp Temperature=523

  3. opt+ts=lvben+CuCl2

  5. 0 2
  6. C -1.83118742 0.40772532 0.00000000
  7. C -1.13209642 1.61890632 0.00000000
  8. C 0.26786658 1.62822232 -0.00000100
  9. C 0.93198558 0.40772532 -0.00000200
  10. C 0.26786658 -0.81277068 0.00000000
  11. C -1.13209542 -0.80345568 0.00000000
  12. H -2.91197542 0.40772432 0.00000100
  13. H -1.66712842 2.55812732 0.00000100
  14. H 0.81628458 2.55783432 0.00000100
  15. H 0.81628558 -1.74238268 -0.00000100
  16. H -1.66712742 -1.74267768 0.00000100
  17. Cl 2.77678258 0.40772532 0.00000000
  18. Cu -0.29897736 -3.33132729 0.29068087
  19. Cl -2.33371097 -3.65168894 0.33075915
  20. Cl 1.73575625 -3.01096565 0.25060259

  22. C H Cl 0
  23. 6-31g(d)
  24. ****
  25. Cu 0
  26. LanL2DZ
  27. ****

  29. Cu 0
  30. LanL2DZ
复制代码


作者
Author: sobereva    时间: 2016-5-26 00:52
6-31G这么烂基组算出来结果没意义
注意正确空行，如果你把文件内容贴出来的时候论坛自动把一些空行去掉了，你就直接上传输入文件，免得令人误会。

如果你只想相对于默认设定改个温度，用temperature=523就完了，没必要用readisotopes写那么一大串。
作者
Author: 升腾    时间: 2016-5-26 20:07
Sobereva说的对，弟子谨遵教诲！
作者
Author: 升腾    时间: 2016-5-26 20:58
liyuanhe211 发表于 2016-5-26 00:34
错误要说明是什么错误。
试试这个，没测试但差不多了，精度仍然不高但显然比你原来的好。改它的动力是“在 ...

嗯嗯，O(∩_∩)O谢谢，刚刚开始接触高斯，遇到各种问题，望各位海涵！



欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3