CyRhmU.jpeg
南方科技大学公共卫生及应急管理学院2025级博士研究生招生报考通知
查看: 2367  |  回复: 7

NPB

金虫 (小有名气)

[求助] 如何确定单三重态能极差已有1人参与

优化S1,T1结构后,下一步如何计算单三重态能极差?如图0.24怎么计算出来的
如何确定单三重态能极差
回复此楼
锲而不舍,金石可镂
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lastzealot

新虫 (著名写手)

【乱世的奸雄】

做两次结构优化,一次是单重态,一次是三重态,能量做差

发自小木虫Android客户端
2楼2016-02-28 12:39:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

NPB

金虫 (小有名气)

引用回帖:
2楼: Originally posted by lastzealot at 2016-02-28 12:39:30
做两次结构优化,一次是单重态,一次是三重态,能量做差

再问一下,这个能量是输出文件最下面一个能量吗
SCF Done:  E(UB3LYP) =  -1601.07440849     A.U. after    7 cycles
锲而不舍,金石可镂
3楼2016-02-28 17:04:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小范范1989

木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
NPB: 金币+20, ★★★★★最佳答案, 感谢你的回帖,刚开始做理论计算,不怎么懂,有机会再求教你,金币都给你,另送红花一枚,笑纳 2016-02-29 00:59:32
搜索: Excited State   1:  会出现如下的显示

Excited State   1:      Singlet-A      2.6047 eV  476.01 nm  f=0.5771  <S**2>=0.000
     231 -> 232       -0.68847
This state for optimization and/or second-order correction.
Total Energy, E(TD-HF/TD-KS) =  -2751.27385282   
Copying the excited state density for this state as the 1-particle RhoCI density.

Excited State   2:      Singlet-A      2.7436 eV  451.91 nm  f=0.2652  <S**2>=0.000
     230 -> 232       -0.68716

Excited State   3:      Singlet-A      2.8319 eV  437.81 nm  f=0.0096  <S**2>=0.000
     229 -> 232       -0.69327
SavETr:  write IOETrn=   770 NScale= 10 NData=  16 NLR=1 NState=    3 LETran=      64.


这个-2751.27385282    就是你要的s1的能量。

如果你的t1也是通过TDDFT的方式优化,那么找法类似。
如果你的t1的UDFT 也就是自选多重读为3的方式,就是文件最后的HF

ps:我猜测楼主是要做延时荧光 TADF。我是理论计算的,我做TADF。这个文章,我有印象,看过。应该是adachi教授的。

» 本帖已获得的红花(最新10朵)

It doesn't matter how slow you are, as long as you're determined to get there, you'll get there.
4楼2016-02-28 18:19:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

NPB

金虫 (小有名气)

引用回帖:
4楼: Originally posted by 小范范1989 at 2016-02-28 18:19:53
搜索: Excited State   1:  会出现如下的显示

Excited State   1:      Singlet-A      2.6047 eV  476.01 nm  f=0.5771  <S**2>=0.000
     231 -> 232       -0.68847
This state for optimiza ...

45|H,-2.2003306863,-3.415322583,-0.2622783912|H,-3.4433252834,-2.74634
87306,0.8039898971|H,-3.634496789,-2.6320231851,-0.9445539024|H,-3.670
4673771,2.5649888071,-0.8439327427|H,-3.4747285891,2.6325439253,0.9060
800456|H,-2.2432944211,3.3456926345,-0.1439553476||Version=IA32W-G09Re
vA.02|State=3-A|HF=-1601.0744085|S2=2.034985|S2-1=0.|S2A=2.000708|RMSD
=4.897e-009|RMSF=1.924e-006|Dipole=-0.3720291,0.0487153,-0.003622|Quad
rupole=10.3531146,-5.428082,-4.9250326,-1.0428136,1.3688025,-1.6679512
|PG=C01 [X(C38H31B1N2)]||@


EDUCATION WITHOUT COMMON SENSE IS A LOAD OF BOOKS ON THE BACK OF AN ASS.
Job cpu time:  0 days 12 hours 41 minutes 13.0 seconds.
File lengths (MBytes):  RWF=    652 Int=      0 D2E=      0 Chk=    312 Scr=      1
Normal termination of Gaussian 09 at*******

最后这堆数里面的这个数值呗,谢谢你的回复
HF=-1601.0744085
锲而不舍,金石可镂
5楼2016-02-29 01:01:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小范范1989

木虫 (著名写手)

引用回帖:
5楼: Originally posted by NPB at 2016-02-29 01:01:50
45|H,-2.2003306863,-3.415322583,-0.2622783912|H,-3.4433252834,-2.74634
87306,0.8039898971|H,-3.634496789,-2.6320231851,-0.9445539024|H,-3.670
4673771,2.5649888071,-0.8439327427|H,-3.4747285891,2 ...

-如果你的t1的UDFT 也就是自选多重读为3的方式,就是文件最后的HF。
如果是上面这个问题,那就是你说的这个HF
It doesn't matter how slow you are, as long as you're determined to get there, you'll get there.
6楼2016-02-29 08:05:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huihuibulaji

木虫 (正式写手)

送红花一朵
引用回帖:
4楼: Originally posted by 小范范1989 at 2016-02-28 18:19:53
搜索: Excited State   1:  会出现如下的显示

Excited State   1:      Singlet-A      2.6047 eV  476.01 nm  f=0.5771  <S**2>=0.000
     231 -> 232       -0.68847
This state for optimiza ...

我也是做理论计算的,我有个疑问,我的S1用了TD-DFT计算,读取的能量和你读取的一样,但是我个人觉得三重态的UB3优化后再做td(triplets)读取结果中的E(TD-HF/TD-KS)我不太明白这两种读取出来的差别在哪,或者说我这样读取正确吗?

发自小木虫Android客户端
7楼2016-02-29 10:49:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小范范1989

木虫 (著名写手)

引用回帖:
7楼: Originally posted by huihuibulaji at 2016-02-29 10:49:25
我也是做理论计算的,我有个疑问,我的S1用了TD-DFT计算,读取的能量和你读取的一样,但是我个人觉得三重态的UB3优化后再做td(triplets)读取结果中的E(TD-HF/TD-KS)我不太明白这两种读取出来的差别在哪,或者说 ...

UB3优化后再做td(triplets)读取结果中的E(TD-HF/TD-KS)我不太明白这两种读取出来的差别在哪,或者说我这样读取正确吗?这个问题我想过,是这样的:
1:优化T1值得一说。有两种做法,一种是用UDFT,即自旋多重度设3直接优化,这样是直接算T1,不牵扯TDDFT计算;另一种是自旋多重度设1同时用TD(triplet),这样每一步就是先计算单重态基态作为参考态,再通过TDDFT计算T1态并优化它。这两种方式原理上都对,结果通常相近,但强烈建议用UDFT方式,因为其耗时只比优化单重态基态增加一点,而用TDDFT优化T1的话计算量则要增加约一个数量级。

不过UDFT方法缺点是只能优化T1,而TDDFT则可以优化任意三重态激发态,而且能够同时给出轨道跃迁信息。
2FT和TDDFT的方法却别,有差距很正常。
It doesn't matter how slow you are, as long as you're determined to get there, you'll get there.
8楼2016-02-29 14:50:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 NPB 的主题更新
信息提示
请填处理意见