计算化学公社

标题: 计算激发态优化结束结构有个碳氢键断裂，为什么？ [打印本页]

作者
Author: 大肥猫 时间: 2021-4-30 14:39
标题: 计算激发态优化结束结构有个碳氢键断裂，为什么？
本帖最后由大肥猫于 2021-4-30 14:44 编辑

在对附件结构进行激发态优化，计算激发能时，激发态优化正常结束，但是得到的结构其中一个碳氢键断裂了（见下图），请问是什么原因呢？或者说明了什么呢？说明在激发态这个断裂的质子比较活泼吗？
输入文件如下：

%nprocshared=14
%mem=15GB
%chk=YHWA-24-search-1.1.chk
# opt pbe1pbe/6-311g(d) TD scrf=(smd,solvent=acetonitrile)

D:\EtOAc-calculate\YHWA-24\cluster.xyz

0 1
C                2.85016500 1.01953700 0.36661200
C                3.06315000 -0.22483500 -0.10265200
C                1.92048800 -1.18136800 -0.26875300
C                0.54587200 -0.71536400 0.05188500
C                0.33669900 0.54438800 0.49528000
C                1.51393700 1.45453600 0.73325500
O                1.34308800 2.55246900 1.24333300
O                2.14330900 -2.30656100 -0.66637400
C                -0.53921600 -1.76495400 -0.14274400
C                -1.68496900 -1.28158800 -0.97303600
C                -2.81627200 -1.33555200 -0.27482100
C                -4.09455500 -0.63385200 -0.55907000
C                -4.01705300 0.78046000 0.07844800
C                -1.00407000 1.14712400 0.83460300
C                -1.71203700 1.81874400 -0.33014500
C                -3.01638200 1.66416400 -0.59933800
C                -0.87218800 2.74679500 -1.15805100
C                -2.53537900 -1.96714100 1.05362300
O                -1.12199800 -2.12819500 1.11305700
O                4.22159300 -0.75929400 -0.45101600
C                5.38273100 0.05588700 -0.32785300
H                3.64668400 1.73624100 0.52731700
H                -0.05613300 -2.65150100 -0.56777900
H                -1.55446700 -0.86601900 -1.96529300
H                -4.26417800 -0.54101500 -1.63588600
H                -4.94955500 -1.17005400 -0.13349500
H                -3.80607700 0.67504400 1.14731500
H                -5.00786900 1.24086500 0.00518900
H                -1.64057700 0.39788000 1.29775600
H                -0.80853200 1.91978800 1.58693400
H                -3.40287600 2.22149700 -1.45391300
H                -0.04545200 2.21367700 -1.64389400
H                -1.46409600 3.22885300 -1.93946400
H                -0.41418900 3.52693100 -0.54048800
H                -3.02530100 -2.94729100 1.15027900
H                -2.87176300 -1.35161400 1.89982200
H                5.53610200 0.35441200 0.71279100
H                6.21430200 -0.56066600 -0.66208100
H                5.30050400 0.94372500 -0.96056700

作者
Author: wzkchem5 时间: 2021-4-30 16:15
仔细读这篇文章 http://sobereva.com/414

作者
Author: sobereva 时间: 2021-4-30 20:51
可以和相同泛函和基组下优化的基态的键长对比，如果有可查觉的增长，暗示激发态的时候这个键被削弱了

作者
Author: 大肥猫 时间: 2021-4-30 21:16

sobereva 发表于 2021-4-30 20:51
可以和相同泛函和基组下优化的基态的键长对比，如果有可查觉的增长，暗示激发态的时候这个键被削弱了

老师，我看了，在基态优化后键长为1.09544，激发态同等计算水平键长为1.13706，是否说明在激发态该位点键长弱化，氢质子反应性更活泼了呢？在实验中，该位点的质子发生了转移。

作者
Author: sobereva 时间: 2021-5-1 22:48

大肥猫发表于 2021-4-30 21:16
老师，我看了，在基态优化后键长为1.09544，激发态同等计算水平键长为1.13706，是否说明在激发态该位点键 ...

可以

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