计算化学公社

标题: 如何计算金属配合物的三重态能级 [打印本页]

作者
Author: 哈哈宝 时间: 2015-2-9 17:45
标题: 如何计算金属配合物的三重态能级
由于实验需要，我希望计算一个含钌过渡金属配合物的三重态能级。
请问各位老师，我的计算方法如下，是否正确，该如何分析？
用B3LYP/6-31g*+LanL2DZ优化了基态几何结构
用UB3LYP/6-31g*+LanL2DZ opt 继续优化激发态几何结构
用td(triplets, nstats=15, root=1) 做发射光谱计算，由于禁阻，振子强度都是零
我观察输出文件中Excited 1，其激发能为1.86eV，而波长在1180 nm，这个数据（1.89eV）就是三重态能级么？可是磷光发光位置似乎不应该到1180左右啊？

希望能得到老师们的解惑，谢谢。

作者
Author: sobereva 时间: 2015-2-9 18:12
如果你是指0 3下用UB3LYP/6-31g*+LanL2DZ opt优化，在此结构基础上，又在0 1下做td(triplets)得到S0和T1的能量差1.86eV，从计算步骤上来说并没错误。

你的单位换算是不是搞错了？1.89eV对应于667nm才对

作者
Author: 哈哈宝 时间: 2015-2-9 19:26
老师，是我换算错了。
请问另外一个同样用这种方法计算的金属配合物，的S0和T1的能量差是1.0eV，对应的波长处于近红外，是否不合理，而T2是在2.0eV，这时，我该选取S0和T2的能量差来做为三重态能级么？求老师指教！

作者
Author: sobereva 时间: 2015-2-9 19:38
可以用M062X算算试试，此方法三重态激发能计算误差较小。没法光凭这个就用T2作为三重态能级。

作者
Author: 哈哈宝 时间: 2015-2-9 19:52
谢谢老师，我再试一下，可以基态用B3，然后激发态用M062X么？

作者
Author: sobereva 时间: 2015-2-9 21:57
不知道你指的基态用B3LYP是什么意思。
如果你是用B3LYP做T1优化，然后再用M062X算激发能，没什么问题。但最好还是一致。

作者
Author: 哈哈宝 时间: 2015-2-10 07:54
明白了，谢谢老师不吝赐教！

作者
Author: blueyangliu 时间: 2015-2-10 10:08
本帖最后由 blueyangliu 于 2015-2-10 10:11 编辑

不知道你用什么软件算的。高斯的话，输出中激发能eV和波长，还有三重态能量都是直接给出的。激发能1.86eV对应于667nm.

作者
Author: 哈哈宝 时间: 2015-2-11 08:35
谢谢blueyangliu，我就是用高斯做的，用sob老师的方法，我用M062X基组做了一下激发能。
只有一个化合物的三重态激发能在 2.3845 eV 519.95 nm，与实验值的波长相差一些，但是能接受。可是另外两个化合物的激发能却在 1.3158 eV 942.24 nm，而 Excited State 2: Triplet-A 2.3129 eV 536.06 nm ，似乎更符合实验上的数据，我又重新优化激发态几何结构了，不知道是不是优化得到的不是能量最低点。正在尝试。

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