计算化学公社

标题: 高斯优化激发态出错 [打印本页]

作者
Author: 小范范1989 时间: 2014-11-30 10:46
标题: 高斯优化激发态出错
高斯优化s1态出错，我也查过，说加上direct。但是我也加了，还是出错，不知道怎么改？谢谢

错误为

作者
Author: 小范范1989 时间: 2014-11-30 10:46
标题: ]高斯优化t1也是出错
我优化t1也是出错，
输入文件：
%chk=test.chk
%mem=2000MB
%nprocs=2
#p b3lyp/6-31* opt td(direct,triplets,nstates=6,root=1)

test

0 1
N                1.55801000 0.19783500 -0.79778700
C                1.09085500 -0.95453600 -1.48926800
C                2.14146600 -1.89758300 -1.54718200
C                -0.13547400 -1.20956500 -2.10922100

错误
NMat=    8 NSing=    0.
Cannot handle 2e integral symmetry, ISym2E=1.
CISAX:  IP=  1 NPass= 1 NMax=    8.
CISAX will form    8 AO SS matrices at one time.
NMat=    8 NSing=    0.
Root    1 has converged.
Root    2 has converged.
Root    3 has converged.
Root    4 has converged.
Root    5 not converged, maximum delta is 0.000011281594252
Root    6 not converged, maximum delta is 0.000006103042082
Excitation Energies [eV] at current iteration:
Root    1 :    1.916757579771638 Change is -0.000000000000918
Root    2 :    2.060965374014521 Change is -0.000000000000801
Root    3 :    2.576906613973796 Change is -0.000000000107453
Root    4 :    2.603692086423475 Change is -0.000000000049297
Root    5 :    2.632340173082165 Change is -0.000000002306130
Root    6 :    2.637826619814619 Change is -0.000000000652924
Iteration 13 Dimension 136 NMult 136
Cannot handle 2e integral symmetry, ISym2E=1.
CISAX:  IP=  1 NPass= 1 NMax=    4.
CISAX will form    4 AO SS matrices at one time.
NMat=    4 NSing=    0.
Cannot handle 2e integral symmetry, ISym2E=1.
CISAX:  IP=  1 NPass= 1 NMax=    4.
CISAX will form    4 AO SS matrices at one time.
NMat=    4 NSing=    0.

作者
Author: 小范范1989 时间: 2014-11-30 10:47
标题: 高斯优化最低三重态出错，而且出错还不相同
输入表头
%mem=2000MB
%nprocs=2
#p b3lyp/6-31* opt td(direct,triplets,nstates=6,root=1)

test

0 1
N                1.50647400 1.80483400 0.04491200
C                1.22289100 1.69315200 1.41808500

输出
Cannot handle 2e integral symmetry, ISym2E=1.
CISAX:  IP=  1 NPass= 1 NMax= 12.
CISAX will form 12 AO SS matrices at one time.
NMat= 12 NSing=    0.
Root    1 not converged, maximum delta is 0.298200414274617
No map to state    2
No map to state    3
No map to state    4
New state    5 was old state    2
Root    5 not converged, maximum delta is 0.672121386957955
No map to state    6
You need to solve for more vectors in order to follow this state.
Error termination via Lnk1e in /home/software/g09/l914.exe at Sat Nov 29 21:23:51 2014.
Job cpu time:  3 days  1 hours  9 minutes 25.1 seconds.
File lengths (MBytes):  RWF= 2932 Int=    0 D2E=    0 Chk= 321 Scr=    1

作者
Author: 小范范1989 时间: 2014-11-30 10:48
标题: 优化t1结构，采用的是电荷为0 ，自旋为3.也是出错
#p b3lyp/6-31* opt

test

0 3
N                1.90437600 -0.13917600 -0.93457000
C                1.38502800 -1.42638600 -0.65539200

错误：
NTrRot= -1 NTRed= 652 NAtoms= 108 NSkip= 334 IsLin=F
Error in internal coordinate system.
Error termination via Lnk1e in /home/software/g09/l103.exe at Fri Nov 28 21:28:37 2014.
Job cpu time:  0 days  5 hours  2 minutes 45.9 seconds.
File lengths (MBytes):  RWF= 1195 Int=    0 D2E=    0 Chk=    30 Scr=    1

