计算化学公社

标题: 使用高斯计算激发态溶剂重组能数值小的问题 [打印本页]

作者
Author: 风起~ 时间: 2023-8-28 20:12
标题: 使用高斯计算激发态溶剂重组能数值小的问题
本帖最后由风起~ 于 2023-8-30 21:20 编辑

各位老师好，我想着利用高斯来计算荧光过程的溶剂重组能
根据《Gaussian中用TDDFT计算激发态和吸收、荧光、磷光光谱的方法》、《乱谈外重组能的计算》两篇帖子

以下所有计算均基于S1激发态的结构，使用态特定模型来计算的溶剂重组能
1.%chk=exts1.chk
#p pbe1pbe/def2svp em=gd3bj td=nstate=5 scrf=(externaliteration,read,solvent=toluene)

None

0 1
xyz coordinates

noneq=write

--link1--
%oldchk=exts1.chk
%chk=exts0.chk
#p scrf=(read,solvent=toluene) pbe1pbe/def2svp em=gd3bj geom=check guess=read

None

0 1

noneq=read

读取最后一个SCF Done得到基态处于非平衡溶剂下能量-1511.99750954

2.%chk=eq.chk
#p pbe1pbe/def2svp em=gd3bj scrf=(solvent=toluene)
读取最后一个SCF Done得到基态处于平衡溶剂下能量-1511.99757493

计算得到的溶剂重组能为(-1511.99750954+1511.99757493)*27.2114=0.00178eV

一：请问一下我以上使用态特定模型计算有问题么。为什么我算出来的值这么的小呀，这是不是哪里有问题

二：态特定模型校正得到的激发2.8493eV，而默认的平衡溶剂模型得到激发能为2.8896eV，两者之差约为0.04eV，这个数值某种意义上也是溶剂重组能，为什么会比0.00178eV大这么多呢，
是不是我在前面的一理解有误呢

作者
Author: sobereva 时间: 2023-8-29 23:41
如置顶的新社员必读贴、论坛首页的公告栏、版头的红色大字非常明确所示，求助帖必须在帖子标题明确体现出此帖内容是求助或提问，并清楚、准确反映出帖子具体内容，避免有任何歧义和含糊性，仔细看http://bbs.keinsci.com/thread-9348-1-1.html。我已把你的不恰当标题 “使用高斯计算激发态溶剂重组能” 改了，以后务必注意

作者
Author: sobereva 时间: 2023-8-29 23:48
跟特定态模型对比便知

第二个任务得到的是平衡溶剂下的能量

作者
Author: 风起~ 时间: 2023-8-30 09:28
本帖最后由风起~ 于 2023-8-30 21:22 编辑

sobereva 发表于 2023-8-29 23:48
跟特定态模型对比便知

第二个任务得到的是平衡溶剂下的能量

谢谢sob老师回复，原帖有几处不当，我重新编辑了一下帖子

我这个就是态特定模型下计算的溶剂重组能，不知道流程有问题么

作者
Author: sobereva 时间: 2023-8-31 00:29

风起~ 发表于 2023-8-30 09:28
谢谢sob老师回复，原帖有几处不当，我重新编辑了一下帖子

我这个就是态特定模型下计算的溶剂重组能， ...

没问题

作者
Author: 风起~ 时间: 2023-9-1 17:15
查了相关文献，激发态的非极性溶剂的外重组能就是非常小的，这个结果在数值上是正常的

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