计算化学公社

标题: orca过渡态问题 [打印本页]

作者
Author: dingliping1108 时间: 2016-5-6 15:06
标题: orca过渡态问题
请问各位老师，我用orca寻找过渡态得到的过渡态（一个虚频）的电子能final single point energy （都是负值)，比两个反应物各自的电子能（负值）都低，两个反应物优化后的电子能分别约-57，-1917，而过渡态的电子能约是-1973，是我的计算过程中有什么错误吗？http://bbs.keinsci.com/forum.php ... =%B9%FD%B6%C9%CC%ACsob老师这个帖子里说是漏掉了中间体。我怎么改正呢？
下面这个是寻找过渡态的命令。

%pal nprocs 4
  end

! RHF OptTS TightSCF SlowConv
! NumFreq
! COSMO(THF)

%basis
Basis _6_31G
end

%geom Calc_Hess true; NumHess true

  end

下面这个是结构优化的命令
%pal nprocs 4
end

! RHF Opt TightSCF SlowConv
! NumFreq
! COSMO(THF)

%basis
Basis _6_31G
end

%geom Calc_Hess true; NumHess true

  end

作者
Author: liyuanhe211 时间: 2016-5-6 16:45
两个反应物的能量是-57和-1917，所以总能量是-1974，你过渡态是-1973，并不比反应物低。

命令写的太啰嗦，例如6-31G可以直接写，基组6-31G太烂，有解析不要用数值。

作者
Author: dingliping1108 时间: 2016-5-7 13:53

liyuanhe211 发表于 2016-5-6 16:45
两个反应物的能量是-57和-1917，所以总能量是-1974，你过渡态是-1973，并不比反应物低。

命令写的太啰嗦 ...

我把Numfreq命令去掉之后计算一直报错，所以就带着他算了，那请问您我用哪个基组比较好呢？

作者
Author: liyuanhe211 时间: 2016-5-7 14:50

dingliping1108 发表于 2016-5-7 13:53
我把Numfreq命令去掉之后计算一直报错，所以就带着他算了，那请问您我用哪个基组比较好呢？

这是练习的时候的一个输入，没用数值并没有报错。不清楚你的是怎么回事。

! B3LYP 6-31G(d) TightSCF OptTs freq
%geom calc_hess true end

*xyz -1 1
C       -1.55384    -0.43723       0.00000
Cl       -0.48384    -0.43723       0.00000
H       -1.91051    -1.20940       0.64919
H       -1.91051    -0.61336    -0.99331
H       -1.91051       0.51106       0.34413
F       -4.54939    -0.44897       0.08390
*

作者
Author: dingliping1108 时间: 2016-5-8 17:36

liyuanhe211 发表于 2016-5-7 14:50
这是练习的时候的一个输入，没用数值并没有报错。不清楚你的是怎么回事。

! B3LYP 6-31G(d) TightSCF ...

我試了，沒有報錯，謝謝liyuanhe211老師

作者
Author: dingliping1108 时间: 2016-5-10 09:20

liyuanhe211 发表于 2016-5-7 14:50
这是练习的时候的一个输入，没用数值并没有报错。不清楚你的是怎么回事。

! B3LYP 6-31G(d) TightSCF ...

  您好，我按您给的例子试了一下，为什么结果里面虚频这么多，这么大呀？
(, 下载次数 Times of downloads: 5)

(, 下载次数 Times of downloads: 3)

0:       0.00 cm**-1
1:       0.00 cm**-1
2:       0.00 cm**-1
3:       0.00 cm**-1
4:       0.00 cm**-1
5:       0.00 cm**-1
6: -731304.34 cm**-1 ***imaginary mode***
7: -588698.46 cm**-1 ***imaginary mode***
8: -55062.11 cm**-1 ***imaginary mode***
9: -250527.65 cm**-1 ***imaginary mode***
  10: -147673.90 cm**-1 ***imaginary mode***
  11: -108016.02 cm**-1 ***imaginary mode***
  12:    6895.77 cm**-1
  13:    10983.89 cm**-1
  14:    42312.88 cm**-1
  15:    46803.39 cm**-1
  16: 146859.94 cm**-1
  17: 205033.29 cm**-1

作者
Author: plus 时间: 2016-5-10 11:45
WARNING!!!!!!!
The optimization did not converge but reached the maximum number of
optimization cycles.
Nevertheless we allow to calculate frequencies and other properties at this point.

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