高斯优化激发态出错

小范范1989

yopoxia 发表于 2014-12-25 21:42
嗯  非常感谢您的热心回复。
不过还是有地方不明白，希望您能继续解答，谢谢！
1、如果我是要算荧光的 ...

对了，还想问你一个问题，你是高斯09，优化？用多核？

yopoxia

sobereva 发表于 2014-12-25 23:42
荧光算的是singlets，root=1。root>1不符合kasha规则，triplet则对应的是磷光
最后一轮激发态优化产生的 ...

嗯  非常安心sob老师的回复。
还有个疑问，希望sob cc能再解释一下。
有的时候最后一轮的Excited State  1所对应的 f 很小，接近于0了，如下举例Excited State  1:      Singlet-A      2.9736 eV  416.95 nm  f=0.0010  

Excited State   2:      Singlet-A      2.9992 eV  413.39 nm  f=0.1000  

Excited State   3:      Singlet-A      3.1772 eV  390.23 nm  f=0.0700  



Excited State   2的 f 相对Excited State 1大，这样的话，是不是可以这样理解，f 太小了，这个对应的只是计算上的荧光最大峰位置，而不是实验的。
实验上的应该是Excited State   2，这个相当于Excited State   1.

yopoxia

小范范1989 发表于 2014-12-26 10:46
对了，还想问你一个问题，你是高斯09，优化？用多核？

嗯 是的 用的是高斯09  
8核  16核都可以

去离子水

小范范1989 发表于 2014-12-26 10:46
对了，还想问你一个问题，你是高斯09，优化？用多核？

1. #P b3lyp/6-31g opt
   
2. 0  3
   
3. 坐标

复制代码

请问一下用这种方法优化出来的T1，输出文件的哪一行对应的是T1的能量

yopoxia

小范范1989 发表于 2014-12-26 08:16
对，sob写的完全对，就是sob的意思。
荧光是从s1发出的，磷光是从t1发出的。
但是也不排除从高激发态发 ...

@sobereva
感谢您的回复。还想再追问几个问题。
您看一下我说的步骤对不。计算苯的荧光，这次就纯文字写了。
一、先高斯画苯的结构图，先算opt，命令#p opt b3lyp/6-31G
二、以第一步优化后的坐标，再计算freq，命令#p freq b3lyp/6-31G
三、以第二步优化后的坐标，计算荧光，命令#p b3lyp/6-31G opt td（singlets,root=1）
这里，有几个问题，
1、第三步还需要加opt关键词么？
2、这样结果里计算的就是荧光波长了？
3、如30楼的例子，Excited state 1的 f 值这么小，是不是这个只是计算上的荧光发射峰位置，不是实验的，实验的应该是Excited state 2，我应该自动过滤掉前面 f 特别小的，而把下面大的作为Excited state 1来处理呢？
4、如果是计算吸收峰位置，是第三步的命令换为#p b3lyp/6-31G  td=（nstates=20），第一步和第二步不变，这样就可以了吧，不需要再加opt了吧
5、还一个问题是#p b3lyp/6-31G opt td（singlets,root=1），用这个默认显示3个td，有必要再增加显示的数目么？应该写什么关键词呢？

yopoxia

sobereva 发表于 2014-12-25 23:42
荧光算的是singlets，root=1。root>1不符合kasha规则，triplet则对应的是磷光
最后一轮激发态优化产生的 ...

sob姐，您能再回答我一下35楼的问题么？非常感谢您。
这几天学到了很多知识~

小范范1989

去离子水 发表于 2014-12-26 14:39
请问一下用这种方法优化出来的T1，输出文件的哪一行对应的是T1的能量

你好，用0  3的方式优化出来的是t1的结构，然而t1的能量应该是输出文件的HF能量。我感觉是这样的。刚刚我说错了
还在问一句，你是要找s1-t1的能差吗？

小范范1989

yopoxia 发表于 2014-12-26 14:27
嗯 是的 用的是高斯09  
8核  16核都可以

