求助，激发态外重组能计算

shyshy · 发表于 Post on 2020-10-11 21:22:36

看了sob老师《乱谈外重组能的计算》的文章，我计算了隐式溶剂模型下S1-T1的外重组能，在4-1时得到的外重组能是负值，下面列出我输入和输出，求解是不是输入有问题？感谢~
1---------------------
%nproc=24
%mem=50GB
%chk=1.chk
#p b3lyp/6-31g*  SCRF(solvent=toluene,read) td(singlets,root=1,nstates=10)
1

0 1

noneq=write
2------------------
%nproc=24
%mem=50GB
%oldchk=1.chk
%chk=2.chk
#p b3lyp/6-31g*  SCRF(solvent=toluene,read)
2
0 3

noneq=read
3------------------
%nproc=24
%mem=50GB
%chk=3.chk
#p b3lyp/6-31g*  SCRF(solvent=toluene,read)
3
0 3

noneq=write
4-------------------
%nproc=24
%mem=50GB
%oldchk=3.chk
%chk=4.chk
#p b3lyp/6-31g*  SCRF(solvent=toluene,read)  td(singlets,root=1,nstates=10)
4
0 1

noneq=read

其中2-3=0.0017eV  4-1=-0.0014eV  ，4-1是负的想请问大家错在哪里了呢？

sobereva · 发表于 Post on 2020-10-12 01:54:54

noneq=write和noneq=read在当前这种线性响应模型（LR）表现溶剂效应的激发态计算时无效（但也不会报错），必须改用特定态模型（scrf里用ExternalIteration的情况）。这里有说：
Gaussian中用TDDFT计算激发态和吸收、荧光、磷光光谱的方法
http://sobereva.com/314（http://bbs.keinsci.com/thread-2413-1-1.html）

shyshy · 发表于 Post on 2020-10-12 10:32:59

sobereva 发表于 2020-10-12 01:54
noneq=write和noneq=read在当前这种线性响应模型（LR）表现溶剂效应的激发态计算时无效（但也不会报错）， ...

明白啦！感谢sob老师~

Kalinite · 发表于 Post on 2021-2-9 00:48:27

sobereva 发表于 2020-10-12 01:54
noneq=write和noneq=read在当前这种线性响应模型（LR）表现溶剂效应的激发态计算时无效（但也不会报错）， ...

sob老师您好，对于用特定态考察S1-T1的情况，第三步计算是优化T1的，但是在外重组能计算用非平衡溶剂效应时，第三步却不优化，这个是为什么呢？

sobereva · 发表于 Post on 2021-2-9 02:51:43

Kalinite 发表于 2021-2-9 00:48
sob老师您好，对于用特定态考察S1-T1的情况，第三步计算是优化T1的，但是在外重组能计算用非平衡溶剂效应 ...

本来计算的目的就是计算外重组能

Kalinite · 发表于 Post on 2021-2-10 02:51:18

sobereva 发表于 2021-2-9 02:51
本来计算的目的就是计算外重组能

冒昧追问，那如果要算的是阳离子VEA过程，就是对VEA过程一样做吗？
感谢sob老师！

sobereva · 发表于 Post on 2021-2-10 04:58:43

Kalinite 发表于 2021-2-10 02:51
冒昧追问，那如果要算的是阳离子VEA过程，就是对VEA过程一样做吗？
感谢sob老师！

是

Osmanthus · 发表于 Post on 2022-3-28 01:52:26

在读过sob老师的博文《乱谈外重组能的计算》(http://sobereva.com/333)后，结合上述的问题仍然有些疑惑：
对于文中给出溶液下外重组能的计算示例，我的理解如下
(1) 优化基态的分子结构，计算基态的分子的电子结构与此时溶剂的状态，并保存此时的溶剂状态
(2) 基于基态的分子结构，读取(1)中保存的溶剂状态，计算激发态的分子的电子结构
(3) 基于基态的分子结构，计算激发态的分子的电子结构与此时溶剂的状态，并保存此时的溶剂状态
(4) 基于基态的分子结构，读取(3)中保存的溶剂状态，计算基态的分子电子结构
E(2)-E(3) 垂直激发前后分子电子结构的改变，所引起溶剂状态改变的能量差
E(4)-E(1) 从分子的FC结构发射回到基态极小点前后分子电子结构的改变，所引起溶剂状态改变的能量差
我的疑惑是在考虑电子态变化前后引发环境分子几何结构弛豫造成的能量变化时，对于发生过程而言，选取垂直发射是否更为合理，因而对于(1)与(4)两步的计算中是否采用激发态极小点与其垂直发射到基态的FC结构？

wzkchem5 · 发表于 Post on 2022-3-28 03:19:16

Osmanthus 发表于 2022-3-27 18:52
在读过sob老师的博文《乱谈外重组能的计算》(http://sobereva.com/333)后，结合上述的问题仍然有些疑惑：
...

外重组能指的是溶剂对重组能的贡献，从定义上就把溶质自己分子结构变化引入的重组能（也就是内重组能）排除在外了。你那种算法算出来的是总重组能

Osmanthus · 发表于 Post on 2022-3-29 00:11:04

补充一个问题，如果按照上述方法算出来的就是总重组能，那么和内外重组能相加得到的总重组能(内重组能通过四点法得到，外重组能按照sob老师的博文计算)有无区别?

Kalinite · 发表于 Post on 2022-3-29 09:16:55

有张配图可能会讨论的更清晰一点。Delayed表示该点是考虑了非平衡溶剂效应，溶剂慢部分未响应的情况。

基于前面大家的讨论，以S1->T1这一过程为例，是否应该如下：
总重组能=E1
内重组能=E1-E4
外重组能=E4
而卢老师博文里提到的电离反应外重组能是E3+E4，这两种情况的区别是什么呢？

xxtnenu · 发表于 Post on 2024-6-21 13:22:13

sobereva 发表于 2020-10-12 01:54
noneq=write和noneq=read在当前这种线性响应模型（LR）表现溶剂效应的激发态计算时无效（但也不会报错）， ...

请问sob老师，在博文333中 2 溶液下外重组能的计算隐式溶剂模型下计算外重组能为什么关键词中没有写隐式溶剂模型呢？以及λext(I)为什么是E2和E3计算得到而不是图中E1和E4计算得到？

sobereva · 发表于 Post on 2024-6-22 02:31:21

xxtnenu 发表于 2024-6-21 13:22
请问sob老师，在博文333中 2 溶液下外重组能的计算隐式溶剂模型下计算外重组能为什么关键词中没有写隐 ...

乱谈外重组能的计算
http://sobereva.com/333

第2节的关键词明明写了scrf

笔误。但不影响最后的外重组能的结果

xxtnenu · 发表于 Post on 2024-6-22 02:58:10

sobereva 发表于 2024-6-22 02:31
乱谈外重组能的计算
http://sobereva.com/333

收到，谢谢老师回复，还请问sob老师，此处的溶剂模型可以用SMD溶剂模型吗？还是用默认即可。

sobereva · 发表于 Post on 2024-6-22 06:04:43

xxtnenu 发表于 2024-6-22 02:58
收到，谢谢老师回复，还请问sob老师，此处的溶剂模型可以用SMD溶剂模型吗？还是用默认即可。

默认的

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[新手求助] 求助，激发态外重组能计算

评分 Rate