计算化学公社

标题: 高斯能计算激发态的过渡态么？ [打印本页]

作者
Author: caocao321 时间: 2014-12-28 18:11
标题: 高斯能计算激发态的过渡态么？
如题，高斯怎么直接计算激发态的过渡态？能计算么？

作者
Author: sobereva 时间: 2014-12-28 18:29
激发态和基态在寻找过渡态的算法上没有差异，关键词也一样。
CIS有二阶解析导数，就照常用过渡态的搜索关键词即可。
TDDFT只有一阶解析导数，需要用opt=(TS,modRedundant,noeigen)或者opt=(TS,gediis,noeigen)避免使用二阶导数，否则得用readfc来读取单独计算的TDDFT的Hessian矩阵。

作者
Author: caocao321 时间: 2014-12-29 10:26

sobereva 发表于 2014-12-28 18:29
激发态和基态在寻找过渡态的算法上没有差异，关键词也一样。
CIS有二阶解析导数，就照常用过渡态的搜索关 ...

谢了sob。

作者
Author: caocao321 时间: 2014-12-30 20:23
本帖最后由 caocao321 于 2014-12-30 20:26 编辑

sobereva 发表于 2014-12-28 18:29
激发态和基态在寻找过渡态的算法上没有差异，关键词也一样。
CIS有二阶解析导数，就照常用过渡态的搜索关 ...

还想再问下，高斯中得到了TDDFT下的过渡态，算出了频率。怎么样才能得到TDDFT下的IRC？貌似算不了的样子。

作者
Author: sobereva 时间: 2014-12-31 05:21
关键词和基态一样。但由于TD没有解析二阶导数，如手册中指出的，此时必须用GradientOnly关键词。如果结果异常的话，用IRC(GradientOnly,euler)代替GradientOnly时默认的DVV方法。

作者
Author: caocao321 时间: 2014-12-31 16:47

sobereva 发表于 2014-12-31 05:21
关键词和基态一样。但由于TD没有解析二阶导数，如手册中指出的，此时必须用GradientOnly关键词。如果结果异 ...

谢谢了。不过我的计算又出现问题了。我把输出文件贴出来，你再帮我看看吧。频率计算出只有一个虚频，可是这个输出文件中给出了warning：more than one imaginary fre.
  ----------------------------------------------------------------------
#p geom=check guess=read irc=(reverse,RCFC,GradientOnly,euler,maxpoint
=100,stepsize=5) td=(singlets,root=1) cam-b3lyp/6-31g(d)

INPUT DATA FOR L123
------------------------------------------------------------------------
GENERAL PARAMETERS:
Rxn path following direction = Reverse
Maximum points per path    = 100
Step size                   = 0.050 bohr
Integration scheme          = Euler
Hessian evaluation          = Gradient only (Hessian at TS only).
---------------------------------------------------------------------
----------------------------------------------------------------------
      ******** Start new reaction path calculation ********
RCFC Option Requested - Data Read From Chk File:
ct-ts-r.chk
Energy From Chk = -759.1681645

Current Structure is TS -> diagonalize the Hessian.
WARNING: TS HAS MORE THAN 1 IMAGINARY FREQUENCY!
Supplied step size of 0.0500 bohr.
Integration on MW PES will use step size of    NaN sqrt(amu)*bohr.
Point Number: 0       Path Number: 1
Calculating another point on the path.
Point Number  1 in REVERSE path direction.
Using Euler Reaction Path Following.
Gradient is small - USING TRANSITION VECTOR.
IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC
Leave Link  123 at Wed Dec 31 10:19:24 2014, MaxMem=  536870912 cpu:    0.0
(Enter /share/apps/gauss/g09/l202.exe)
                     Input orientation:
---------------------------------------------------------------------
Center    Atomic    Atomic          Coordinates (Angstroms)
Number    Number    Type          X          Y          Z
---------------------------------------------------------------------
   1       6          0          NaN       NaN       NaN
   2       6          0          NaN       NaN       NaN
   3       7          0          NaN       NaN       NaN
   4       7          0          NaN       NaN       NaN
   5       1          0          NaN       NaN       NaN
   6       1          0          NaN       NaN       NaN
   7       7          0          NaN       NaN       NaN
   8       7          0          NaN       NaN       NaN
   9       6          0          NaN       NaN       NaN
   10       6          0          NaN       NaN       NaN
   11       6          0          NaN       NaN       NaN
   12       6          0          NaN       NaN       NaN
   13       6          0          NaN       NaN       NaN
   14       6          0          NaN       NaN       NaN
   15       6          0          NaN       NaN       NaN
   16       6          0          NaN       NaN       NaN
   17       6          0          NaN       NaN       NaN
   18       6          0          NaN       NaN       NaN
   19       6          0          NaN       NaN       NaN
   20       6          0          NaN       NaN       NaN
   21       1          0          NaN       NaN       NaN
   22       1          0          NaN       NaN       NaN
   23       1          0          NaN       NaN       NaN
   24       1          0          NaN       NaN       NaN
   25       1          0          NaN       NaN       NaN
   26       1          0          NaN       NaN       NaN
   27       1          0          NaN       NaN       NaN
   28       1          0          NaN       NaN       NaN
   29       1          0          NaN       NaN       NaN
   30       1          0          NaN       NaN       NaN
