计算化学公社

 找回密码 Forget password
 注册 Register
Views: 16854|回复 Reply: 17
打印 Print 上一主题 Last thread 下一主题 Next thread

[新手求助] 荧光的计算问题

[复制链接 Copy URL]

67

帖子

0

威望

265

eV
积分
332

Level 3 能力者

本帖最后由 落花飘雨679 于 2017-3-11 11:30 编辑

请教各位老师,我在计算荧光时,做激发td时,激发到了S2,然后在此基础是上做了发射td (root=2) opt ,发现从S1回到基态,但是此s1的能量高与s2的能量,请问是不是说明出现了势能面的交叉?然后我在此基础上继续做了td(root=1)  opt ,结果为s2的振子强度最大,请问老师这是不是能够说明此荧光 其实是从s2发射的?这是我计算的分子,优化完没有出虚频后发现它的振动方式主要是红色的硅上烷基的振动,不知道该怎么选择角度进行扫描,以前也进行过的方面的尝试,但是太耗时了,所以就不知道该如何找到另一个极小点的结构,不能与实验上的双峰相对应,望老师们给出意见,感谢.

评分 Rate

参与人数
Participants 1
eV +2 收起 理由
Reason
kechengbiao + 2 欢迎讨论

查看全部评分 View all ratings

1295

帖子

0

威望

6289

eV
积分
7584

Level 6 (一方通行)

2#
发表于 Post on 2017-3-9 20:33:19 | 只看该作者 Only view this author
1:是不是这个情况 http://bbs.keinsci.com/forum.php?mod=viewthread&tid=4931
2:或者是不是分子有两个构型?
https://www.x-mol.com/groups/fan_jianzhong

strive for greatness

67

帖子

0

威望

265

eV
积分
332

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2017-3-9 21:13:18 | 只看该作者 Only view this author
小范范1989 发表于 2017-3-9 20:33
1:是不是这个情况 http://bbs.keinsci.com/forum.php?mod=viewthread&tid=4931
2:或者是不是分子有两个 ...

对,就是这种情况,我觉着分子应该是有两个构型的,因为实验中给出的是双峰,而我计算的结果是只有一个峰,请问我应该怎样找着两个构型?

67

帖子

0

威望

265

eV
积分
332

Level 3 能力者

4#
 楼主 Author| 发表于 Post on 2017-3-9 21:15:45 | 只看该作者 Only view this author
落花飘雨679 发表于 2017-3-9 21:13
对,就是这种情况,我觉着分子应该是有两个构型的,因为实验中给出的是双峰,而我计算的结果是只有一个峰 ...

请教老师,对于这种交叉的情况应该这样处理,认为是s1发射的还是s2发射的呢?感谢您的回复

1295

帖子

0

威望

6289

eV
积分
7584

Level 6 (一方通行)

5#
发表于 Post on 2017-3-9 21:17:18 | 只看该作者 Only view this author
落花飘雨679 发表于 2017-3-9 21:13
对,就是这种情况,我觉着分子应该是有两个构型的,因为实验中给出的是双峰,而我计算的结果是只有一个峰 ...


我的感觉哈,不一定对
1:看看哪个角度容易变化,做一个柔性扫描。如果得到两个最低点,然后在做优化。
2:sob老师的那个构象搜索:http://bbs.keinsci.com/forum.php?mod=viewthread&tid=577
个人见解,不一定准确哈
https://www.x-mol.com/groups/fan_jianzhong

strive for greatness

67

帖子

0

威望

265

eV
积分
332

Level 3 能力者

6#
 楼主 Author| 发表于 Post on 2017-3-9 21:25:02 | 只看该作者 Only view this author
小范范1989 发表于 2017-3-9 21:17
我的感觉哈,不一定对
1:看看哪个角度容易变化,做一个柔性扫描。如果得到两个最低点,然后在做优化。
...

好的,谢谢您,我先看看学习一下

909

帖子

37

威望

5527

eV
积分
7176

Level 6 (一方通行)

7#
发表于 Post on 2017-3-9 22:58:52 | 只看该作者 Only view this author
落花飘雨679 发表于 2017-3-9 21:25
好的,谢谢您,我先看看学习一下

你这个不要扫描了 刚性结构。
另外你描述的不清楚,什么叫做“做激发td时,激发到了S2”
什么叫做“发现从S1回到基态”
最好是把log文件贴上来才好分析。

67

帖子

0

威望

265

eV
积分
332

Level 3 能力者

8#
 楼主 Author| 发表于 Post on 2017-3-10 08:12:16 | 只看该作者 Only view this author
ggdh 发表于 2017-3-9 22:58
你这个不要扫描了 刚性结构。
另外你描述的不清楚,什么叫做“做激发td时,激发到了S2”
什么叫做“发 ...

这是做激发时的.log文件,可以看到主要激发到了s2,
Excited State   1:      Singlet-A      3.4364 eV  360.80 nm  f=0.0068  <S**2>=0.000
     176 ->179         0.55860
     177 ->181         0.14190
     178 ->180         0.34449
     178 ->181        -0.12843
This state for optimization and/or second-order correction.
Total Energy, E(TD-HF/TD-KS) =  -2919.29994700   
Copying the excited state density for this state as the 1-particle RhoCI density.

Excited State   2:      Singlet-A      3.4422 eV  360.19 nm  f=0.8481  <S**2>=0.000
     175 ->179         0.26560
     177 ->179         0.17188
     178 ->179         0.60975

