Gaussian中用TDDFT计算激发态和吸收、荧光、磷光光谱的方法

xiaoming

sobereva 发表于 2016-1-13 18:14
我也读不懂你的意思。non-pcm指什么，步骤具体怎么弄的都没交代清楚。
是否考虑非平衡溶剂有可能结果差 ...

老师，您好，，是这样的，我之前 是用溶液中基态优化的几何结构，做一个TD计算，然后读Excited State   1:  2.98ev， 认为这个能量就是吸收光谱，Total Energy, E(TD-HF/TD-KS) 就是改激发态下的总能量。   我现在看了帖子，，认为我之前这样做是没有考虑非平衡溶剂化的。。   按照文中的3步，我做了吸收， 得到After PCM corrections, the energy is  -649.440045444     a.u.，，减去第一步中的总能量，得到3.0eV是考虑了非平衡溶剂化后的吸收光谱能量。  这对吧？  但是我想问第3步计算 最后 给出来的  Excited State   1:    3.2390 eV   Total Energy, E(TD-HF/TD-KS) =  -649.432379141 ，，这个能量是什么？ 为什么不认为这个是垂直激发的能量？   我用这种3步计算的吸收光谱的方法， 输入文件没有density，  但是最后也输出了CI 密度分布  和 偶极矩，这是激发态的电荷分布和偶极矩么？

sobereva

xiaoming 发表于 2016-1-13 18:55
老师，您好，，是这样的，我之前 是用溶液中基态优化的几何结构，做一个TD计算，然后读Excited State   1 ...

你说的第3步是哪个第三步？最好说明对应文中那一节哪一步，并且关键词明确列出，要不然容易混淆。

xiaoming

sobereva 发表于 2016-1-14 02:41
你说的第3步是哪个第三步？最好说明对应文中那一节哪一步，并且关键词明确列出，要不然容易混淆。

2.5节中提到的高斯7步中的前3步，算溶剂环境下的吸收光谱

sobereva

xiaoming 发表于 2016-1-14 15:51
2.5节中提到的高斯7步中的前3步，算溶剂环境下的吸收光谱

所以说要忘记所谓的“七步”

不写density得到的是基态的密度、基态的原子电荷、基态的偶极矩

帖子里没提到的量都不用管
E(TD-HF/TD-KS) =  -649.432379141 至多只能说是激发态能量，显然不是激发能。这在非平衡溶剂计算时什么也不是

mengyingbian

sob老师，算T1 T2 T3时输出文件中哪个值是T1 T2 T3的值呀

healthyworld

太牛了

杨小狗

sobereva 发表于 2015-12-23 07:18
TDDFT优化T1，最后一次输出的这个量就是磷光发射能。

不知道别人问没问过，我在这里有两个问题等待老师有空解答。
1.请问您：如何判断一个分子有没有三重态？是比如氧原子P轨道数那种的分析吗？
2.并行提交gaussian作业时，Keywords加上out=wfn 最后一空行下边应该怎么写输出路径呢？
感谢您在百忙之中的阅读。

sobereva

mengyingbian 发表于 2016-5-10 21:45
sob老师，算T1 T2 T3时输出文件中哪个值是T1 T2 T3的值呀

基态若是单重态，你写TD(triplet)，输出的最低的三个态就是T1、T2、T3

sobereva

杨小狗 发表于 2016-5-11 20:15
不知道别人问没问过，我在这里有两个问题等待老师有空解答。
1.请问您：如何判断一个分子有没有三重态？ ...

1 任何一个分子都有三重态
2 这和并行提交没有关系啊，不管你怎么提交，输入文件该怎么写还是怎么写。至于具体路径，根据系统实际情况来写，不确定的话可以找个很小体系试试。

ter20

Sob老师，关于T1的计算有几个问题想请教
1. 如您所说，除了通过UDFT 基态自旋多重度为3优化得到T1，还可以用TDDFT优化T1
关键字：# B3LYP/6-31G* TD(triplet) Opt
自旋多重度： 0 1
a) 这里为何要将自旋多重度设为1？是因为体系的基态为singlet吗？
b) 在计算Quintet等高自旋态时，这里是不是需要测试选取不同的自旋多重度?即选取不同的参考态

