在单三态能量差ΔEST以及磷光的计算中都会涉及第一激发三重态T1的几何结构优化。本文单独谈一下T1的结构优化。T1的结构优化常用以下两种方式:(1)将自旋多重度设为3,使用UDFT(非限制性DFT)来优化;(2)使用TD-DFT来优化,使用TD(triplet)关键词。以下我们通过两个例子来对比两种方法的优化结果。
一、联苯分子
优化基态的输入文件:
代码语言:javascript复制#p pbe1pbe/6-311G** opt
Title Card Required
0 1
C -2.09908937 0.62433012 0.00000000
C -0.69667737 0.62433012 0.00000000
C 0.00686463 1.82694412 0.00000000
C -0.67921737 3.04124012 -0.00276100
C -2.07387037 3.04933912 -0.00410600
C -2.78183137 1.84928112 -0.00129900
C -2.83952037 -0.63611388 0.00125500
C -2.40815137 -1.71412788 -0.78523500
C -3.11629937 -2.91403288 -0.78348600
C -4.25927837 -3.05309588 0.00359000
C -4.69372537 -1.98593688 0.78943600
C -3.99048837 -0.78310388 0.78885600
H -0.15162637 -0.33154288 0.00768300
H 1.10672063 1.81686712 0.00327200
H -0.12224837 3.98928512 -0.00385900
H -2.61816537 4.00512212 -0.00849600
H -3.88213737 1.86004112 -0.00794900
H -1.51013437 -1.60658188 -1.41200100
H -2.77215737 -3.75282488 -1.40621400
H -4.81617637 -4.00118288 0.00450100
H -5.59347037 -2.09282188 1.41300000
H -4.33345437 0.05459012 1.41453700
UDFT的输入文件:
代码语言:javascript复制#p pbe1pbe/6-311G** opt
Title Card Required
0 3
坐标略
TD-DFT的输入文件:
代码语言:javascript复制#p pbe1pbe/6-311G** opt td(triplet)
Title Card Required
0 1
坐标略
以下是优化得到的结构:
优化得到的基态结构中,两个苯环的夹角为39.9度,而TD-DFT和UDFT优化得到的夹角均为0度,两者一致。由于基态和激发态结构的差异,也使得联苯分子的光谱中能观察到较大的Stokes位移。
二、CPDT分子
结构如下图所示,来自文献Adv. Funct. Mater. 2020, 1910301,计算时将侧链简化为甲基。
UDFT的输入文件为
代码语言:javascript复制#p opt scrf=(solvent=chloroform) upbe1pbe/def2svp
[No Title]
0 3
C -4.27441500 -2.27361200 0.21562800
C -4.10289800 -0.95436000 0.10922600
C -2.79048500 -0.70090600 0.12183500
C -2.14232600 -1.86237600 0.22902900
C -3.01121800 -3.09411600 0.30036000
S -1.88629700 0.44048500 0.04446400
C -0.61228700 -0.31622600 0.12773500
C -0.82805200 -1.64431600 0.23486000
C -5.57527700 -2.56104800 0.19944200
C -6.26627400 -1.40498800 0.14989700
S -5.36607300 -0.23337600 -0.00310900
C 2.12610300 1.83386600 -0.03217000
C 2.77201200 0.67073300 0.06512100
C 4.08528200 0.91730200 0.03911000
C 4.26089600 2.23576700 -0.07562200
C 2.99868600 3.06120500 -0.13403100
S 5.34096700 0.17939700 0.10500000
C 6.23816000 1.34527700 -0.07352500
C 5.56378200 2.51127700 -0.11882400
C 0.81167800 1.61697000 -0.01821800
C 0.59615700 0.28877600 0.09084500
S 1.86967900 -0.46986400 0.16691000
C -11.61443982 -3.95076415 0.47232452
C -11.74659167 -5.25155927 0.78659067
C -10.67005589 -6.03358835 0.96782128
C -9.43913809 -5.51529723 0.83879395
C -9.33877111 -4.21464541 0.53173444
C -10.39108086 -3.40946319 0.34457975
C -8.19255853 -3.51959142 0.36820696
C -9.93097513 -2.18443152 0.00367941
O -7.07639588 -3.97938324 0.48541112
C -10.71255268 -1.18734349 -0.46568022
C 11.44754589 4.08455644 -0.72201735
C 11.51138846 5.38375359 -1.06290925
C 10.40002021 6.07974045 -1.35378676
C 9.20086856 5.47928980 -1.30274563
C 9.16555626 4.18624786 -0.95285511
C 10.25592151 3.46792143 -0.65635830
C 8.05932177 3.41897423 -0.83957488
C 9.86537643 2.19872795 -0.40735157
O 6.91948431 3.78555373 -1.02920097
C 10.70233435 1.14170220 -0.31668287
C 7.58408600 1.28441400 -0.05746600
C -7.60432000 -1.30125400 0.02166100
C 8.50757700 2.20642103 -0.40364757
C -8.57765221 -2.23207623 0.12233890
C -2.82666200 -4.01372800 -0.90841400
C -2.88360000 -3.82588100 1.63827600
C 2.84919900 3.98827900 1.07396800
C -12.02833839 -1.23426811 -0.56497771
N -13.18408882 -1.23188865 -0.66374030
C -10.25704015 -0.09407965 -1.04879248
N -9.88727401 0.86568539 -1.58492031
C 12.02001076 1.20810835 -0.33844275
N 13.17924977 1.23077795 -0.36457461
C 10.28528364 -0.10565884 -0.45628795
N 9.93893211 -1.19190004 -0.67409753
C 2.83540800 3.78478200 -1.47241500
H -0.08610600 -2.45000000 0.30979100
H -5.96666200 -3.58571600 0.20814200
H 5.98625000 3.52068900 -0.20939500
H 0.06936200 2.42246500 -0.09140100
H -12.53928772 -3.37659812 0.33764705
H -12.75453667 -5.68781782 0.89485477
H -10.79637682 -7.09915697 1.22255937
H -8.55153820 -6.15074752 0.98723032
H 12.39436279 3.57509994 -0.50283235
H 12.49196712 5.88784718 -1.11036999
H 10.47196716 7.14376164 -1.63494916
H 8.28530791 6.04433016 -1.53953759
H 7.96068829 0.27756903 0.20434873
H -7.94074014 -0.30174767 -0.31832409
H -1.81845500 -4.48644200 -0.89486400
H -3.58769600 -4.82673000 -0.90641300
H -2.92786100 -3.45445400 -1.86637900
H -1.88002500 -4.29827700 1.73660700
H -3.65199200 -4.62733600 1.72530800
H -3.01834300 -3.13122300 2.49850200
H 3.61331300 4.79793800 1.04752000
H 1.84289000 4.46507000 1.08426900
H 2.97335600 3.43398800 2.03213600
H 3.59848200 4.58824200 -1.58317600
H 2.95047200 3.08515500 -2.33156000
H 1.82835100 4.25409100 -1.54790400
TD-DFT的输入文件为
代码语言:javascript复制#p opt td=triplet scrf=(solvent=chloroform) pbe1pbe/def2svp
[No Title]
0 1
结构略
两者优化得到的结构为
RMSD值为0.01 Å,差别非常小。两种方法优化都是37圈收敛,而UDFT使用的时间为3 h 30 min,TD-DFT所用时间为12 h 35 min,差别还是非常大的。计算在两颗Intel Platinum 8163 CPU(共48核)上完成。
从以上两例可以看出,UDFT和TD-DFT优化得到的结构基本是一致的,但是TD-DFT的优化时间要比UDFT长很多。一般计算中使用UDFT优化即可。