对于错误，我已经看过帖子，说加上OPT=Cartesian，可以解决，但是，我还是出错。

作者
Author: sobereva 时间: 2014-11-30 11:52
都是同一类问题，请发在同一个帖子里，不要分成四个单独发，否则看帖和回答都比较麻烦。我已经进行合并了

作者
Author: sobereva 时间: 2014-11-30 12:09
1、3L：看此贴http://bbs.keinsci.com/forum.php?mod=viewthread&tid=299
2L：输出中并没有出错信息。如果还在算就继续等着，如果已经停了就重算，还不行就换版本或者改方法，或者直接给高斯中文客服发邮件
4L：改用opt=cartesian

作者
Author: 小范范1989 时间: 2014-11-30 15:35

sobereva 发表于 2014-11-30 11:52
都是同一类问题，请发在同一个帖子里，不要分成四个单独发，否则看帖和回答都比较麻烦。我已经进行合并了

嘿嘿，麻烦sob了。我以后会注意的。:lol

作者
Author: 小范范1989 时间: 2014-12-7 15:40

sobereva 发表于 2014-11-30 12:09
1、3L：看此贴http://bbs.keinsci.com/forum.php?mod=viewthread&tid=299
2L：输出中并没有出错信息。如果 ...

sob你好，我出现这样的错误，我按照你说的，设置nstates改大，但是还是出错，不管是单重态还是三重态。谢谢sob
输入%chk=test.chk
%mem=2GB
%nprocs=2
#p b3lyp/6-31* opt td(direct,triplets,nstates=20,root=1)

test

0 1
N                1.50647400 1.80483400 0.04491200
C                1.22289100 1.69315200 1.41808500
C                2.22262000 0.88887400 2.02102800

输出
New state    1 was old state    3
Root    1 not converged, maximum delta is 0.571065717764176
New state    2 was old state    10
Root    2 not converged, maximum delta is 0.624581561156494
New state    3 was old state    1
Root    3 not converged, maximum delta is 0.020512414872001
New state    4 was old state    2
Root    4 not converged, maximum delta is 0.289690504679012
New state    5 was old state    4
Root    5 not converged, maximum delta is 0.436164850671245
New state    6 was old state    5
Root    6 not converged, maximum delta is 0.513899117437082
No map to state    7
No map to state    8
No map to state    9
No map to state    10
New state    11 was old state    6
Root    11 not converged, maximum delta is 0.455626045840689
No map to state    12
New state    13 was old state    8
Root    13 not converged, maximum delta is 0.309748276085905
No map to state    14
You need to solve for more vectors in order to follow this state.
Error termination via Lnk1e in /home/software/g09/l914.exe at Sun Dec  7 05:49:12 2014.
Job cpu time:  9 days 11 hours 37 minutes 46.1 seconds.
File lengths (MBytes):  RWF= 6110 Int=    0 D2E=    0 Chk= 844 Scr=    1

作者
Author: sobereva 时间: 2014-12-7 15:43
又是三重态又是root=1，直接用0 3计算就行了，不需要TDDFT，又快又不会出现这些问题

作者
Author: 小范范1989 时间: 2014-12-7 15:56

sobereva 发表于 2014-12-7 15:43
又是三重态又是root=1，直接用0 3计算就行了，不需要TDDFT，又快又不会出现这些问题

嗯，好的，谢谢sob。我用两种方式都做了，确实0  3 出来的快多了，也没有错误。
但是，我优化s1的时候还是这样的错误。
我也是按照，nstates扩大，都到20了。
%chk=test.chk
%mem=2GB
%nprocs=2
#p b3lyp/6-31* opt=(MaxStep=10) td(nstates=20,singlets,root=1,direct)

test

0 1
N                1.55801000 0.19783500 -0.79778700
C                1.09085500 -0.95453600 -1.48926800

还是说you need tosolve.....

作者
Author: sobereva 时间: 2014-12-7 16:35
nstates再设大一些。maxstep不仅设小，还应该加上notrust以免自动在优化中更新置信半径，否则步长仍可能超过maxstep的值。