偶，你要是方便，能加qq吗？我们这边计算资源太缺了，有很多测算不了，你要是方便能帮助算一算吗？没事，要是你们也紧张，就算了，我再等等就行，哈哈

去离子水

小范范1989 发表于 2014-12-26 14:49
你好，用0  3的方式优化出来的是t1的结构，然而t1的能量应该是输出文件的HF能量+第一个td的能量。
理解 ...

不太理解，你就拿你在http://emuch.net/bbs/viewthread.php?tid=8072030 第二种方法.log举例好了

小范范1989

去离子水 发表于 2014-12-26 14:54
不太理解，你就拿你在http://emuch.net/bbs/viewthread.php?tid=8072030 第二种方法.log举例好了

不好意思，刚刚我的不对。我改了。不好意思

小范范1989

yopoxia 发表于 2014-12-26 14:42
@sobereva
感谢您的回复。还想再追问几个问题。
您看一下我说的步骤对不。计算苯的荧光，这次就纯文字写 ...

1：第三步是优化s1，我感觉，要加opt。
2：你第三步优化出来的文件，默认有三个td、第一个td对应的波长就是荧光波长。
3：这个问题我也发现了，我还不知道怎么弄。麻烦sob大神回答后，通知我一声，谢谢你。
4：对，计算吸收波长，就是相当于基于s0的结构做td。你说的这个就是
5：#p b3lyp/6-31G opt td（singlets,root=1，nstates=10），这个默认3个，是nstates=3，要想显示几个就把nstates=n。

sobereva

yopoxia 发表于 2014-12-26 14:42
@sobereva
感谢您的回复。还想再追问几个问题。
您看一下我说的步骤对不。计算苯的荧光，这次就纯文字写 ...

计算荧光，对这么简单的体系只需要一步，直接写#p b3lyp/6-31G* opt td就够了，默认就是singlet,root=1。最后一步的第一激发能就是荧光对应的垂直发射能

4 对。但freq那步是多余的，freq不会对坐标进行任何改变。
5 通常没必要增加，只会徒增时间。nstates数越大每一步激发态优化的耗时也会越长。

PS：6-31G太烂，计算光谱的底限是6-31G*

sobereva

yopoxia 发表于 2014-12-26 14:26
嗯  非常安心sob老师的回复。
还有个疑问，希望sob cc能再解释一下。
有的时候最后一轮的Excited State ...

这没法称为“计算上的荧光最大峰位置”。
就算振子强度很小，或者哪怕为0，由于振动耦合，还是能从S1回到基态的。只不过强度会很弱。
通常是从S1返回S0，但碰上当前情况，也可能S2返回S0的情况不可忽略，因此实验上最大峰与哪个激发态的辐射跃迁相对应，不好说。

yopoxia

小范范1989 发表于 2014-12-26 15:09
1：第三步是优化s1，我感觉，要加opt。
2：你第三步优化出来的文件，默认有三个td、第一个td对应的波长 ...

嗯  sob姐已经回复我了。
您看一下吧。
还有  您给我回复的，我还有个疑问。这个基于s0  s1的优化做td，选择用s0 还是s1，是靠的计算的关键词，就是说，我弄好了初期的opt freq。然后如果我写“#p b3lyp/6-31G*  td=（nstates=20）”就是基于s0的td，是吸收能；如果做好初期opt freq，我写“#p b3lyp/6-31G* opt td（singlets,root=1）”就是基于s1的td，是计算的发射。是这个样子么？

yopoxia

sobereva 发表于 2014-12-26 16:13
这没法称为“计算上的荧光最大峰位置”。
就算振子强度很小，或者哪怕为0，由于振动耦合，还是能从S1回 ...

嗯  非常感谢sob老师！
我还有个疑问。这个基于s0还是s1的优化做td，选择用s0 还是s1，是靠的计算的关键词么？
就是说，我弄好了初期的opt freq。然后如果我写“#p b3lyp/6-31G*  td=（nstates=20）”就是基于s0的td，是吸收能；如果做好初期opt freq，我写“#p b3lyp/6-31G* opt td（singlets,root=1）”就是基于s1的td，是计算的发射。是这个样子么？