关于光敏剂激发态计算的能量的问题

蛋蛋

各位老师好，在计算光敏剂催化单电子转移的能量的时候，我需要计算激发态的能量，拿到激发态的吉布斯自由能计算反应热。在学习计算的过程中，我还是不能确认自己的计算步骤对不对，我将自己的学习过程总结如下，请各位老师指正。（1）这个是我的光敏剂的结构Ru*，实验人提供的最大吸收波长是452nm。


（2）在不知道是哪个激发态起作用的时候，先对基态结构进行优化，再默认计算第一激发态，算了十个态。（疑问：想着基态优化过了，我只想看看激发态的能量就没有加opt ,但是这样找到的Excited State10这个态设置对吗？）
输入：

输出：



（3）从输出文件看到Excited State10的吸收峰 449nm，接近Ru*光敏剂的452nm, 接下来对第十个激发态重点优化，输入文件如下：（疑问：不知道这样的设置是不是算出来的刚好是第十个激发态的能量）



（4）等第十个态的计算完毕，都取check文件再换大基组做单点计算，输入如下：



不知道我的计算过程是否正确，这样拿到的吉布斯自由能再经溶剂矫正是否可以使用？请各位老师指点。

biogon

先不说别的，T10振子强度是0凭什么认为是它发光，而且发光也不可能从这么高激发态来的，2之后完全就是瞎算

tiandikuoyuan

triplet振子强度默认是0，一般认为吸收主要是单线态吸收；你这个应该直接优化T1结构，钌配合物敏化的时候应该是处于最低激发三重态上，个人观点。

蛋蛋

tiandikuoyuan 发表于 2020-12-25 16:12
triplet振子强度默认是0，一般认为吸收主要是单线态吸收；你这个应该直接优化T1结构，钌配合物敏化的时候应 ...

我最开始是这么算的，然后有个师兄说不一定是T1,要我找到跟它最大吸收波长对应的态，拿到按个态的能量才可以。请问你知道具体怎么做嘛，教教我

wzkchem5

你的分子画了氢原子吗？是把氢原子全部隐藏了还是没有画氢原子？
另外初始结构的坐标单位（埃/Bohr）对吗？

蛋蛋

biogon 发表于 2020-12-25 15:55
先不说别的，T10振子强度是0凭什么认为是它发光，而且发光也不可能从这么高激发态来的，2之后完全就是瞎算

我后面对激发态进行优化的时候贞子强度就不是0 了。请问具体应该是什么步骤啊，麻烦你说的详细一点，我好再去试一下

蛋蛋

wzkchem5 发表于 2020-12-26 09:03
你的分子画了氢原子吗？是把氢原子全部隐藏了还是没有画氢原子？
另外初始结构的坐标单位（埃/Bohr）对吗 ...

我画了，都隐藏了，我怕看起来乱

蛋蛋

蛋蛋 发表于 2020-12-26 10:01
我画了，都隐藏了，我怕看起来乱

这个是我从文献里粘贴出来的坐标，应该是没问题的

蛋蛋

wzkchem5 发表于 2020-12-26 09:03
你的分子画了氢原子吗？是把氢原子全部隐藏了还是没有画氢原子？
另外初始结构的坐标单位（埃/Bohr）对吗 ...

坐标应该没啥问题，是我从文献里粘贴出来的

wzkchem5

蛋蛋 发表于 2020-12-26 10:03
坐标应该没啥问题，是我从文献里粘贴出来的

粘贴出来也要确认一下单位是对的，如果文献的单位是bohr，你输入文件的单位是埃，就错了

蛋蛋

wzkchem5 发表于 2020-12-26 13:38
粘贴出来也要确认一下单位是对的，如果文献的单位是bohr，你输入文件的单位是埃，就错了

计算有机反应机理的SI里所有的坐标都是高斯导出来的呀，单位都是一致的 （埃），不会特意改变单位的。