---------------------------------------------------------------------
               Distance matrix (angstroms):
                  1       2       3       4       5
   1  C 0.000000
   2  C 0.000000 0.000000
   3  N 0.000000 0.000000 0.000000
   4  N 0.000000 0.000000 0.000000 0.000000
   5  H 0.000000 0.000000 0.000000 0.000000 0.000000
   6  H 0.000000 0.000000 0.000000 0.000000 0.000000
   7  N 0.000000 0.000000 0.000000 0.000000 0.000000
   8  N 0.000000 0.000000 0.000000 0.000000 0.000000
   9  C 0.000000 0.000000 0.000000 0.000000 0.000000
10  C 0.000000 0.000000 0.000000 0.000000 0.000000
11  C 0.000000 0.000000 0.000000 0.000000 0.000000
12  C 0.000000 0.000000 0.000000 0.000000 0.000000
13  C 0.000000 0.000000 0.000000 0.000000 0.000000
14  C 0.000000 0.000000 0.000000 0.000000 0.000000
15  C 0.000000 0.000000 0.000000 0.000000 0.000000
16  C 0.000000 0.000000 0.000000 0.000000 0.000000
17  C 0.000000 0.000000 0.000000 0.000000 0.000000
18  C 0.000000 0.000000 0.000000 0.000000 0.000000
19  C 0.000000 0.000000 0.000000 0.000000 0.000000
20  C 0.000000 0.000000 0.000000 0.000000 0.000000
21  H 0.000000 0.000000 0.000000 0.000000 0.000000
22  H 0.000000 0.000000 0.000000 0.000000 0.000000
23  H 0.000000 0.000000 0.000000 0.000000 0.000000
24  H 0.000000 0.000000 0.000000 0.000000 0.000000
25  H 0.000000 0.000000 0.000000 0.000000 0.000000
26  H 0.000000 0.000000 0.000000 0.000000 0.000000
27  H 0.000000 0.000000 0.000000 0.000000 0.000000
28  H 0.000000 0.000000 0.000000 0.000000 0.000000
29  H 0.000000 0.000000 0.000000 0.000000 0.000000
30  H 0.000000 0.000000 0.000000 0.000000 0.000000
                  6       7       8       9       10
   6  H 0.000000
.........
Small interatomic distances encountered:    2    1    3    1    3    2    4    1    4    2
Small interatomic distances encountered:    4    3    5    1    5    2    5    3    5    4
Small interatomic distances encountered:    6    1    6    2    6    3    6    4    6    5
Small interatomic distances encountered:    7    1    7    2    7    3    7    4    7    5
Problem with the distance matrix.
Error termination via Lnk1e in /share/apps/gauss/g09/l202.exe at Wed Dec 31 16:40:53 2014.
Job cpu time:  0 days  0 hours  0 minutes  1.0 seconds.
File lengths (MBytes):  RWF=    9 Int=    0 D2E=    0 Chk=    75 Scr=    1

作者
Author: sobereva 时间: 2014-12-31 17:01
先确认chk里的TS结构是否确实已经在cam-b3lyp/6-31g(d) td=(singlets,root=1)下优化好了。并且建议去掉geom=check和guess=read，直接写坐标，尽量简化关键词多尝试才容易找出原因。
RCFC没必要写。

作者
Author: caocao321 时间: 2015-1-1 10:30

sobereva 发表于 2014-12-31 17:01
先确认chk里的TS结构是否确实已经在cam-b3lyp/6-31g(d) td=(singlets,root=1)下优化好了。并且建议去掉geom ...

确实是关键词给的太多了，删了RCFC，给了几何坐标就可以算了。谢谢sob了，受益匪浅。

作者
Author: Regina 时间: 2019-4-23 10:28
请问一下，寻找过渡态的时候需要固定键长吗？还是说直接寻找就可以，我直接寻找找到的不是我想要的，但是固定键长寻找就会有两个虚频

作者
Author: sobereva 时间: 2019-4-24 03:13

Regina 发表于 2019-4-23 10:28
请问一下，寻找过渡态的时候需要固定键长吗？还是说直接寻找就可以，我直接寻找找到的不是我想要的，但是固 ...

这主贴关系不密切的问题不要在回复里提问
简单问题有专门的简单量化问题答疑专帖

固定键长必然得不到严格的过渡态结构

作者
Author: Regina 时间: 2019-4-24 08:53

sobereva 发表于 2019-4-24 03:13
这主贴关系不密切的问题不要在回复里提问
简单问题有专门的简单量化问题答疑专帖

好的，谢谢老师

作者
Author: 上善若水 时间: 2022-8-16 21:59

sobereva 发表于 2014-12-31 17:01
先确认chk里的TS结构是否确实已经在cam-b3lyp/6-31g(d) td=(singlets,root=1)下优化好了。并且建议去掉geom ...

请问您卢老师，1）使用Gaussian 09在激发态过渡态走IRC(gradientonly,euler)时，我看您的帖子帖子说成功率极低，能否用readfc方法替代这个关键词呢？
2）如果用TS方法这个read的CHK文件是否是摆好的结构的激发态的chk文件呢？
3）激发态过渡态的计算，是否可以理解为反应物为S1（或者其他激发态）的分子到S1态的产物的S1过渡态呢？

作者
Author: sobereva 时间: 2022-8-17 17:42

上善若水发表于 2022-8-16 21:59
请问您卢老师，1）使用Gaussian 09在激发态过渡态走IRC(gradientonly,euler)时，我看您的帖子帖子说成功 ...

1、2 尽量用G16结合calcfc。就算readfc，你也得算激发态Hessian，耗时还是很高
3 就是你研究的那个激发态势能面上的一阶鞍点

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