作者
Author: 小范范1989 时间: 2014-12-7 17:11

sobereva 发表于 2014-12-7 16:35
nstates再设大一些。maxstep不仅设小，还应该加上notrust以免自动在优化中更新置信半径，否则步长仍可能超 ...

谢谢sob。thanks:lol

作者
Author: 小范范1989 时间: 2014-12-8 22:35

sobereva 发表于 2014-12-7 16:35
nstates再设大一些。maxstep不仅设小，还应该加上notrust以免自动在优化中更新置信半径，否则步长仍可能超 ...

sob，再问您个问题。还是激发态的问题。谢谢你
优化第二激发态后计算的结果发现第二激发态的能量比第一激发态能量还低，正常吗，下面分别是优化第一和第二激发态的结果。

优化激发态（root=2）：
Excited State 1:    Singlet-A    2.4155 eV  513.29 nm  f=0.3072  <S**2>=0.000
   152 ->156       0.10882
   153 ->156       0.10052
   154 ->156       0.29460
   155 ->156       -0.62377

Excited State 2:    Singlet-A    1.6508 eV  751.06 nm  f=0.0014  <S**2>=0.000
   153 ->156       0.21356
   154 ->156       0.58364
   155 ->156       0.31673

Excited State 3:    Singlet-A    2.6231 eV  472.67 nm  f=0.1146  <S**2>=0.000
   148 ->156       -0.13364
   153 ->156       0.65309
   154 ->156       -0.22779
，这个出来的第二个激发能比第一个还小。这是怎么理解？优化的是第二激发态。
下面是我的理解
在基于s2的结构做td，出来的第一个td，也就是2.4155ev，指的是，从基态势能面，s2的结构处，直接td到s2势能面的最低处的激发能。
而这个1.6508ev，指的是，基态势能面，s2的结构处，直接td，到s1势能面上的激发能（注意，不是s1时能面的的最低处，而是，垂直激发）。

作者
Author: sobereva 时间: 2014-12-8 23:38
这是因为你选定对root=2进行优化，后来势能面发生交叉，能量高低顺序发生改变，目前的结构下曾经的root=2已经不是激发能第二高的态了，但是优化一直是对这个态来优化的。

作者
Author: flyskyzff 时间: 2014-12-11 09:37

小范范1989 发表于 2014-12-7 15:56
嗯，好的，谢谢sob。我用两种方式都做了，确实0 3 出来的快多了，也没有错误。
但是，我优化s1的时候还 ...

请问优化t1,你这里说的“直接0 3计算就行了，不需要TDDFT”的意思是不是route部分只要# b3lyp/6-31* opt，然后电荷为0 ，自旋为3就可以了？

作者
Author: 小范范1989 时间: 2014-12-11 09:45
本帖最后由小范范1989 于 2014-12-11 09:48 编辑

flyskyzff 发表于 2014-12-11 09:37
请问优化t1,你这里说的“直接0 3计算就行了，不需要TDDFT”的意思是不是route部分只要# b3lyp/6-31* opt ...

对，也就是使用DFt的方式，0 3
http://emuch.net/bbs/viewthread.php?tid=8072030

作者
Author: 小范范1989 时间: 2014-12-16 15:52

sobereva 发表于 2014-12-7 15:43
又是三重态又是root=1，直接用0 3计算就行了，不需要TDDFT，又快又不会出现这些问题

sob好，之前发的帖子，是优化激发态的，按照你的方法，我优化出来了最低三重态，使用的0 3这样的方法，
但是，我有个问题，就是我优化出来后，然后td（50-50），结果如下，我不知道那个是激发态的三重态，哪个是单重态。谢谢sob
同时，我想的方法是，要是这样不能解决，我们可以在做完td（50-50）之后，在重新做td。这样，两个不同的，就是t1激发态s1或者是s2 s3....的。行得通吧？
谢谢sob

作者
Author: sobereva 时间: 2014-12-16 17:25
你看<S**2>就能判断，见http://bbs.keinsci.com/forum.php?mod=viewthread&tid=111
貌似你是在0 3下做的TDDFT，因此造成<S**2>不为整数。不建议这样，TDDFT应当在0 1下做。

