计算化学公社

 找回密码 Forget password
 注册 Register
Views: 11854|回复 Reply: 4
打印 Print 上一主题 Last thread 下一主题 Next thread

[新手求助] 氯苯在CuCl2表面的催化反应计算出错

[复制链接 Copy URL]

27

帖子

0

威望

127

eV
积分
154

Level 3 能力者

跳转到指定楼层 Go to specific reply
楼主
请问大神:我想对一个氯苯与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

3097

帖子

29

威望

1万

eV
积分
17098

Level 6 (一方通行)

2#
发表于 Post on 2016-5-26 00:34:53 | 只看该作者 Only view this author
本帖最后由 liyuanhe211 于 2016-5-26 03:17 编辑

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

  2. opt+ts=lvben+CuCl2

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

  19. C H Cl 0
  20. 6-31g(d)
  21. ****
  22. Cu 0
  23. LanL2DZ
  24. ****

  25. Cu 0
  26. LanL2DZ
复制代码


6万

帖子

99

威望

5万

eV
积分
120165

管理员

公社社长

3#
发表于 Post on 2016-5-26 00:52:05 | 只看该作者 Only view this author
6-31G这么烂基组算出来结果没意义
注意正确空行,如果你把文件内容贴出来的时候论坛自动把一些空行去掉了,你就直接上传输入文件,免得令人误会。

如果你只想相对于默认设定改个温度,用temperature=523就完了,没必要用readisotopes写那么一大串。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

27

帖子

0

威望

127

eV
积分
154

Level 3 能力者

4#
 楼主 Author| 发表于 Post on 2016-5-26 20:07:35 | 只看该作者 Only view this author
Sobereva说的对,弟子谨遵教诲!

27

帖子

0

威望

127

eV
积分
154

Level 3 能力者

5#
 楼主 Author| 发表于 Post on 2016-5-26 20:58:45 | 只看该作者 Only view this author
liyuanhe211 发表于 2016-5-26 00:34
错误要说明是什么错误。
试试这个,没测试但差不多了,精度仍然不高但显然比你原来的好。改它的动力是“在 ...

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

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2025-8-16 18:41 , Processed in 0.151328 second(s), 20 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list