计算化学公社

标题: 关于 Ir(ppy)3金属络合物激发态的计算结果问题 [打印本页]

作者
Author:
Marc0    时间: 2023-8-16 15:01
标题: 关于 Ir(ppy)3金属络合物激发态的计算结果问题
使用高斯计算金属络合物 Ir(ppy)3 的激发态
输入文件如下(已预先完成结构优化)
%nprocshared=10
%mem=20GB
%chk=Irppy3_singlet_TD20.chk
#p td(50-50, nstates=10) b3lyp/genecp geom=check guess=read

Ir(ppy)SINGLET TD20

0 1

C H N
6-31G*
****
Ir
SDD
****

Ir
SDD

输出文件末尾显示"Normal termination...",并得到以下结果:

Excited State   1:      Triplet-A      2.6455 eV  468.67 nm  f=0.0000  <S**2>=2.000
     125 ->131        -0.10206
     126 ->133        -0.12261
     127 ->132        -0.12260
     128 ->132        -0.17893
     128 ->133         0.14106
     129 ->132        -0.14107
     129 ->133        -0.17897
     130 ->131         0.55852
This state for optimization and/or second-order correction.
Total Energy, E(TD-HF/TD-KS) =  -1540.60048920   
Copying the excited state density for this state as the 1-particle RhoCI density.

Excited State   2:      Triplet-E      2.6556 eV  466.88 nm  f=0.0000  <S**2>=2.000
     126 ->131        -0.12867
     128 ->131         0.20321
     128 ->132         0.13433
     128 ->133         0.10787
     129 ->131        -0.17730
     129 ->132         0.10788
     129 ->133        -0.13432
     130 ->133         0.51857

Excited State   3:      Triplet-E      2.6556 eV  466.88 nm  f=0.0000  <S**2>=2.000
     127 ->131        -0.12867
     128 ->131        -0.17730
     128 ->132        -0.10790
     128 ->133         0.13434
     129 ->131        -0.20321
     129 ->132         0.13431
     129 ->133         0.10786
     130 ->132         0.51857

Excited State   4:      Triplet-A      2.8973 eV  427.93 nm  f=0.0000  <S**2>=2.000
     128 ->132         0.45154
     129 ->133         0.45153
     130 ->131         0.19209

Excited State   5:      Singlet-A      2.9066 eV  426.55 nm  f=0.0059  <S**2>=0.000
     130 ->131         0.69620


我的问题是:
1. 输入文件(例如计算档次的选择、关键词设置)是否合理?
2. 计算结果中,Excited State 1和2 分别出现了 Triplet-A 和 Triplet-E,A 和 E 的含义是什么?
3. 发现 Excited State 1-3 的结果十分接近,我能否认为,这三个 Excited State 共同对应于 T1;而 Excited State 4 对应于 T2?

以上,希望各位老师不吝赐教!谢谢!




作者
Author:
wzkchem5    时间: 2023-8-16 15:48
1. 问题不大,但B3LYP未必是最适合的泛函,可能PBE0乃至某些范围分离泛函更适合。Ir光催化剂属于成熟的体系,有benchmark文章,可以查一下文献推荐什么泛函。
2. 不可约表示
3. 不是的,如果不加SOC的话,T1的每个简并态的能量只输出一次,而不是3次。所以你的计算结果表明最低的三个三重态的能量相近,每个三重态又包含3个microstate,一共9个microstate
作者
Author:
Marc0    时间: 2023-8-16 16:29
wzkchem5 发表于 2023-8-16 15:48
1. 问题不大,但B3LYP未必是最适合的泛函,可能PBE0乃至某些范围分离泛函更适合。Ir光催化剂属于成熟的体系 ...

谢谢wzk老师。
还有个小问题:Excited state 2 和 3 的能量数值完全一致该如何理解呢?有没有可能是它们之差太小了(比如一个是2.65557eV,一个是2.65564eV),四舍五入后都变成了2.6556eV?
此外,根据您的解释,我应该把Excited state 2 和 3分别看作 T2 和 T3,对吗?
作者
Author:
wzkchem5    时间: 2023-8-16 17:28
Marc0 发表于 2023-8-16 09:29
谢谢wzk老师。
还有个小问题:Excited state 2 和 3 的能量数值完全一致该如何理解呢?有没有可能是它们 ...

因为空间简并性。E表示是二维不可约表示,所以是二重简并的
对,这是一种方法,另一种方法是把excited state 2、3记为T2的两个简并分量,这两种方法都对,但是同一篇文章里采用的方法必须统一




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3