orca过渡态问题

dingliping1108

请问各位老师，我用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

liyuanhe211

两个反应物的能量是-57和-1917，所以总能量是-1974，你过渡态是-1973，并不比反应物低。

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

dingliping1108

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

命令写的太啰嗦 ...

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

liyuanhe211

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
*

dingliping1108

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

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

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

dingliping1108

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

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

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

2016-5-10 09:18 上传 Uploaded

点击下载
Click to download

test.out (806.44 KB, 下载次数 Times of downloads: 3)

2016-5-10 09:18 上传 Uploaded

点击下载
Click to download

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

plus

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.