Excited State   3:      Singlet-A      3.6544 eV  339.27 nm  f=0.3497  <S**2>=0.000
     176 ->179        -0.34442
     177 ->180        -0.16779
     177 ->181         0.14359
     178 ->180         0.42326
     178 ->181        -0.36210
然后在此结构下做发射,优化的root=2,得到的结果文件为:
Excited State   1:      Singlet-A      2.8177 eV  440.02 nm  f=1.0015  <S**2>=0.000
     175 ->179        -0.24642
     176 ->179        -0.13546
     178 ->179        -0.63734

Excited State   2:      Singlet-A      2.4442 eV  507.25 nm  f=0.0000  <S**2>=0.000
     177 ->179        -0.69504
This state for optimization and/or second-order correction.
Total Energy, E(TD-HF/TD-KS) =  -2919.31666766   
Copying the excited state density for this state as the 1-particle RhoCI density.

Excited State   3:      Singlet-A      3.3739 eV  367.48 nm  f=2.2248  <S**2>=0.000
     174 ->181        -0.19307
     175 ->179        -0.12502
     177 ->180         0.64167
可以看到s1的振子强度最大,但是s1的能量高于s2,应该是出现了势能面的交叉吧,然后我又在此结构下优化了root=1,得到的结果为:
Excited State   1:      Singlet-A      2.4440 eV  507.29 nm  f=0.0000  <S**2>=0.000
     177 ->179         0.69504
This state for optimization and/or second-order correction.
Total Energy, E(TD-HF/TD-KS) =  -2919.31666767   
Copying the excited state density for this state as the 1-particle RhoCI density.

Excited State   2:      Singlet-A      2.8175 eV  440.05 nm  f=1.0015  <S**2>=0.000
     175 ->179         0.24641
     176 ->179         0.13549
     178 ->179         0.63734

Excited State   3:      Singlet-A      3.3739 eV  367.48 nm  f=2.2248  <S**2>=0.000
     174 ->181        -0.19308
     175 ->179        -0.12498
     177 ->180         0.64177
发现又变成了从s2发射,老师这种现象我应该认为是从s2发射的吗,而且实验中给出的双峰我应该通过什么方法找到?

6万

帖子

99

威望

6万

eV
积分
125140

管理员

公社社长

9#
发表于 Post on 2017-3-10 11:00:01 | 只看该作者 Only view this author
落花飘雨679 发表于 2017-3-10 08:12
这是做激发时的.log文件,可以看到主要激发到了s2,
Excited State   1:      Singlet-A      3.4364 eV ...


这些信息只是给出各激发态的特征,不叫“激发到了”某个态,除非是在现实中只给此分子照射相应频率的光。提问时只需要对计算得到的数据进行说明,而不要加入额外的物理意义上的猜想以免让别人混淆。

这种牵扯到交叉的问题,如果搞不清楚情况,就把步长上限固定在某个值(不固定往往也没事),然后把每一步的基态、S1、S2能量全提出来,绘制成“优化步数-能量”的曲线图,一看就知道怎么回事。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

909

帖子

37

威望

5527

eV
积分
7176

Level 6 (一方通行)

10#
发表于 Post on 2017-3-10 13:42:08 | 只看该作者 Only view this author
从结果看比较明显的是HOMO到LUMO是S2亮态,HOMO-1到LUMO是S1暗态,不过你的两次优化,其实优化的都是那个暗态。我建议你拿出最后一次优化的结构,优化S2(也就是亮态)。看看是个什么结果。另外实验上来看的话 这个分子荧光强不强?

67

帖子

0

威望

265

eV
积分
332

Level 3 能力者

11#
 楼主 Author| 发表于 Post on 2017-3-10 14:11:02 | 只看该作者 Only view this author
ggdh 发表于 2017-3-10 13:42
从结果看比较明显的是HOMO到LUMO是S2亮态,HOMO-1到LUMO是S1暗态,不过你的两次优化,其实优化的都是那个暗 ...

嗯嗯,好的,我优化一下试试,实验上的两个峰都挺强的

909

帖子

37

威望

5527

eV
积分
7176

Level 6 (一方通行)

12#
发表于 Post on 2017-3-10 14:14:55 | 只看该作者 Only view this author
落花飘雨679 发表于 2017-3-10 14:11
嗯嗯,好的,我优化一下试试,实验上的两个峰都挺强的

两个峰都强是什么意思?难道是双荧光?

67

帖子

0

威望

265

eV
积分
332

Level 3 能力者

13#
 楼主 Author| 发表于 Post on 2017-3-10 14:15:39 | 只看该作者 Only view this author
sobereva 发表于 2017-3-10 11:00
这些信息只是给出各激发态的特征,不叫“激发到了”某个态,除非是在现实中只给此分子照射相应频率的光 ...

好,谢谢老师回复

67

帖子

0

威望

265

eV
积分
332

Level 3 能力者

14#
 楼主 Author| 发表于 Post on 2017-3-10 14:17:20 | 只看该作者 Only view this author
ggdh 发表于 2017-3-10 14:14
两个峰都强是什么意思?难道是双荧光?

老师,这个分子在实验中就是双荧光,计算时只能计算出一个,不知道另外一个怎么找到

909

帖子

37

威望

5527

eV
积分
7176

Level 6 (一方通行)

15#
发表于 Post on 2017-3-10 14:19:38 | 只看该作者 Only view this author
落花飘雨679 发表于 2017-3-10 14:17
老师,这个分子在实验中就是双荧光,计算时只能计算出一个,不知道另外一个怎么找到

有意思,两个波长分别是多少,什么溶剂。考虑过excimer没有?

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-20 21:47 , Processed in 0.168739 second(s), 21 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list