2. 在计算自由基体系，例如单自由基（doublet）时，TD关键字不需要指定Singlet、Triplet，直接TD(Root=n,Nstates=n)就好了吧？（这个问题很Silly...）

3. 您提到的根据该发射能，通过Multiwfn作出的荧光/磷光发射谱其实就是直接对一个激发能数据展宽来绘制谱图对吧？我的理解是必须通过TD计算得到的log文件，将除了S1之外的振子强度设为0然后才好绘图对吧？如果只有发射能数据，除非手动做一个log文件，没法直接用Multiwfn绘制发射谱对吧

sobereva

ter20 发表于 2016-6-28 14:14
Sob老师，关于T1的计算有几个问题想请教
1. 如您所说，除了通过UDFT 基态自旋多重度为3优化得到T1，还可以 ...

1 基态是几重态就设几。TD里面设的是基于这个基态波函数要算的激发态
2 对。而且基态是开壳层的时候，本身也没法在TD里设定要算几重态，只能是从算出来的一批态里面找自己要的。
3 对。但是Multiwfn也可以从普通文本文件里读取信息来作图，不管你用什么程序做的计算，只要把激发能、振子强度都按照指定格式写成一个文本文件，Multiwfn就能从里面读取数据来作图。格式要求见手册3.13.2节靠后部分。

ter20

sobereva 发表于 2016-6-29 03:05
1 基态是几重态就设几。TD里面设的是基于这个基态波函数要算的激发态
2 对。而且基态是开壳层的时候， ...

谢谢Sob老师，受益匪浅！

ter20

还有两个问题请教Sob老师，求指教
1. 如果基态是Singlet，在做激发态优化时，我用了# B3LYP/6-31G* TD(50-50)，那优化结果是Singlet还是Triplet怎么判断？
2. 在TD计算的输出部分会有相关的轨道跃迁贡献，这些轨道是对应基态的轨道还是激发态的轨道？换句话说，如果我们绘制轨道跃迁贡献图，需要在做TD计算时加上Density关键字，然后用相应的chk文件里的信息来绘制轨道吗？还是不加Density直接绘制轨道？

ter20

又打扰Sob老师了，另外两个问题是：
1. 如果做TD优化，# B3LYP/6-31G* TD(root=2) Opt，得到了第二激发态的结构，那么该结果文件中的
Excited State   1:      Singlet-B2     5.3446 eV  231.98 nm  f=0.1525  <S**2>=0.000
      19 -> 23         0.18851
      20 -> 22         0.68003

Excited State   2:      Singlet-A1     6.6881 eV  185.38 nm  f=0.0000  <S**2>=0.000
      19 -> 22         0.55398
      20 -> 23        -0.42789

Excited State  3:      Singlet-A2     6.8138 eV  181.96 nm  f=0.0000  <S**2>=0.000
      18 -> 22         0.70308
这样的信息，我的理解这就是在第二激发态Minimum结构下的垂直激发计算结果，这和直接用该结构做的TD单点计算（# B3LYP/6-31G* TD(root=1)）应该是一样的对吧？

2. 如果激发态可能是三重态（并不确定），做TD优化时，是不是最好分开做，先做Singlet，再做Triplet？不然感觉混在一起比较头大

sobereva

ter20 发表于 2016-6-29 11:50
还有两个问题请教Sob老师，求指教
1. 如果基态是Singlet，在做激发态优化时，我用了# B3LYP/6-31G* TD(50- ...

1 切勿这么写，否则自己都不知道到底优化的哪个。而且只需要优化一种自旋的激发态，同时算出两种自旋，也多费时间。
2 TD没有所谓的激发态轨道，轨道就是基态的。
不知道你说的轨道跃迁贡献图具体指什么，反正用不着density关键词。