作者
Author: 小范范1989 时间: 2014-12-16 19:50

sobereva 发表于 2014-12-16 17:25
你看就能判断，见http://bbs.keinsci.com/forum.php?mod=viewthread&tid=111
貌似你是在0 3下做的TDDFT， ...

谢谢sob。我还是有点迷惑，我是在 03的情况下，优化出来了t1态，而不是使用0 1的方式优化出来t1，因为用0 1的方式，优化一直出错，而且也不准确，所以是在0 3的状态下优化出来的t1。
现在我在 03 优化出来的基础上，做的td（50-50），出现上面的结果，
按照您刚刚说的，TDDFT应该在 01下面做，那我应该怎么进行下一步？
是不是在 03 优化出来的基础上，然后采用0 3优化出来的坐标，然后设置电荷为0，自旋为1，在此基础上做td？
谢谢sob，

作者
Author: sobereva 时间: 2014-12-16 20:03
是的

作者
Author: 小范范1989 时间: 2014-12-16 20:39

sobereva 发表于 2014-12-16 20:03
是的

谢谢sob。我发现只要有不会的问题，来问sob，一定会得到解答。
sob不但知识面很广，而且为人和热情细心，谢谢sob。

作者
Author: 小范范1989 时间: 2014-12-18 15:11

sobereva 发表于 2014-12-16 20:03
是的

sob好，还是这个优化激发态的问题，我这不是也发过帖子，问过好多，也修改了好多，但是还是出问题，还请sob再次指点。
我的输入文件是：
%chk=test.chk
%mem=2GB
%nprocs=2
#p b3lyp/6-31g* opt td(nstates=10,singlets,root=1,direct)

test

0 1
N                1.55801000 0.19783500 -0.79778700
C                1.09085500 -0.95453600 -1.48926800
C                2.14146600 -1.89758300 -1.54718200
C                -0.13547400 -1.20956500 -2.10922100
C                1.96005300 -3.11900500 -2.20189300
C                -0.29383200 -2.43297000 -2.76144400
H                -0.93851000 -0.48827500 -2.08290200
C                0.73480500 -3.38468600 -2.80530000
H                2.76816500 -3.84409000 -2.24412800
.....................（还有很多坐标）

错误时：
NMat= 20 NSing= 20.
New state    1 was old state    7
New state    2 was old state    8
New state    3 was old state    1
New state    4 was old state    6
New state    5 was old state    4
New state    6 was old state    3
New state    7 was old state    9
New state    8 was old state    5
New state    9 was old state    10
New state    10 was old state    2
Excitation Energies [eV] at current iteration:
Root    1 :    5.982060062346076
Root    2 :    6.383931700600978
Root    3 :    8.901901716469364
Root    4 :    9.700349738000840
Root    5 : 11.258633326878070
Root    6 : 11.308304626971730
Root    7 : 11.713305846440200
Root    8 : 12.248164103617310
Root    9 : 13.290151060085080
Root    10 : 13.512389675544020
Iteration    2 Dimension 40 NMult 40
Cannot handle 2e integral symmetry, ISym2E=1.
CISAX:  IP=  2 NPass= 1 NMax= 20.
CISAX will form 20 AO SS matrices at one time.
NMat= 20 NSing= 20.
Cannot handle 2e integral symmetry, ISym2E=1.
CISAX:  IP=  2 NPass= 1 NMax= 20.
CISAX will form 20 AO SS matrices at one time.
NMat= 20 NSing= 20.
New state    1 was old state    2
Root    1 not converged, maximum delta is 0.433400373447380
New state    2 was old state    1
Root    2 not converged, maximum delta is 0.528876387304794
No map to state    3
You need to solve for more vectors in order to follow this state.
Error termination via Lnk1e in /home/software/g09/l914.exe at Thu Dec 21 13:17:32 2006.
Job cpu time:  3 days  8 hours 11 minutes 34.4 seconds.
File lengths (MBytes):  RWF= 3198 Int=    0 D2E=    0 Chk= 442 Scr=    1

