请选择 进入手机版 | 继续访问电脑版

计算化学公社

 找回密码
 现在注册!
查看: 153|回复: 4

[Gaussian/gview] 求助:振动分辨的电子光谱计算出错

[复制链接]

29

帖子

0

威望

96

eV
积分
125

Level 2 能力者

发表于 2019-8-11 09:24:18 | 显示全部楼层 |阅读模式
本帖最后由 liaorongbao 于 2019-8-11 09:43 编辑

Sob老师的博文http://sobereva.com/223
第一步:进行激发态计算如下
%rwf=SSexc.rwf
%nosave
%mem=30GB
%nprocshared=30
%chk=SSexc.chk
# maxdisk=500GB
#P wB97XD/6-31+G* TD opt int=ultrafine freq=saveNM  scrf=(solvent=acetonitrile)

it is already opted at wB97XD/6-31+G* scrf=(solvent=acetonitrile)

0 1
C                 -2.27482700   -0.67684600    0.13841000
C                 -3.37666400   -1.59428000   -0.14148200
O                 -3.28443300   -2.60572800   -0.82584500
O                 -4.54114700   -1.22585300    0.43084400
N                 -2.75174000    1.28735400    1.74718600
C                 -2.53455900    0.41319300    1.00708300
C                 -0.99077400   -1.03422800   -0.30606700
H                 -0.94839000   -2.06160800   -0.65625000
C                  0.21341600   -0.30230200   -0.27605500
C                  1.48777100   -1.01108100   -0.15095100
C                  2.64770800   -0.27603200    0.25484700
C                  2.57907100    1.12093600    0.43987100
H                  3.46787500    1.65683200    0.76410800
C                  1.42349200    1.84596100    0.10329000
C                  0.25068800    1.15188900   -0.33580600
C                  1.62037100   -2.39257000   -0.34678900
H                  0.77380300   -2.99082200   -0.66271900
C                  2.84169700   -3.04258500   -0.16119900
H                  2.90575900   -4.11281100   -0.32936000
C                  3.96852200   -2.32658700    0.24640500
H                  4.91540900   -2.83378100    0.40025200
C                  3.86831300   -0.95902000    0.45594700
H                  4.73791600   -0.39144500    0.77580500
C                  1.42741600    3.26019000    0.11079100
H                  2.30077300    3.78010100    0.49487400
C                  0.35144000    3.97950600   -0.39128600
H                  0.37747600    5.06439200   -0.39343500
C                 -0.74817700    3.29992600   -0.91885600
H                 -1.57640800    3.85214300   -1.35083900
C                 -0.79329700    1.90438600   -0.88937600
H                 -1.64525200    1.39991000   -1.33231800
H                 -5.21177300   -1.88447100    0.18691800

第二步:进行基态计算如下
%rwf=SSUVgr.rwf
%nosave
%mem=30GB
%nprocshared=30
%chk=SSUVgr.chk
# maxdisk=500GB
#P wB97XD/6-31+G* OPT freq=FC scrf=(solvent=acetonitrile) nosymm

it is already opted at wB97XD/6-31G* scrf=(solvent=acetonitrile)

0 1
分子结构说明

SSexc.chk


第三步:打开基态计算后的结果,发现出错情况如下:
     ==================================================
                 Duschinsky Transformation
     ==================================================

     ==================================================
              Calculations of Band Intensities
     ==================================================

Energy of the 0-0 transition:  19884.5994 cm^(-1)


Notes about the overlap integrals description:
  - DeltaE is the relative energy of the transition (wrt 0-0)
  - TDMI^2 is the square of the transition dipole moment integral
  - Intensity is the line intensity

ERROR: The Franck-Condon factor corresponding to the overlap integral
        between both vibrational ground states is too small:
         |<0'|0">|^2 =  0.13404E-08
Error termination via Lnk1e in /soft/g09/l718.exe at Sun Aug 11 17:48:54 2019.
Job cpu time:       1 days  3 hours 31 minutes 50.3 seconds.
File lengths (MBytes):  RWF=   1187 Int=      0 D2E=      0 Chk=     16 Scr=      1
为什么出错?|<0'|0">|^2 =  0.13404E-08是不是意味着没有紫外吸收光谱?如果0-0跃迁不行的话可以接着给出其他跃迁呀





SSUVgr.gjf

2.12 KB, 下载次数: 0

基态输入文件

1万

帖子

25

威望

2万

eV
积分
44706

管理员

公社社长

发表于 2019-8-11 10:08:57 | 显示全部楼层
这部分代码开发者的文档里的答案
This error indicates that it is very likely that the underlying theory used to generate the spectrum, the Franck-
Condon principle, is unsuitable for the case of interest. This is a typical issue when the equilibrium geometry of the
initial and final states are noticeably different. This can happen when one of the structures is planar in one electronic
state and not planar in the other electronic state.
Another indicator of this problem is the Duschinsky matrix (which can be displayed with the option “PRTMAT=1”
in the ReadFCHT section) that shows the projection of the normal modes of one state on those of the other state. In
this case, the Duschinsky matrix is rather diffuse, that is to say that many normal coordinates of the final state are
expressed by a combination of a rather large number of normal modes of the initial state.
If the calculations must be carried anyway, this error can be bypassed with the keyword ForcedFCCalc. However,
please note that the resulting spectrum will be likely unreliable and the corresponding calculations should be
considered as tests.

评分

参与人数 1eV +4 收起 理由
liaorongbao + 4 正解

查看全部评分

北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社QQ群,1号:18616395,2号:466017436。达5000人,专门交流理论、计算化学。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

29

帖子

0

威望

96

eV
积分
125

Level 2 能力者

 楼主| 发表于 2019-8-11 11:28:42 | 显示全部楼层
社长真细心!这都发现了。俺用 ForcedFCCalc看看可行,反正就是走个形式

8

帖子

0

威望

178

eV
积分
186

Level 3 能力者

发表于 7 天前 | 显示全部楼层
社长对于该博文的前面几句话,请教社长,该方法在使用G16时会有更好的方法吗?如何实现?

29

帖子

0

威望

96

eV
积分
125

Level 2 能力者

 楼主| 发表于 7 天前 | 显示全部楼层
comboy2008 发表于 2019-8-14 17:41
社长对于该博文的前面几句话,请教社长,该方法在使用G16时会有更好的方法吗?如何实现?

感觉没读懂你问的是什么问题,尽管很关注。不知社长能否读懂
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

手机版|北京科音自然科学研究中心|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949-1号 )

GMT+8, 2019-8-21 20:45 , Processed in 0.166420 second(s), 28 queries .

快速回复 返回顶部 返回列表