计算化学公社

标题: 咨询一下溶剂化效应下的荧光问题 [打印本页]

作者
Author: 落花飘雨679 时间: 2016-11-2 15:20
标题: 咨询一下溶剂化效应下的荧光问题
老师，我用PBE0泛函算溶剂化下的荧光，算激发td时，发现激发到了第一激发态，然后对root=1 td opt,发现也是在第一激发态发射，但是实验给出的是两个峰，不知道第二个峰所对应的激发态应该怎样找到，希望老师指点，谢谢老师了。下面是我的td和td-opt的结果：
td:
Excited State 1:    Singlet-A    2.7899 eV  444.41 nm  f=0.4083  <S**2>=0.000
   152 ->156       0.13160
   154 ->155       0.69045
This state for optimization and/or second-order correction.
Total Energy, E(TD-HF/TD-KS) =  -2682.03111434
Copying the excited state density for this state as the 1-particle RhoCI density.

Excited State 2:    Singlet-A    2.9626 eV  418.50 nm  f=0.0000  <S**2>=0.000
   152 ->155       0.59468
   154 ->156       -0.37035

Excited State 3:    Singlet-A    3.0183 eV  410.77 nm  f=0.0000  <S**2>=0.000
   152 ->157       -0.12354
   153 ->155       0.68909
td-opt:
Excited State 1:    Singlet-A    2.3783 eV  521.32 nm  f=0.4698  <S**2>=0.000
   152 ->156       0.11429
   154 ->155       0.69572
This state for optimization and/or second-order correction.
Total Energy, E(TD-HF/TD-KS) =  -2682.03873392
Copying the excited state density for this state as the 1-particle RhoCI density.

Excited State 2:    Singlet-A    2.5030 eV  495.34 nm  f=0.0000  <S**2>=0.000
   152 ->155       0.70335

Excited State 3:    Singlet-A    2.6246 eV  472.39 nm  f=0.0000  <S**2>=0.000
   152 ->157       0.10341
   153 ->155       0.69507

Excited State 4:    Singlet-A    2.8613 eV  433.31 nm  f=0.0000  <S**2>=0.000
   154 ->156       0.69395

Excited State 5:    Singlet-A    3.0743 eV  403.30 nm  f=0.4370  <S**2>=0.000
   153 ->156       0.67290
   154 ->157       0.16274

作者
Author: 小范范1989 时间: 2016-11-2 16:00
也就是实验室双荧光，你这里计算的只有一个是吧？
1：我看你计算的文件，td和opt td。激发态之间的能量间隔不是很大，荧光应该是从S1发出的，高激发态发光可能性比较小。
2：所以，你既然找到了一个和实验符合的荧光，另外一个峰我感觉也是从s1发出来的。那么就是另外一个构型。
我感觉你可以这样试试：比如说针对某一个角度进行一下柔性扫描。可能会发现两个或者几个结构，在发现的结构在opt td。
说白了：也就是分子可能有其他构型发光。你只是找到了一个而已。
个人观点，不一定对。望指点。

作者
Author: 落花飘雨679 时间: 2016-11-2 16:19

小范范1989 发表于 2016-11-2 16:00
也就是实验室双荧光，你这里计算的只有一个是吧？
1：我看你计算的文件，td和opt td。激发态之间的能量间 ...

对，实验是双荧光，请问一下：“针对某个方向进行柔性扫描”应该怎样做，刚刚开始接触，很多地方都不懂，还望前辈指教，感激不尽

作者
Author: 小范范1989 时间: 2016-11-2 17:16

落花飘雨679 发表于 2016-11-2 16:19
对，实验是双荧光，请问一下：“针对某个方向进行柔性扫描”应该怎样做，刚刚开始接触，很多地方都不懂， ...

比如说D-A之间的夹角了等等，我是这个意思

作者
Author: 落花飘雨679 时间: 2016-11-2 18:19

小范范1989 发表于 2016-11-2 17:16
比如说D-A之间的夹角了等等，我是这个意思

我查了一下，就是说我要选择一个二面角进行扫描，但是我算的这个分子对称性很高，我应该选择哪种类型的二面角呢

作者
Author: bomsaude 时间: 2016-11-2 18:34