当时按照您说的，这样的错误 You need to solve for more vectors in order to follow this state，就是一个增加nstates，但是我都到20了，还在增加吗？而且，我看还有错误是不是不收敛？
Root    1 not converged, maximum delta is 0.433400373447380。
你说，是不是我应该重新调节一下我的分子构型？（但是，我的s1的初始构型，是优化的s0的构型）。或者是使用什么方法使其能计算出来？
谢谢sob，这个优化激发态，我已经弄了快2个月了，还是弄不出来，而且我还没有计算激发态的频率。这样耗时间太长了。所以一次一次的来麻烦sob。望谅解。谢谢sob

作者
Author: sobereva 时间: 2014-12-18 18:56
nstates再加到40还不行就不是这个原因了。
可以换换泛函、基组再试。
也可以先用比如CIS优化，如果能顺利收敛，然后再切换回TDDFT继续优化。
也可以减小步长上限，或者用别的几何优化方法试试。
总之没有万能办法，各种路子都试试。

作者
Author: 小范范1989 时间: 2014-12-18 19:00

sobereva 发表于 2014-12-18 18:56
nstates再加到40还不行就不是这个原因了。
可以换换泛函、基组再试。
也可以先用比如CIS优化，如果能顺利 ...

嗯，好的，晓得了，谢谢sob:lol

作者
Author: 小范范1989 时间: 2014-12-19 10:40
本帖最后由小范范1989 于 2014-12-19 10:41 编辑

回复错了sob，在我的另外一个帖子，sorry。看着好几个您的回帖，一下弄乱了，sorry

作者
Author: yopoxia 时间: 2014-12-25 13:11
本帖最后由 yopoxia 于 2014-12-25 13:13 编辑

您好，最近我也在学习高斯。想请教您几个问题。
我看了一些帖子，都说步骤如下，您看我理解的对不
1、先在高斯里画出结构，然后比如用这个计算“#p opt freq B3LYP/6-31G”，得到优化后的原子坐标。
2、再用第一步优化后的坐标，写计算命令“#p B3LYP/6-31G td=(nstates=10)”
后面的我就不懂了。。。。还有root什么的，s1 t1就晕了。。。。第一激发态的结构是怎么得到的。。。
您能帮我解答一下么。。。具体的步骤是什么。。。

作者
Author: 小范范1989 时间: 2014-12-25 15:16

yopoxia 发表于 2014-12-25 13:11
您好，最近我也在学习高斯。想请教您几个问题。
我看了一些帖子，都说步骤如下，您看我理解的对不
1、先 ...

你好，很高兴和你讨论问题。
1：这个root=n，我的感觉，就是对那个态感兴趣，比如说，我要优化第一单重激发态，那么就要写（singlets，root=1）。如果说对第二个三重态感兴趣，（triplets，root=2）
2：s1代表第一单重激发态，是分子的一种状态，和基态差不多意思。同样t1代表第一三重激发态。
3：做法：#p b3lyp/6-31g opt 这样优化出来基态坐标。然后使用基态坐标，该表头为#p b3lyp/6-31g opt td（singlets,root=1）这样就是优化s1。
#p b3lyp/6-31g opt td（triplets,root=1）这样是t1。
4：如果说计算出错，可以参考我的那几个帖子。有解决的方式，比如说加入direct等等

作者
Author: yopoxia 时间: 2014-12-25 21:42

小范范1989 发表于 2014-12-25 15:16
你好，很高兴和你讨论问题。
1：这个root=n，我的感觉，就是对那个态感兴趣，比如说，我要优化第一单重 ...

嗯  非常感谢您的热心回复。
不过还是有地方不明白，希望您能继续解答，谢谢！
1、如果我是要算荧光的话，是选择哪个激发态呢？是（singlets，root=1）还是（singlets，root=2）还是其他或者（triplets，root=1）（triplets，root=2）
2、您说的第三条我好像明白了。我下面说例子，您看我理解的对不
比如我想算水的。
高斯里面画好，得到gif如下%chk=H2O.chk
%mem=6GB
%nprocshared=8
#p opt b3lyp/6-31G

H2O

0 1
O                X1 Y1 Z1
H                X2 Y2 Z2
H                X3 Y3 Z3

然后再频率计算计算，取优化好的坐标，再计算，计算所用gjf内容如下
%chk=H2O.chk
%mem=6GB
%nprocshared=8
#p freq b3lyp/6-31G

