计算化学公社

标题: NTO分析的一个问题的咨询。 [打印本页]

作者
Author: 小范范1989 时间: 2015-4-26 09:24
标题: NTO分析的一个问题的咨询。
我看了sob的帖子：http://hi.baidu.com/sobereva/item/16a9e1d833658019d68ed0c3
我说说我遇到的问题：
1：背景：我想做s1到s0的跃迁，按说经过弛豫，应该是一个电子在lumo上，一个电子在homo上。所以我应该看优化s0结构后的lumo和homo。
2：为了近似的研究，我看了s0到s1的跃迁，基于s0做了td，但是这个td之后，我发现s0-s1的跃迁大部分是homo-1到lumo的跃迁。
3：为了把s0-s1的跃迁和s1-s0的跃迁的轨道统一起来，我是不是可以做NTO分析？也就是都看最高的占据NTO和最低的未占据NTO，这样，我来分析s0-s1跃迁或者是s1-s0跃迁，涉及到的NTO都是一样的？我不确定，所以来咨询大家一下。
4：我基于s0做了NTO得到了下面的文件，
Alpha  occ. eigenvalues -- 0.00008 0.00015 0.00016 0.00028 0.00042
Alpha  occ. eigenvalues -- 0.00192 0.00418 0.99263
Alpha virt. eigenvalues -- 0.99263 0.00418 0.00192 0.00042 0.00028
通过学习sob的帖子，这个我知道怎么看，但是不太理解为什么还有这个输出
Excited State 1:    Singlet-A    4.8054 eV  258.01 nm  f=0.0000  <S**2>=0.000
   42 -> 49       0.10386
   44 -> 49       0.61486
   44 -> 54       -0.27473
这个再怎么办？看NTO的时候，这个就可以不看了吗？
谢谢大家的指点。

作者
Author: blueyangliu 时间: 2015-4-26 10:30
s0-s1的跃迁和s1-s0的跃迁的轨道统一起来?这是什么意思

作者
Author: 小范范1989 时间: 2015-4-26 10:47

blueyangliu 发表于 2015-4-26 10:30
s0-s1的跃迁和s1-s0的跃迁的轨道统一起来?这是什么意思

也就是这个意思：
1：so到s1的跃迁，在这里是homo-1到lumo。
2：s1到s0的跃迁，我自己认为的是应该是lumo到homo的跃迁。
所以，我想能不能表示成一种跃迁，比如说s0-s1的跃迁是H-L。s1-s0的跃迁也是L-H，这样最好，但是这里是H-1到L。所以我就想能不能通过NTO来做，因为NTO表示的不都是最高的占据NTO到最低的为占据NTO？
这样不就能表示成一样的形式了？我表达能力不好，加上我的思路可能乱，说的不清楚，见谅、

ps：咨询您个问题：

我想咨询一个问题，高斯做NTO的时候，用CIS是这样写：cis/6-31+G* pop(saveNTO,NTO) density=transition=2，
如果说用b3lyp，是不是需要加上td的选项？就是b3lyp/6-31+G* td（nstates=6，singlet，root1） pop(saveNTO,NTO) density=transition=2，？谢谢指点

作者
Author: blueyangliu 时间: 2015-4-26 14:13
加上td的

作者
Author: 小范范1989 时间: 2015-4-26 15:48

blueyangliu 发表于 2015-4-26 14:13
加上td的

嗯，我已经算上了，谢谢你的回复。

作者
Author: sobereva 时间: 2015-4-26 19:09
都表示成H-L就是人为作假了。
做NTO分析后虽然都可以近似表示为本征值最大的占据NTO和非占据NTO的跃迁，但是S1-S0时和S0-S1时NTO的形状会明显的不同。

以NTO来考察的话就不用管诸如42 -> 49这些基于MO的表示了。

TDDFT时就是那么写

作者
Author: 小范范1989 时间: 2015-4-27 08:21

sobereva 发表于 2015-4-26 19:09
都表示成H-L就是人为作假了。
做NTO分析后虽然都可以近似表示为本征值最大的占据NTO和非占据NTO的跃迁，但 ...

嗯，对，sob指点的有道理，谢谢sob

作者
Author: 小范范1989 时间: 2015-5-21 12:03

sobereva 发表于 2015-4-26 19:09
都表示成H-L就是人为作假了。
做NTO分析后虽然都可以近似表示为本征值最大的占据NTO和非占据NTO的跃迁，但 ...

sob老师好，我想咨询您一下。做NTO的时候，
1：做s1的NTO是不是这样写 b3lyp/6-31+G* td（nstates=6，singlet，root=1） pop(saveNTO,NTO) density=transition=1
2：如果我想做t1的NTO是不是这样写 b3lyp/6-31+G* td（nstates=6，triplet，root=1） pop(saveNTO,NTO) density=transition=1
下面的咨询多重度是不是写 1？
谢谢sob，我对这个T1的NTO不太会写，所以来咨询sob老师。谢谢sob老师指点

作者
Author: sobereva 时间: 2015-5-21 19:39

小范范1989 发表于 2015-5-21 12:03
sob老师好，我想咨询您一下。做NTO的时候，
1：做s1的NTO是不是这样写 b3lyp/6-31+G* td（nstates=6， ...

你说的都对

作者
Author: 小范范1989 时间: 2015-5-22 07:56

sobereva 发表于 2015-5-21 19:39
你说的都对

谢谢sob老师的指点。

作者
Author: 石头 时间: 2016-12-28 09:47

小范范1989 发表于 2015-5-21 12:03
sob老师好，我想咨询您一下。做NTO的时候，
1：做s1的NTO是不是这样写 b3lyp/6-31+G* td（nstates=6， ...

您好，请问您是基于优化的S1态和T1态构型进行这个NTO 计算的是吗？那么T1态的优化是采用 0 3还是TD计算的呢？

作者
Author: 小范范1989 时间: 2016-12-28 10:09

石头发表于 2016-12-28 09:47
您好，请问您是基于优化的S1态和T1态构型进行这个NTO 计算的是吗？那么T1态的优化是采用 0 3还是TD计算的 ...

基于S0的结构做的。
b3lyp/6-31+G* td（nstates=6，triplet，root=1） pop(saveNTO,NTO) density=transition=1

作者
Author: 石头 时间: 2016-12-28 11:37

小范范1989 发表于 2016-12-28 10:09
基于S0的结构做的。
b3lyp/6-31+G* td（nstates=6，triplet，root=1） pop(saveNTO,NTO) density=trans ...

也就是说您这里看的是基于S0构型下，S0到S1跃迁和S0到T1跃迁的NTO。那么如果我要看基于S1构型下和基于T1构型下,S1到S0跃迁和T1到S0跃迁，就要基于优化的S1和T1几何构型了。这样的话计算方法还是一样是吗？感谢您的回复！

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