计算化学公社

标题: S1，S0与T1非常靠近的体系在高斯里面该用什么方法计算 [打印本页]

作者
Author: BinWang 时间: 2025-9-2 15:45
标题: S1，S0与T1非常靠近的体系在高斯里面该用什么方法计算
本帖最后由 BinWang 于 2025-9-2 15:45 编辑

各位老师：
   大家好！我用高斯计算一个三重态上的抽取氢质子的反应：AH（+1价离子）和Cl-(氯离子)在三重态T1抽取氢质子的反应，T1上的IRC是平滑的（见下图的红色曲线），我把IRC上的点用Sob老师的程序提取出来以后在S0和S1态下计算单点能，评估这条IRC与激发态的交叉情况（见下图）。蓝色曲线是S0上的单点能（直接提取SCF Done这一行），紫色曲线是S1上的单点能（直接提取SCF Done这一行），可以看出大部分S0和S1的点是重合的，可能是S0的波函数不稳定。黑色曲线是我提取S1的能量点（直接提取 Total Energy, E(TD-HF/TD-DFT) = 这一行）。从这个图来看，我目前的体系S0和S1非常靠近。我现在用的方法是uwb97xd/genecp（Ir，SDD，其他主族元素6-31++G***）。我目前有如下几个问题：
   1：这种情况下S1上的单点能该提取”SCF Done“这一行，还是”total Energy, E(TD-HF/TD-DFT) = “这一行？按理说该是”total Energy, E(TD-HF/TD-DFT) = “这一行，但我这个结果显示”SCF Done“这一行更合理。
   2：我选用的泛函是不是有问题？这种体系不适合用wb97xd泛函计算？
   3：无论计算S0还是S1，都用uwb97xd，计算S1的时候是不是该用wb97xd，避免开壳层的计算？
另外补充一下，我之所以选择算T1上的IRC，是因为AH（+1价离子）的T1态比相应的S0态稳定2 kcal/mol。

作者
Author: wal 时间: 2025-9-2 16:09
你完全搞错了。你的情况直接做TD单点即可，此任务输出中SCF Done能量为S0能量，E(TD-HF/TD-DFT)为S1能量。

作者
Author: BinWang 时间: 2025-9-2 17:39

wal 发表于 2025-9-2 16:09
你完全搞错了。你的情况直接做TD单点即可，此任务输出中SCF Done能量为S0能量，E(TD-HF/TD-DFT)为S1能量。

，非常感谢老师的回复！如果从S1的td单点通过“SCF Done”提取出来的是基态能量，”total Energy, E(TD-HF/TD-DFT) = “这一行提取出来的是激发态能量。这样来看的话，我这个体系前一段的反应更有可能在激发态S1进行，后一段的反应在基态S0进行。

作者
Author: BinWang 时间: 2025-9-2 17:43

BinWang 发表于 2025-9-2 17:39
，非常感谢老师的回复！如果从S1的td单点通过“SCF Done”提取出来的是基 ...

作者
Author: wal 时间: 2025-9-2 17:44

BinWang 发表于 2025-9-2 17:39
，非常感谢老师的回复！如果从S1的td单点通过“SCF Done”提取出来的是基 ...

要判定这个的话你得搜一个S1势能面的TS，另外S1/S0圆锥交叉没法用一般TDDFT计算，如果你的反应过程中间发生了圆锥交叉，那你可能要换软件了

作者
Author: BinWang 时间: 2025-9-2 17:45

wal 发表于 2025-9-2 17:44
要判定这个的话你得搜一个S1势能面的TS，另外S1/S0圆锥交叉没法用一般TDDFT计算，如果你的反应过程中间发 ...

好的，非常感谢老师

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