计算化学公社

标题: 求助关于使用TrEsp进行激子耦合计算 [打印本页]

作者
Author:
yangjunfang    时间: 2025-12-1 15:21
标题: 求助关于使用TrEsp进行激子耦合计算
本帖最后由 yangjunfang 于 2025-12-1 15:22 编辑

老师好,我想请教一下关于舒勇TrEsp进行激子耦合计算,我的计算步骤如下:
先将二聚体的两个单体根据该关键词进行计算,#p td IOp(9/40=4) cam-b3lyp/6-311+g**,得到激发态,计算S0->S1的TDM。
然后将fchk输入multiwfn中,根据下面步骤产生TDM.fch
18->9->1->log文件->2->y->y产生两个单体的TDM.fch

然后将TDM.fch输入multiwfn产生每个单体的chg
7->12->5->3->1这样就产生了TrEsp电荷

然后将两个chg文件合并,继续载入multiwfn,通过步骤7->-2得到激子耦合结果如下图,这个很不合理,请问我是哪里有问题嘛



作者
Author:
sobereva    时间: 2025-12-2 06:57
舒勇 语义不明

chg文件里的坐标明显离谱,载入Multiwfn进主功能0一看便知
(, 下载次数 Times of downloads: 5)

记得恰当使用nosymm避免调整到标准朝向下
谈谈Gaussian中的对称性与nosymm关键词的使用
http://sobereva.com/297http://bbs.keinsci.com/thread-1355-1-1.html


作者
Author:
yangjunfang    时间: 2025-12-2 10:12
sobereva 发表于 2025-12-2 06:57
舒勇 语义不明

chg文件里的坐标明显离谱,载入Multiwfn进主功能0一看便知

老师,您好,这个坐标确实和初始的坐标不一致,然后我换了和初始坐标一致的坐标,计算的结果 Electrostatic interaction energy: 13466044.64 kJ/mol  3218464.17 kcal/mol,感觉也是不合理的
因为我这个需要考虑对称性,所以没使用nosymm
作者
Author:
wal    时间: 2025-12-2 16:22
你这两个单体互相连一点偏移都没有,优化出来的无虚频结构真是这样吗?可以上传一下结构优化的log么
作者
Author:
yangjunfang    时间: 2025-12-3 15:50
wal 发表于 2025-12-2 16:22
你这两个单体互相连一点偏移都没有,优化出来的无虚频结构真是这样吗?可以上传一下结构优化的log么

我这个分子没有优化
作者
Author:
wal    时间: 2025-12-3 16:14
本帖最后由 wal 于 2025-12-3 16:37 编辑
yangjunfang 发表于 2025-12-3 15:50
我这个分子没有优化

没优化那你算出什么结果都在意料之内 (, 下载次数 Times of downloads: 6)




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3