落花飘雨679 发表于 2016-11-2 10:19
我查了一下，就是说我要选择一个二面角进行扫描，但是我算的这个分子对称性很高，我应该选择哪种类型的二 ...

对于单个的有机分子，比如几十个原子的分子，Spartan构像扫描还是挺好的，软件可参考yuanhe同学的帖子（http://ccc.keinsci.com/forum.php ... p;highlight=spartan）

作者
Author: 落花飘雨679 时间: 2016-11-2 18:43

bomsaude 发表于 2016-11-2 18:34
对于单个的有机分子，比如几十个原子的分子，Spartan构像扫描还是挺好的，软件可参考yuanhe同学的帖子（h ...

好的，谢谢前辈指点

作者
Author: 落花飘雨679 时间: 2016-11-2 18:57

小范范1989 发表于 2016-11-2 17:16
比如说D-A之间的夹角了等等，我是这个意思

请问如果做柔性扫描，坐标一定要用内坐标形式吗

作者
Author: 冰释之川 时间: 2016-11-4 08:39

落花飘雨679 发表于 2016-11-2 18:57
请问如果做柔性扫描，坐标一定要用内坐标形式吗

不一定的

作者
Author: 落花飘雨679 时间: 2016-11-4 11:56

冰释之川发表于 2016-11-4 08:39
不一定的

前辈能不能提供一下有关这方面的资料，真的不知道该怎么办了，谢谢前辈了

作者
Author: 冰释之川 时间: 2016-11-4 12:36

落花飘雨679 发表于 2016-11-4 11:56
前辈能不能提供一下有关这方面的资料，真的不知道该怎么办了，谢谢前辈了

举个例子：

%mem=2000MB
%nprocshared=4
#p cam-b3lyp/6-31g(d) opt=modredundant

SO_R

0 1
C                -4.50397196 -0.75194846 0.10260019
C                -4.00537288 -1.53817228 1.13526414
C                -2.68971631 -1.36896666 1.57307320
C                -1.89683235 -0.41598133 0.96407877
C                -2.40908537 0.36750331 -0.07116942
C                -3.71236081 0.21622226 -0.51562741
C                -0.47689694 0.01851284 1.22005191
C                -0.14738800 0.84877200 -0.04661900
N                -1.44071808 1.29042861 -0.50475390
C                0.79254000 2.00468300 0.21586900
N                1.98914887 1.88824934 0.62580199
C                2.60570900 0.62031100 0.15047400
C                1.84835200 -0.40515300 -0.42404700
C                2.41973500 -1.63950300 -0.70144300
C                3.75255000 -1.85895800 -0.37498400
C                4.51355600 -0.85373800 0.21841600
C                3.94024900 0.38516200 0.46878200
O                0.54607000 -0.18793200 -0.76711200
H                -5.52683865 -0.89210854 -0.23369007
H                -4.63820976 -2.28630309 1.60111859
H                -2.29563737 -1.98387029 2.37726960
H                -4.10657288 0.82784778 -1.32128637
H                0.37678242 3.00052272 0.06071824
H                1.80915600 -2.40934300 -1.16075000
H                4.19868800 -2.82670700 -0.58197700
H                -0.40994231 0.65398290 2.11174093
H                0.24174725 -0.79638158 1.33462941
H                -1.46359610 1.56618881 -1.47944368
H                4.50965500 1.19888000 0.90555200
H                5.55266200 -1.03390400 0.47259400

* 8 10 * R
9 8 10 11 S 100 1.0 // 这是二面角扫描，扫描100步，步长为1°

作者
Author: 落花飘雨679 时间: 2016-11-4 13:41

冰释之川发表于 2016-11-4 12:36
举个例子：

%mem=2000MB

请问老师，下面这个* 8 10 * R代表什么意思

作者
Author: 冰释之川 时间: 2016-11-4 14:05

落花飘雨679 发表于 2016-11-4 13:41
请问老师，下面这个* 8 10 * R代表什么意思

这一行好像可以不用写，你去掉试试看，很久以前的任务了……

作者
Author: 落花飘雨679 时间: 2016-11-4 14:12

冰释之川发表于 2016-11-4 14:05
这一行好像可以不用写，你去掉试试看，很久以前的任务了……

好的，谢谢老师

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