氯苯在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

liyuanhe211

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

1. #p opt=(calcfc,noeigentest,ts) b3lyp/genecp Temperature=523
   
2. 
   
3. opt+ts=lvben+CuCl2
   
4. 
   
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
    
21. 
    
22. C H Cl 0
    
23. 6-31g(d)
    
24. ****
    
25. Cu 0
    
26. LanL2DZ
    
27. ****
    
28. 
    
29. Cu 0
    
30. LanL2DZ
    

复制代码

sobereva

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

如果你只想相对于默认设定改个温度，用temperature=523就完了，没必要用readisotopes写那么一大串。

升腾

Sobereva说的对，弟子谨遵教诲！

升腾

liyuanhe211 发表于 2016-5-26 00:34
错误要说明是什么错误。
试试这个，没测试但差不多了，精度仍然不高但显然比你原来的好。改它的动力是“在 ...

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