H2O

0 1
O                X1a Y1a Z1a
H                X2a Y2a Z2a
H                X3a Y3a Z3a

这样一来，优化opt和频率freq计算就做好了。
然后再用这个最新的坐标，再计算，gjf内容如下写
%chk=H2O.chk
%mem=6GB
%nprocshared=8
#p b3lyp/6-31g opt td（singlets,root=1）

H2O

0 1
O                X1b Y1b Z1b
H                X2b Y2b Z2b
H                X3b Y3b Z3b

我的疑问是，上面计算的结果，就是最后一轮最大f对应的就是发射峰的位置了么？不是激发能？这个（singlets，root=1）还是（singlets，root=2）还是其他或者（triplets，root=1）（triplets，root=2）是怎么确定写哪个呢？

作者
Author: sobereva 时间: 2014-12-25 23:42

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

荧光算的是singlets，root=1。root>1不符合kasha规则，triplet则对应的是磷光
最后一轮激发态优化产生的第一激发能就是荧光最大峰位置的近似。f对应的是强度而非峰的位置。

作者
Author: 小范范1989 时间: 2014-12-26 08:16

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

对，sob写的完全对，就是sob的意思。
荧光是从s1发出的，磷光是从t1发出的。
但是也不排除从高激发态发出荧光。
但是，一般都是从s1发错的。
所以，我们算荧光发射的时候，都是td（singlets，root=1），
singlets代表单重态，root=1代表第一个，所以加起来就是s1。
那么，最后一步，你在做#p b3lyp/6-31g opt td（singlets,root=1）的时候，计算结果会有3个td（默认的）
找到第一个td（一定是文件的最后的那个，可以从文件最低往上找），对应的能量，转化为波长，就是荧光波长，看看下面的图

作者
Author: 小范范1989 时间: 2014-12-26 10:46

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

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

作者
Author: yopoxia 时间: 2014-12-26 14:26

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.

作者
Author: yopoxia 时间: 2014-12-26 14:27

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

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

作者
Author: 去离子水 时间: 2014-12-26 14:39

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

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

复制代码

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

作者
Author: yopoxia 时间: 2014-12-26 14:42
本帖最后由 yopoxia 于 2014-12-26 14:45 编辑

小范范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，有必要再增加显示的数目么？应该写什么关键词呢？

作者
Author: yopoxia 时间: 2014-12-26 14:46

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

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

作者
Author: 小范范1989 时间: 2014-12-26 14:49
本帖最后由小范范1989 于 2014-12-26 14:59 编辑

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

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

作者
Author: 小范范1989 时间: 2014-12-26 14:51

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

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

作者
Author: 去离子水 时间: 2014-12-26 14:54

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

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

作者
Author: 小范范1989 时间: 2014-12-26 15:00

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

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

作者
Author: 小范范1989 时间: 2014-12-26 15:09

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。

作者
Author: sobereva 时间: 2014-12-26 15:24

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*

作者
Author: sobereva 时间: 2014-12-26 16:13

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

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

作者
Author: yopoxia 时间: 2014-12-26 16:22
本帖最后由 yopoxia 于 2014-12-26 16:25 编辑

小范范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，是计算的发射。是这个样子么？

作者
Author: yopoxia 时间: 2014-12-26 16:29

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，是计算的发射。是这个样子么？

作者
Author: sobereva 时间: 2014-12-26 16:32

yopoxia 发表于 2014-12-26 16:29
嗯非常感谢sob老师！
我还有个疑问。这个基于s0还是s1的优化做td，选择用s0 还是s1，是靠的计算的关键 ...

是的

freq做不做无所谓，不是必需的。

作者
Author: yopoxia 时间: 2014-12-26 16:34

sobereva 发表于 2014-12-26 16:32
是的

freq做不做无所谓，不是必需的。

好的，谢谢sob老师。
那以后就不做freq了，做freq用大基组跑不动，还费时间。
再次感谢sob老师。:P

作者
Author: 小范范1989 时间: 2014-12-26 17:02

yopoxia 发表于 2014-12-26 16:22
嗯 sob姐已经回复我了。
您看一下吧。
还有您给我回复的，我还有个疑问。这个基于s0 s1的优化做td ...

偶，我的理解是，做吸收，是基于s0.的结构做td
做发射，基于优化出来的s1结构做td。（其实没必要做td了，应为，优化激发态，默认3个td。优化s1文件中的前三个td和你基于s1结构做td前三个是一样的。你可以看看）

作者
Author: yopoxia 时间: 2014-12-26 17:38

小范范1989 发表于 2014-12-26 17:02
偶，我的理解是，做吸收，是基于s0.的结构做td
做发射，基于优化出来的s1结构做td。（ ...

刚才sob姐也回复我了。我说的那个是对的。
你说的这个。。。我没理解。。。。

作者
Author: 小范范1989 时间: 2014-12-26 18:50

yopoxia 发表于 2014-12-26 17:38
刚才sob姐也回复我了。我说的那个是对的。
你说的这个。。。我没理解。。。。

只要会就可以，可能我的理解方式不对，哈哈。

作者
Author: yopoxia 时间: 2014-12-26 22:11

小范范1989 发表于 2014-12-26 18:50
只要会就可以，可能我的理解方式不对，哈哈。

http://bbs.keinsci.com/forum.php ... &extra=page%3D1
这个问题，能帮我看一下么？

作者
Author: ter20 时间: 2015-7-15 10:44
求问，TD里的Direct有什么用？我查了手册，CIS里有子关键字Direct，TD里没见着，所以求问
另外，加了Direct关键字后是不是会加快计算速度？谢谢！

作者
Author: sobereva 时间: 2015-7-15 11:46

ter20 发表于 2015-7-15 10:44
求问，TD里的Direct有什么用？我查了手册，CIS里有子关键字Direct，TD里没见着，所以求问
另外，加了Direc ...

CIS和TD计算时都需要MO积分。
CIS在g09中有三种方式处理积分，direct,MO,AO。默认是MO，先要做AO->MO积分变换并储存到硬盘上，速度和资源瓶颈主要在积分变换，之后迭代求解CI矩阵就很快了。AO是完全不储存MO积分，而是在迭代过程中每次用到MO积分的时候重新算AO，这样虽然比起SCF部分计算无任何额外资源消耗，但是在davidson迭代过程中颇慢。direct是介于二者中间，每次迭代也重算积分，但也储存一部分，所以davidson迭代虽慢，但比AO快得多得多。如果硬盘又快又有地方，MO总耗时最低；如果硬盘不够或者慢，用direct是最好选择。
TD就是用direct方式做的，没有其它可选余地

作者
Author: ter20 时间: 2015-7-15 11:48

sobereva 发表于 2015-7-15 11:46
CIS和TD计算时都需要MO积分。
CIS在g09中有三种方式处理积分，direct,MO,AO。默认是MO，先要做AO->MO积 ...

明白了，多谢Sob大神！

作者
Author: lijiaqi 时间: 2018-5-15 15:57

小范范1989 发表于 2014-12-18 15:11
sob好，还是这个优化激发态的问题，我这不是也发过帖子，问过好多，也修改了好多，但是还是出问题，还请s ...

你好，对于激发态优化，我出现了与你相同的结果，不收敛，赠加nstate或者减小基组都不可以，请问怎么改啊？谢谢指教

作者
Author: 小范范1989 时间: 2018-5-15 16:10

lijiaqi 发表于 2018-5-15 15:57
你好，对于激发态优化，我出现了与你相同的结果，不收敛，赠加nstate或者减小基组都不可以，请问怎么改啊 ...

http://sobereva.com/164
http://sobereva.com/61
尝试一下试试。
不过激发态优化确实很费劲，特别是大分子。当年研二的时候，就优化了几个激发态。但是，调节的过程中，学到了很多东西。
你可以看一下。

作者
Author: lijiaqi 时间: 2018-5-15 20:15

小范范1989 发表于 2018-5-15 16:10
http://sobereva.com/164
http://sobereva.com/61
尝试一下试试。

好的，谢谢。原先我用b3lyp,pbe0优化激发态可以出来结果，这次试的lc-wpbe方法，计算结果都没有收敛，我去试试，谢谢

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